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Abstract 


The properties of several atmospheric models, which have been proposed for the purpose 
of numerical computation of prognostic maps, are discussed, and a new version of the 
“21/,-dimensional model” is given. It is pointed out that the “barotropic” model is unable 
to explain adequately the growth of wave perturbations, and that the advective model 
gives too strong instability for short waves, whereas the two-layer model and the contin- 
uous 21/,-dimensional model seem to describe almost correctly the behaviour of waves 


in a baroclinic current. 


I. Introduction 


Looking back at the results of the research 
on numerical weather forecasting carried out 
after the last war, we may say that a first stage 
of the research program is now finished, since 
it is definitely established that the two-dimen- 
sional, non-divergent model is able to explain 
a considerable part of the changes that we 
observe on the 500 mb maps. We have thus 
arrived at a serviceable first approximation, 
which enables us to carry on the research 
work with more confidence than before. 
The next stage, the advance from this bridge- 
head towards more accurate methods for 
numerical weather forecasting, has already 
begun. The goal is to find a system of prog- 
nostic equations which is realistic enough to 
give forecasts of practical value, and at the 
same time simple enough to permit computa- 
tions of such forecasts with a reasonable 
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amount of labour. Since this is a compromise, 
there may not be just one solution, suitable 
for all purposes. We may perhaps need several 
methods, of various degrees of accuracy and 
mathematical complexity; which one to use 
will depend on the nature of the forecast in 
question, as well as on the capacity of the 
available calculating equipment. 

The derivation of suitable prognostic equa- 
tions, which fulfil the conditions mentioned 
above, is one of the two main problems of 
numerical weather forecasting. We may 
hope to arrive at such equations by testing a 
variety of intelligent proposals. A discussion 
of some of the models that have been proposed, 
together with the author’s version of the “2 %- 
dimensional model”, is given below. The 
second main problem of numerical forecast- 
ing, to overcome the mathematical difficulties 
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involved in the process of numerical inte- 
gration, will not be discussed here. 


2. Quasi-geostrophic equations in 
three dimensions 


So far, most of the prognostic equations 
that have been considered have been based 
on the quasi-static and the quasi-geostrophic 
approximations. The quasi-static approxima- 
tion is valid for motions whose horizontal 
extent is large in comparison with their 
vertical extent, and is generally believed to be 
legitimate for all large-scale motions in the 
atmosphere. More questionable is the legiti- 
macy of the quasi-geostrophic approximation, 
which has been discussed by CHARNEY (1948). 
If we study solutions of linearized equations, 
we find that this approximation, together 
with the quasi-static approximation, elimi- 
nates all sound waves and ordinary gravita- 
tional-inertial waves whose perturbation veloc- 
ities are completely non-geostrophic. On 
the other hand, the same equations describe 
almost correctly the behaviour of waves 
with nearly geostrophic perturbation veloci- 
ties, such as the gravitational-inertial waves 
in a baroclinic current discussed by CHARNEY 
(1947) and Eapy (1949), as well as Rossby- 
waves. Since most large-scale disturbances 
connected with weather seem to belong to 
the group with almost geostrophic perturba- 
tion velocities, we have reason to hope that 
the quasi-geostrophic approximation will serve 
our purpose. On the other hand, we know 
from synoptic experience that large non- 
geostrophic winds sometimes exist in the upper 
troposphere; in particular, we must expect 
severe difficulties when applying the geo- 
strophic approximation to sharp upper ridges, 
whose importance for the dynamics of the 
flow has been stressed by J. BJERKNES (1951). 
We must therefore expect the quasi-geostroph- 
ic approximation to cause considerable errors 
in some cases. | 

The quasi-static approximation enables us 
to use pressure (p) instead of height (+) as a 
vertical coordinate, which is convenient for 
several reasons. This means that partial differ- 
entiation with respect to horizontal coordi- 
nates (x, y) and time (f) will be performed at 
constant pressure, such that 0/dx and olay 
denote differentiation along an isobaric sur- 
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face, and 9/9t means rate of change as measured 
by an observer which follows an isobaric 
surface in its vertical motion. 

It is convenient to use vector notation v 
for the horizontal wind velocity (with com- 
ponents u and v), and to introduce a hori- 
zontal del-operator V7 (with components 
9/9x and d/dy) to express differentiation with 
respect to horizontal coordinates (at constant 
p). Let further k denote a vertical unit vector, 
pointing upward; w individual rate of change 
of pressure, « specific volume, @ potential 
temperature, f Coriolis parameter and g 
acceleration of gravity. 

The hydrostatic equation may then be 
written 


(1) OZ a 


The equation of continuity assumes the simple 
form 


a) 
(2) we 


Neglecting frictional forces, we derive from 
the horizontal equation of motion the vor- 
ticity equation 


dË 3 0 
6) F+vVU+D+0R+ 


d 
+k-Vox mt U+O ve 
where 
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is the “isobaric” vorticity, which differs from 
the vertical vorticity component in the differ- 
entiations being taken along isobaric instead 
of level surfaces. 

Elimination of \7- v between (2) and (3) 
gives the vorticity equation in the form 


a a 

(s) 3 7 il ah hao a 
IV Ow 
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Since the amount of heating and cooling 
is not known, we are at the present time 
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unable to take this effect into account, and 
will have to assume that the potential tem- 
perature is conserved. The equation expressing 
this assumption is 


ua Ka o 00 
(6) 


cites DIR pq 9 


since differentiation at constant pressure gives 
oln O/dt= In «/dtand 7 nO= V7 In«. Making 


use of the hydrostatic equation (1) in the first 


term, this equation may also be written 
(7) Ce __ ao 90 _ 
N pag gO a 


The quasi-geostrophic method consists in 
substituting the geostrophic wind for the true 
wind in the terms which are not critically 
affected by this substitution. Thus the geo- 
strophic wind approximation cannot be used 
in the divergence term, or in the expression 
for the work of the horizontal pressure force 
(if the energy equations are considered). 
The simplest and most logical way to intro- 
duce the geostrophic approximation is given 
by CHARNEY (1948). His method consists in 
substituting the geostrophic wind for the true 
wind everywhere in the thermodynamic 
energy equation and in the vorticity equation, 
after having eliminated the horizontal diver- 
gence by means of the equation of continuity. 


From the geostrophic wind formula 
(8) ve Ek x iz 
a8 


we find the following approximate formula 
for the geostrophic vorticity 


(9) Fey ho aes =). 


From these equations and the hydrostatic 
equation (1) we find 
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Using Charney’s method, we now introduce 
the geostrophic wind approximation into 
equations (s) and (7), and obtain 


dz f(f+%,)0w I 


er Gen Scale ah a ze 
= —k-VzxV(f+Q,) 
à Oz a JO I 


(13) Viet TO ee ap nes Witte 


This may be considered as a system of linear 
differential equations for the two unknown 
functions Jz/dt and w, whose coefficients can 
be determined if the field z (x, y, p) of the 
contour height is known at the instant con- 
sidered. When the absolute vorticity and the 
static stability are both positive, which corre- 
sponds to normal conditions in the atmos- 
phere, the system is of the elliptic type. 
Together with the boundary conditions at 
the top and the bottom of the atmosphere 


@ S060 at p = ©; 
(14) = Fang So ake ground, assum- 
fe Lie ed to be level 


(p = po © 1000 mb) 


and suitable conditions for the lateral bound- 
aries, the system (12, 13) determines dz/dt 
and w throughout the region considered. The 
solution may be worked out numerically, 
and may be used to extrapolate the field 
Z(x, y, p) over a small interval of time. Iterating 
this procedure, we obtain an integration in 
small steps of time. This method, based on 
essentially the same system of equations, has 
been discussed previously by the author 
(ELIASSEN, 1949). 

In order to solve the system (12, 13) nu- 
merically, we may eliminate one of the 
unknown functions, whence we obtain a 
Poisson-type differential equation in three 
dimensions for the other, which may be solved 
e.g. by relaxation. Such an equation for the 
height (or pressure) tendency has been discuss- 
ed by CHARNEY (1949). 

It is interesting to note that if the third 
term on the left of (12) is neglected (which is 
often done), then the system (12), (13) is 
equivalent to the equilibrium equations for a 
loaded elastic system consisting of a number 
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of horizontal elastic membranes exposed to 
weights of variable magnitude, and connected 
with vertical springs. Interpreting dz/dt as 
the vertical displacement of the membranes 
and © as the tension in the vertical springs, we 
see that the equation expressing equilibrium 
between the vertical forces has the form (12), 
whereas (13) expresses a linear stress-strain 
relationship for the vertical springs. The 
boundary conditions (14) also have obvious 
interpretations in this elastic model. Since we 
have an immediate feeling of how this model 
will behave under the influence of a given 
distribution of weights, and for given values 
of the constant term in Hooke’s law (right- 
hand side of (13)), it may serve to give us an 
idea of what the solutions of our meteorologi- 
cal problem will look like. 

Computations of forecasts based on the 
three-dimensional, quasi-geostrophic equations 
involve a large amount of computational 
work. To the author’s knowledge, no such 
forecasts have yet been computed. So far, 
the tests have been based on simplified 
equations, which reduce the problem to one 
in two dimensions in space. This is attained by 
making certain assumptions, based on empir- 
ical knowledge, of the variation with height 
of the field of motion. 


3. The barotropic model 


The simplest model of this kind is the so- 
called barotropic model. The suggestion that 
the large-scale motion of the atmosphere may 
be considered as two-dimensional and non- 
divergent was first put forward by RossBy 
(1939). CHARNEY (1949) has given a more 
precise formulation of the physical basis for 
this hypothesis. He uses the vorticity equation 
(5) in the simplified form 


x L ~ 00 
CREATED 


=. 


with the simplified boundary conditions 
(16) 


Furthermore, he makes use of the observed 
fact that as far as the large-scale motion is 
concerned, the magnitude of the wind velocity 
tends to change with height in a typical 
fashion, whereas the streamline pattern stays 


w@=0 at p=o and at p= py. 
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about the same from level to level. This 
assumption makes it possible to carry out an 
integration of (15) with respect to p from o 
to po, thus eliminating w. The resulting 
equation expresses that at a certain inter- 
mediate level, which is found to be in the 
vicinity of 500 millibars, vorticity changes 
are due to advection only: 


(17) 2 - —v-V (f+ ¢) at p » 500 millibars. 


The motion at this level may therefore be 
considered as horizontal and non-divergent. 
When the geostrophic assumption is intro- 
duced, a two-dimensional Poisson equation 
for the height tendency is obtained, which 
may be used as a prognostic equation’. 

Several numerical tests have been carried 
out by comparing the changes computed from 
this equation with the observed ones (CHARNEY, 
ELIASSEN, 1949; CHARNEY, FJORTOFT, VON 
NEUMANN, 1950; BOLIN, CHARNEY, 1951; 
STAFF MEMBERS, UNIVERSITY OF STOCKHOLM, 
1952). The results are very promising indeed; 
the discrepancies, although important, are 
nothing more than we would have to expect 
from such a crude model. 

The properties of horizontal, two-dimen- 
sional flow, which is governed by equation (17), 
has been subject to many theoretical studies. 
RossBy (1939) investigated horizontal, non- 
divergent wave perturbations on a uniform 
zonal current, and found that these waves 
move with the speed 


(18) c= De 


where U is the velocity of the current, f= 
dfidy, and k the wave number. These waves 
are stable for all wave lengths. CRAIG (1945) 
and Hortanp (1950) have found solutions 
representing two-dimensional non-divergent 
waves with finite amplitude on a spherical 
earth. These theoretical motions have a great 
resemblance with the observed wave patterns 
on the 500 mb maps. 

The most apparent defect of the barotropic 
model seems to be connected with the growth 


1 Concerning the details, the reader is referred to 
CHARNEY’S paper (1949) and also to CHARNEY and Ertas- 
SEN (1949), CHARNEY, FJORTOFT, VON NEUMANN (1950). 
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of wave perturbations. Unstable waves do 
exist in the barotropic model, but the only 
source of energy for the growing perturba- 
tions is the kinetic energy of the basic flow. 
Investigations made by CHARNEY (1947), 
EADY (1949) and Fyortorr (1950), of wave 
perturbations in a baroclinic, zonal current, 
based on the three-dimensional equations, 
have shown the existence of unstable waves 
that grow at the expense of the potential 
energy of the system; and there is rather 
strong evidence that most waves in the wester- 
lies develop as a result of instability of this 
kind. Such unstable waves do not exist in the 
barotropic model, since transformations of 
potential into kinetic energy are excluded a 
priori in this model. We must conclude that 
the barotropic model is not capable of an 
adequate explanation of the growth of wave 
perturbations, and we must therefore look 
for more realistic models. 


4. The advective model 


It is possible to take into account important 
‚effects related to transformation of potential 
into kinetic energy, i.e. to the baroclinicity, 
without using the complicated three-dimen- 
sional equations. So called 2 %-dimensional 
models have been designed, which in a crude 
manner take into account the three-dimen- 
sional structure of the atmosphere and the 
vertical motion, although they are still two- 
dimensional, mathematically speaking. In these 
models, two horizontal fields (e.g. one 
contour field and one temperature field) are 
being used to define the state and motion of 
the atmosphere. A better name would there- 
fore perhaps be “two-parametric models”. 

The simplest model of this type is the 
advective model. It is based on the assumption 
that temperature changes are caused by hori- 
zontal advection only. Moreover, one assumes 
that the thermal wind has a similar variation 
in magnitude with height in all places, whereas 
its direction is independent of height. In other 
words, one assumes that the wind field can 
be represented by the formula 


(19) v(x, 7, p,1) = V(x, y,t) + A(p) vr(x, 7.0) 


where the bar denotes the mean value with 
respect to pressure, A (p) represents the varia- 
tion of thermal wind with height, and vr is the 
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thermal wind measured between two selected 
isobaric levels. Introducing the thickness h 
between these isobaric levels, and the mean 
contour height z, the geostrophic assumption 
gives 


Vr =F kx Th, 


and the corresponding vorticities 


en) 


Kix v2 


Bae run 


J ie 


When these expressions are introduced into 
the approximate vorticity equation (15), and 
the mean value with respect to pressure is 
taken, one obtains 


(2) = 
dt 

where the last term gives the contribution of 

the thermal wind field to the mean vorticity 


change. The advective hypothesis gives fi- 
nally 


(23) 


oh a 
Een Vie 


The two last equations are seen to contain the 
two horizontal fields z and h only, and may be 
integrated numerically in steps of time to 
give the change of these fields. 

In a form suitable for numerical integration, 
the advective model is described in a paper 
by Carney, Fyortorr and von NEUMANN 
(1950), (see also Fyortorr, 1951, and CHAR- 
NEY, 1951). The assumption that tempera- 
ture changes are due to horizontal advec- 
tion only has actually been used earlier by 
RossBy (1942) and by J. BJERKNES and Horm- 
BOE (1944) in the study of long waves in a 
baroclinic zonal current, although without 
the derivation of the corresponding formula 
for the wave velocity. If we linearize (22) and 
(23) for the case of small perturbations on a 
baroclinic zonal current without horizontal 
shear, we find for sine waves of wave number 
k and of infinite lateral extent the velocity of 
propagation 


C= a - VE)-F Ur 
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where U is the mean wind speed and Ur 
the thermal wind of the current. This formula 
has been derived by Fyorrorr (1950) and by 
SUTCLIFFE (1951). À similar formula is given 
by BERSON (1950). 

Equation (24) gives stable or unstable 
waves, according as the wave length is greater 
or smaller than a certain critical wave length, 
the magnitude of which increases with the 
square root of the thermal wind. This is in 
agreement with results derived by CHARNEY 
(1947) from the three-dimensional equations. 
Thus the advective model seems to be more 
realistic than the barotropic model, which in 
this case gives stable waves only. For wave 
lengths considerably larger than the critical 
wave length, one of the values (24) of the 
wave speed is in approximate agreement 
with the Rossby formula (18), whereas the 
other is considerably larger. This second 
value, however, corresponds to waves with a 
nodal surface in the middle of the atmosphere, 
roughly at the level where the barotropic 
model applies. These waves are therefore 
eliminated a priori in the barotropic model. 
We may conclude that the barotropic model 
is in agreement with the advective model 
as far as waves considerably longer than the 
critical wave length are concerned. For waves 
comparable with, or shorter than the critical 
wave length, the results derived from the 
two models differ appreciably (compare 
SUTCLIFFE 1951). : 

The unstable waves derived from the advec- 
tive model have the property that their 
growth-rate increases indefinitely as the wave 
length approaches zero. This “ultra-violet 
catastrophe” has the unagreeable consequence 
that we are unable to solve the linearized 
initial value problem by means of Fourier 
analysis; most likely, this problem has no 
meaning. We may find a solution by dropping 
the shortest wave components, which is 
automatically done if the continuous functions 
are replaced by their values in a grid of points. 
However, such a procedure will not eliminate 
the difficulty, because the solution will depend 
critically on the wavelength of the shortest 
waves kept in the computation, i.e. on the 
grid size. We must expect to find the same 
difficulty also if we apply the equations of 
the advective model in non-linear form. 

The extreme instability for very short 
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waves, obtained from the advective model, is 
not supported by the theory based on the 
three-dimensional equations. This instability 
is just an error, introduced by the advective 
assumption, which is not legitimate for short 
waves. As we shall see, the error is quite un- 
necessary, and the advective model is there- 
fore not a particularly good model to use for 
the purpose of numerical forecasting. 

We shall now proceed to discuss some two- 
parametric models where also the vertical 
transport of potential temperature 1s taken 
into account. 


5. Phillips’ two-layer model 


Puitups (1951) has designed a model con- 
sisting of two homogeneous and incompress- 
ible fluid layers separated by a quasi-horizontal 
interface and bounded by two rigid, horizon- 
tal planes. The density of the upper layer is 
somewhat smaller than that of the lower, 
so that the system is statically stable. The 
system is baroclinic if the isobaric surfaces 
intersect the interface; and since each layer is 
autobarotropic, the baroclinicity is concen- 
trated in the interface. The interface is per- 
mitted to slope and to move up or down, 
corresponding to vertical motion and hori- 
zontal divergence in each layer. 

From the hydrostatic relation and the dynam- 
ical boundary condition at the interface one 
derives a linear relationship between the 
contour heights z(x, y, f) and 2’ (x, y, t) of 
the two layers and the height h (x, y, t) of 
the interface. The state of the system is thus 
defined by two two-dimensional fields, e. g. 
the fields of z and h. Application of the 
geostrophic potential vorticity equation to 
the two layers gives for these two fields 
differential equations of the form 


oh oh 

(as) az D (x, y) ae (x, y) 
Oz 

(26)my RES = r(x, }) 


where the coefficient functions a and q can 
be computed from knowledge of the fields 
h and z, and r can be computed from know- 
ledge of h, z and dh/ot.x Together with 


! For further details see PmırLip’s paper (1951). 
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suitable lateral boundary conditions, these 
equations may be used to determine dh/dt and 
dz/dt when z and h are known at the instant 
considered. The tendencies thus computed 
may be used to extrapolate z and h through a 
finite interval of time; and iteration of this 
procedure gives a numerical integration in 
steps of time. It is important to note that 
although the motion of the model is three- 
dimensional, the forecast problem is mathe- 
matically two-dimensional in space, since 
the vertical coordinate (p) does not occur as an 
independent variable. 

Phillips has applied his model to study the 
behaviour of perturbations on a uniform 
zonal current with different wind speeds in 
the two layers. For long waves, his results are 
in agreement with Charney’s results and with 
the advective model: waves longer than a 
critical wave length are stable, and shorter 
waves unstable. However, as the wave length 
decreases, the growth rate of the unstable 
waves does not increase indefinitely, as in the 
advective model, but instead increases to a 
maximum and then falls down to zero at 
another critical wave length. Still shorter 
waves are stable in Phillips’ model. This part 
of the frequency curve is in agreement with 
the results obtained by Eapy (1949) from the 
three-dimensional perturbation equations for 
small wave perturbations on a baroclinic 
current. The whole frequency curve for the 
two-layer model corresponds closely to the 
frequency curve for a continuous medium 
derived by Fyorrorr (1950) from energy 
considerations. This is a strong indication that 
the two-layer model in many respects will 
behave like the three-dimensional atmos- 

here. 
Phillips has also applied his model to an 
actual weather situation and computed tend- 
encies and vertical velocities which look 
very promising. 

A disadvantage of the two-layer model is its 
geometrical dissimilarity with the real at- 
mosphere, which makes it difficult to define a 
proper correspondence between the states of 
the model and the states of the atmosphere. 
Phillips establishes such a correspondence by 
requiring the frequency curve, derived from 
the model for waves on a zonal current, to be 
quantitatively correct in its main properties. 
But even with the best possible choice of 
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correspondence, one would expect the model 
to reproduce atmospheric motion with a 
certain amount of distortion. For the purpose 
of numerical forecasting, a. more life-like 
model would therefore probably give better 
results. 


6. Continuous two-parametric model 
with vertical stability 


The continuous two-parametric model with 
static stability is a natural generalization of 
the barotropic and advective models, and it 
is therefore not surprising that several re- 
search workers seem to have arrived at this 
model more or less independently. Thus 
Eady has proposed a very simple and elegant 
version of this model which he presents in an 
article in this issue of Tellus (cf pp 157—167). 

In the author’s version of the 2 %-dimension- 
al model, described below, it has been attempt- 
ed to attain the closest possible correspond- 
ence with atmospheric conditions, in order 
to minimize errors due to geometrical dissim- 
ilarity between the atmosphere and the 
model. The author’s model is somewhat more 
complicated than Eady’s, but does not differ 
much from Eady’s model in its main proper- 
ties. 

In the two-parametric model, the tempera- 
ture distribution along the vertical is represent- 
ed by one parameter. In other words, the 
possible ascent curves of our model are restric- 
ted to a one-dimensional mannifold of curves. 
Moreover, it will be convenient to have the 
temperature (T) as a linear function of the 
parameter defining the particular ascent curve. 
We will therefore make the following assump- 
tion for the vertical temperature distribution: 


(27) TT (p)=Tn (p) +t (p) (h—h). 


Here h is the parameter which defines the 
temperature curve, and h, a constant included 
for the sake of convenience. The function 
Tn (p) is chosen as a “normal” ascent curve 
(fig. 1), and t(p) is chosen positive in the 
troposphere and negative in the stratosphere, 
as shown in fig. 2. Formula (27) then gives, 
for different values of h, a family of ascent 
curves as shown in fig. 3. These curves may be 
thought of as a set of homologes for different 
air masses; they are in accordance with real 
conditions in so far as a warm troposphere 
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Fig. 1. ‘“‘Normal’’ ascent curve Th: 


goes together with a cold stratosphere and 
vice versa. Since the result is unchanged if r is 
reduced and at the same time h and h, in- 
creased by the same factor, it follows that we 
may always choose h to represent the thick- 
ness of a prescribed isobaric layer, e. g. from 
1,000 to 500 mb; h„ will then be the corre- 
sponding thickness for the “normal” tem- 
perature curve T,. Assuming h to be a func- 
tion of x, y and f, formula (27) will represent a 
three-dimensional temperature field. 

The corresponding field of specific volume is 


(28) 


er 
a (0 77.8) = RE + 


HER u [h (0% ys t) — hu]; 


(R, gas constant), and the contour field 


(29) z(x,y, p, t) = 2 (x, y, 9) + H(p) + 
Ap ER 


Fig. 3. Ascent curves for different values of h. 


ARNT ELIASSEN 


D NE RE oe 2 A 

Anns aan) 

40 80x10 54m 
7 


Fig. 2. Temperature deviation T per unit thickness. 


Here the bar denotes as before the mean 
value with respect to pressure, and the func- 
tions H(p) and A(p) are defined by 


(30) Nele), H=o, 
oP 
(31) are SL), A=o. 
SP 
The three-dimensional field z (x, y, p, f) thus 


depends on the two horizontal fields of z and h. 


From (29) we find the field of geostrophic 
wind 


(31) v y,p,t) =V (x, 9,4) + A (p) vr (x,y, 4), 


and the corresponding geostrophic vorticity 
(32) (pp) =E (x 7,0 + A(p) or (x, y, 0), 
where V, vr, &, and Cy are defined by (20) 
and (21). Fig. 4 shows the function A(p). 
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Fig. 4. The function A(p). 
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We will also have to make a certain assump- 
tion about the distribution with height of 
vertical motion, and put 


(33) © (x, y, p, )=B(p)-w (x y, #) 


where B (p) is chosen as a function of the 
type shown in fig. 5. The approximate 
boundary conditions (16) are automatically 
fulfilled when B is chosen so that B (0) = B (po) 
= 0. 

We now insert the formulae (28), (31) and 
(33) into the thermodynamic equation (6). 
Approximating «, when undifferentiated, as 
well as the static stability factor, by their 
normal values, we obtain 


3 dh =! 
(34) + (+ ¥-Vh) + 


where ©, is the potential temperature corre- 
sponding to the “normal” temperature T,. 
As a rule, we cannot attain this equation to be 
exactly fulfilled at all levels; a reasonable 
requirement is, however, that the mean value 
of its left hand side shall vanish. This gives 


B dO, 
O, dp 


w=oO 


> 


dh = 
CC) 


where we have put 


T HABidOwWebiat 
0) 9=(r) LG) - 
The constant y can be evaluated when the 
functions T,, tT and B have been chosen. 


The vorticity equation shall be used in the 
form 


a ac 
FERIEN CE 
Ow 
4 Rz , 
bite ha O 


where only the term involving Vo of 
the “exact” equation (5) has been dropped. 
It is even possible to take this term into 
account, but that would complicate the compu- 
tations considerably. Introduction of (31), 


(32) and (33) into (37) gives 


B(p) 
Fig. 5. The function B (p). 
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As in the case of the thermodynamic equa- 
tion (34), we cannot hope to have the vorticity 
equation exactly satisfied at all levels. Our 
model is thus inconsistent, a fact which, 
however, will not disturb a theoretical meteor- 
ologist. We have the freedom to impose 
two conditions on the left-hand side of the 
vorticity equation, because we need two 
prognostic equations for the two horizontal 
fields z and h. What we have to do is to 
choose these conditions so as to minimize the 
error in the vorticity equation. 

For fixed values of x, y and t, the left-hand 
side of (38) is a function of p, which will be 
denoted by M(p). To this function we de- 
fine a linear regression function L(p), such 


that (L—M)? = a minimum. We now require 
L(p) =o for all values of p. This gives us 
the two wanted conditions 


(39) M =o and pM—p-M=o. 


When we here insert for M the left-hand side 
of (38), carry out the details of the calculation 
and eliminate w by means of (35), we obtain 
the two prognostic equations 
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Equations (40) and (41) represent, when 
the geostrophic formulae (20) and (21) are 
inserted, two linear differential equations for 
the tendencies 97/9t and 9h/9t, of the same 
type as the prognostic equations (25), (26) 
obtained for the two-layer model. The sys- 
tem (40), (41) may be integrated numerically 
in steps of time. At each step, we will first 
have to determine dh/dt from (41); having 
done that, 2z/9t may be obtained from (40). 
Both equations may be solved by relaxation. 
If desired, the contour height at any partic- 
ular level (e.g. soo mb) may be used instead 
of the mean height z, the relation between 
them being given by (29). This has the ad- 
vantage that the two dependent variables h 
and z (500) can then be compared directly 
with the synoptic maps. 

Each step will involve the evaluation of 
the fields of five advection terms, which have 
the form of Jacobian determinants. Their 


number may, however, be reduced to three. 
Let 


(43) h=(v+ Avr: V(fHT+Aln), 


Ji = (WE ANT) 7 Ft Ae 


denote the vorticity advection at two levels p, 
and p,, which are chosen such that A, = 
A(p,) and A, = A(p,) are the roots of the 
equation 


(44) 
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Then we have 
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(as) VG 2) aver in riggs 
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which we may introduce into (40) and (41). 

The coefficients of the prognostic equations 
(40), (41) may be computed from (42) and 
(36), when the functions 7,, A and B have 
been chosen. If we adopt the values given in 
figs. 1, 4 and 5, which are based on esti- 
mates made by the author’s colleague, Mr. 
Smebye, we find 


a = 0.65, b 000 


(47) g=.0.89 - 102 10,7 — 0.70 10-000 
Ar =.0.85. 


However, there is no unique way of defining 
the functions A and B, and the above values 
of the coefficients may therefore not be the 
best possible ones. 

To the author’s knowledge, all two-para- 
metric models designed so far give the same 
prognostic equations (40), (41), only with 
different numerical values of the coefficients. 
This applies to the model proposed by Eady, 
as well as to Phillips’ two-layer model, when 
suitably interpreted. For the purpose of nu- 
merical forecasting, it is therefore not necessary 
to pay too much attention to the details of 
the structure of the model; the best approach 
is perhaps to determine the coefficients of 
(40), (41) empirically so as to get the least 
error in the computed forecast. 


7. Waves in a baroclinic zonal current. 


We shall finally investigate on the basis 
of the continuous two-parametric model the 
speed of propagation and growth-rate of 
wave perturbations in a baroclinic zonal 
current. Let the wind velocity of the basic 


current at any level be given by U+ A(p) Ur, 


where U and Ur are constants. For sim- 
plicity, we shall let our coordinate system 
move with the mean wind, so that we may 


assume U = 0. Assuming the perturbations = 
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‘and h of the mean height and thickness to be 
functions of x and t only, we find from (40) 
and (41) the perturbation equations: 


0° = Oz 03 h 
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These equations may be satisfied by a trigono- 
metric wave solution of the form: Zz, h ~ exp 


i (kx+vt), provided 


IR 2 
vm + Voi—vi 


v= 


(st) 
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2 2 
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(53) 


When » > vi we obtain stable waves 
with the speed of propagation c = — v/k. 
When 1% < v2, the waves are unstable; their 
speed of propagation is c= — B/k® + vy/k, 
and their growth-rate »; = y,— vr. For the 
values of a and b given in (47), the speed of 
propagation and the growth-rate are shown 
in fig. 6 as functions of the wave length L 
and the thermal wind Ur. The quantities L, 
y;, c and Ur are in the diagrams represented 
as the non-dimensional quantities ky L/2z, 
y;/¥9, ¢/Uy, and Ur/U,, respectively, where 


ay p B 


i ata 
0 be 0 ko? 


This has the advantage that the diagrams are 
valid for all latitudes and for all values of the 
coefficient r in (50), which involves the static 
stability. In latitude 45°, and with the value 
of r given in (47), we find 
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Fig. 6. Upper part: Growth rate v,/v, of unstable waves 


as a function of wave-length K,L/2n, for different 
values of thermal wind UT/U,. 

Lower part: Speed of propagation c/U, of stable and 
unstable waves as a function of wave-length K,L/22, 
for different values of thermal wind UT/U,. 
Dashed curve: speed of propagation of Rossby-waves 
(U 0). 
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We notice that the waves are unstable in an 
intermediate band of wave lengths, in agree- 
ment with the findings of Fyorrorr (1950) 
and Pauzzrs (1951). The growth-rate is a 
maximum when k,L/27 & 1.5, which with 
the numerical values given above corresponds 
to a wave length of about 3,300 km. The 
wave length of maximum instability is seen 
to increase in proportion to the square root 
of the gtatic stability, and will also vary in 
inverse proportion to the Coriolis parameter; 
this is in good agreement with synoptic ex- 
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perience. Thus, for less stable air, or in higher 
latitudes, we will obtain considerably smaller 
values of the wave length of maximum in- 
stability. The values of the growth-rate and 
the speed of propagation obtained from the 
diagrams seem to be reasonable. 

Thus the results obtained from the two- 
parametric model are plausible, as far as the 
behaviour of wave perturbations are concern- 
ed. This gives us some reason to hope that 
this model will give good results also when 
the equations (40), (41) are applied in non- 
linear form. Should these equations prove 
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insufficient, there is still the possibility of 
increasing the number of parameters in the 
model. 
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Note on Weather Computing and the so-called 


24-dimensional Model 


By E. T. EADY 
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(Manuscript received 21 July 1952) 


Abstract 


The problem of setting up a system of equations by means of which the future large- 
scale flow pattern might be computed, using an electronic calculating machine, in a rea- 
sonable time and with maximum accuracy and efficiency, is discussed. Equations deter- 
mining the development of a two-parameter, two-dimensional ‘model’, capable of repre- 
senting the major features of motion of a baroclinic atmosphere, are derived. Tests of 
the model are made on some representative problems, the accurate three-dimensional 


solutions of which are known. 


Now that actual computation of certain 
aspects of the weather forecast has become a 
practical possibility, interest in certain simpli- 
fied “models” of atmospheric motion has 
increased. It seems at first paradoxical that 
this should be so. For although one necessary 
prerequisite for weather computing, an ade- 
quate network of upper-air data, has only 
recently become available, the main cause of 
optimism regarding the possibility of actually 
carrying out the necessary calculations in the 
time available has been the development of 
large-memory high-speed electronic com- 
puting machines. These machines work so 
very much faster than human computers that 
it is natural to suppose that at last we can 
forget about crude “models” and compute 
changes in the “actual” atmosphere. Practical 
experience in attempting to design a computa- 
tion scheme is disillusioning. It is not merely 
that certain approximations (equivalent to 
replacing the “actual” atmosphere by a “mod- 


el” but resulting in very little change in 
those aspects of the motion related to weather) 
are highly desirable from the point of view of 
computational simplicity, speed and stability. 
The main difficulty, which appears to result 
as the combined effect of 4 dimensions (3 
space and 1 time) and non-linear equations, 
is the very rapid increase in computation time 
with increase in representation of “detail”. 

It is worth while to attempt to express the 
problem in general terms. At a given instant 
the weather situation may be regarded as re- 
presented by certain values of a number of 
parameters which may, but need not, be the 
values of pressure, temperature etc. at points 
on a 3-dimensional grid. Just as the definition 
on a television screen increases with increase 
in the number of spots so the accuracy with 
which the weather situation may be represented 
increases with increase in the number of 
parameters. The equations and boundary 
conditions of motion together with the ther- 
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mal and continuity equations determine the 
rate of change of these parameters and there- 
fore the forecast situation, but since the 
equations involve details either not represented 
or inaccurately represented by the parame- 
ters some approximations are inevitable. The 
forecast will therefore contain not only errors 
of detail due to the use of a finite number of 
parameters but these parameters themselves 
will necessarily be to some extent in error. 
It is important to recognize these two kinds of 
error which are unavoidable even when the 
computations are made with complete accu- 
racy. To distinguish them from a third type 
of error to be discussed later they will be 
referred to as “physical” or “model” errors. 
They arise simply and solely because we 
cannot include all the accurate data. Whether 
or not such errors are serious depends partly 
on the number of parameters, partly on the 
manner in which they are chosen. If the param- 
eters are judiciously chosen it may be pos- 
sible to represent with sufficient accuracy both 
the features in which we are most interested 
and those primarily responsible for changes in 
them, with only a comparatively small num- 
ber of parameters. The problem is to discover 
the most efficient representation, that which 
enables us to calculate the important features 
as directly as possible, short-circuiting irrel- 
evant detail. Ideally any significant variation 
in the initial values of any one or more of the 
parameters should correspond to a significant 
difference in the forecast. An advantage of 
looking at the problem from this point of 
view is that it suggests the right sense of pro- 
portion. Clearly there is no point in computing 
with more independent well-defined param- 
eters than are determined by the initial 
data. On the other hand the forecast will 
contain the same number of parameters and 
therefore the same amount of detail: it will 
be as good as one obtained by more elaborate, 
time-wasting methods. 

This abstract, mathematical approach gives 
us the general idea and is a good practical 
guide in the later stages of design of a com- 
puting scheme. In the early stages, however, 
the guidance is too vague and a different, 
physical approach which, though basically 
equivalent, is more specific, is preferable. 
Let us commence with an example. Ibis ais 
well known that if we are interested in motion 


JB, Wb, Je eID) NC 


on a grand scale we may usefully replace the 
“actual” 3-dimensional atmosphere by a 2- 
dimensional “model”, the “barotropic model”. 
There is a more or less close relation between 
the horizontal component of the grand-scale 
motion at about 500 mb in the atmosphere 
and the motion computed for the model. 
Here the model is a physically possible hydro- 
dynamic system different from the atmosphere 
but behaving, in many important respects, 
in a similar manner. We might have set up 
this model directly from physical reasoning. 
On the other hand if we had integrated the 
vorticity equation along the vertical and then 
neglected certain terms we should have ob- 
tained the same final equations and it might 
be regarded as an accident that these equations 
correspond to a simple hydrodynamical sys- 
tem. From the point of view of our earlier 
abstract approach we should have suppressed 
all the parameters representing variation in 
the vertical. Setting up a model is equivalent 
to making certain approximations and since 
all computers have to make approximations 
there is no inherent defect in the use of mod- 
el. On the contrary, if approximations 
correspond to setting up a model we can be 
sure that the former are at least self-consistent. 
It does not follow, however, that the most 
suitable approximations must correspond to a 
physically possible system. The representation 
to be described does not accurately correspond 
to any physical model though it does approxi- 
mately do so. If the word “model” had not 
already been used rather vaguely it might be 
preferable to use some other word but since 
one advantage of the representation is the 
simplicity of physical interpretation and since 
the epithet “2 %-dimensional” is already pic- 
turesque it is convenient to call this represen- 
tation a “model”. It represents an improve- 
ment on the barotropic model in so far as it 
contains a very crude representation of varia- 
tion in the vertical, so crude as to be considered, 
playfully, as worth only half a dimension! 
Such a crude representation, with only 
two parameters along each vertical, may be 
given some preliminary justification. In the 
first place the barotropic model, with only 
one parameter along the vertical, has already 
had some success. Theoretical reasoning also 
suggests that we may have a good sense of 
proportion if we use many more parameters 
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to describe horizontal variation than vertical 
variation. For, making much less stringent 
approxinations, we may set up the equation 
for a 3-dimensional model in the form: (Eapy, 
1949, p. 47) 


d | 2p BER ae 1: don: 
dla op © a(Vbz)2 | dy ax’ 
| dp 
5 [2 35 
Vb = es ee (1) 


Here p is the pressure deviation from a standard 
value, x and y are horizontal and z is the 
a A à I 00 
vertical co-ordinate. (or — —, where @ is 
dz IZ 


potential temperature) is the static stability, 


g the acceleration of gravity, f the coriolis 
factor. Thus Vb is a pure number. A typical 


value in middle latitudes would be Vb = roo. 
Now although (1) is inaccurate for disturb- 
ances which are too intense, the form of 


the equation strongly suggests that Yb is a 
scale factor, representing the ratio of dynam- 
ically equivalent distances in the horizontal 
and vertical respectively. Thus the distance 
from ground to tropopause (say 10 km) is 
dynamically equivalent to a horizontal dis- 
tance of about 1,000 km. If therefore we 
devote 2 parameters to representation along 
the vertical we should use 1 parameter in 
about soo km in each horizontal direction. 
Since we need to consider the field over many 
thousands of kilometers in the horizontal 
it is clear that most of our parameters should 
describe horizontal variation. Unless the data 
are sufficiently detailed and accurate to deter- 
mine with some accuracy more than one 
parameter per 500 km x 5oo km square we 
are not justified in attempting to use a more 
accurate model. Although we cannot say 
without more investigation just how many 
parameters at a given (upper) level are well- 
determined by the data this number cannot be 
greater than the number of radio-sondes. It 
may, owing to inaccuracy and bad spacing of 
the sondes be much less. Hence it would 
appear that models of the kind to be described 
are not much, if at all, cruder than the best 
justified by the available data. Any attempt to 
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represent much more detail would be a waste 
of time in weather computing. (This is not 
true of course for calculations designed to give 
theoretical insight.) 

Instead of looking at the problem from the 
point of view of the amount of detail justified 
by the initial data we may look at it from the 
point of view of the computation time. 
Suppose we have an? parameters where a=1 
for the barotropic model and a=2 for the 
present model. # is the number of parameters 
representing variation in each horizontal 
direction. Now with linear equations it is 
frequently possible to choose the parameters 
in such a manner that the variation of each is 
independent of all the others. With non- 
linear equations it is, however, not usually 
possible to separate the variations in this 
manner and in general the variation of each 
parameter depends on all the other parame- 
ters. It is evident that for accurate calculations 
the number of arithemetical operations will 
increase much faster than #2. It looks as if 
the number of operations should vary as n4 
(as it does when the relevant differential 
equations are solved by means of a Green’s 
function) but it is possible with both kinds 
of model so to arrange the work that the 
number of operations varies only as n?. Nev- 
ertheless the time required for the computa- 
tion of time-variations, and therefore for each 
small time-step in preparing the forecast, 
increases rapidly as more detail is included. 
The increase in time taken to calculate the 
complete forecast is still larger because with 
more detail smaller time-steps have to be 
taken. Without going into details it is evident 
that we must be content with quite moderate 
values of # and it is important that we should 
choose an efficient parametric representation. 
To give an idea of the magnitudes involved 
it may be mentioned that the writer estimates 
that, using the most rapid method known to 
him the value n =13 is about the largest that 
(corresponding to a= 2) a high-speed electronic 
computer could cope with in order to prepare a 
forecast within a few hours. It is important to 
note that this is based on the assumption that 
the method must not lead to serious errors 
in the forecast parameters due to computational 
faults. Computational methods and errors 
associated with them will not be discussed 
here as they wili form the subject of a sub- 
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sequent paper. It will suffice to note that most 
of the methods at present in use suffer from 
very serious computational errors: the number 
of parameters whose variation, over a rea- 
sonably long time interval, is reasonably 
accurately computed (as compared with the 
true variation of the model) is much smaller 
than the number of grid points and often very 
small indeed. Thus whatever method we use 
the amount of true information is not very 
large and it is in fact actually harmful to try to 
include more detail than we can cope with. 

In barotropic motion (or, more generally, in 
motion without horizontal temperature gra- 
dients) there is no variation of motion with 
height so that only two space dimensions are 
involved. In baroclinic motion not only is 
there vertical motion but even the horizontal 
motion is a function of three dimensions. 
However, there is considerable synoptic evi- 
dence that for large-scale motion the major 
variations in the horizontal motion are describ- 
ed by supposing the thermal gradient inde- 
pendent of height. Hence we can construct 
‘a model of the horizontal part of the motion 
with two 2-dimensional fields, the field of 
the mean motion (averaged with regard to 
height) and the field of relative motion, the 
mean thermal wind (also averaged with 
regard to height). Alternatively by adding 
and subtracting a suitable multiple (either 
constant or at least independent of the fields 
of motion) of the relative motion to the mean 
motion we can express the model in terms of 
the motion at two representative levels. Now 
in horizontal barotropic motion the motion is 
(to a very close approximation) non-diver- 
gent so that in this case the mean motion 
(here the actual motion at any level) may be 
represented by a (scalar) stream function ¥. 
In the case of the model with both mean and 
relative motion we shall find that each of 
these may be regarded as approximately non- 
divergent so that now in addition to the stream 
function for the mean-motion, #, we have the 
stream-function for the relative motion, ®. 
These will be regarded as (at a given time) 
independent, though of course their changes 
are interrelated. The idea of using a model 
of this type appears to be due to SUTCLIFFE 
(1947) who has described one in which there 
are two representative levels. The method of 
derivation of the equations given below is 
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in some respects rather similar to Sutcliffe’s 
though there are important differences in the 
way approximations are made, in the ancilliary 
assumptions and in the presentation of the 
results. The most important difference is the 
inclusion of the effect of vertical motion on 
the Y-field. 

The definition of the Y and © fields serves 
merely to describe the motion succintly. In 
order to obtain a working model we have to 
derive sufficient (i.e. two) partial differential 
equations, with boundary conditions, which 


determine ae = in terms of Y and © only, 
(or in terms of quantities which can be com- 
puted when Y and ® are given). In order to 
do this it will be necessary to make further 
assumptions or postulates some of which 
may be regarded as fairly plausible, others 
somewhat arbitrary. The final test of these 
assumptions is the closeness with which the 
model simulates the behaviour of the atmos- 
phere and for this reason a “test” on a problem 
which has been solved in three-dimensions, 
with much less stringent assumptions, will be 
included. Some of the assumptions made are 
not absolutely necessary but have the advan- 
tage of simplifying the presentation. By making 
more complicated assumptions it may be 
possible to improve the fidelity of the model. 
Alternatively a similar result may be obtain- 
able if suitable (empirically determined) 
weighting factors are included. Actual use 
of the model will indicate what type of modi- 
fication is most effective. 

We shall suppose that the model represents 
motion on a large scale and that this motion 
is quasi-geostrophic, that is to say that at all 
or almost all points the geostrophic formula 
gives a fairly good approximation to the 
velocity field and also the curl of the velocity 
field. The horizontal divergence cannot of 
course be computed ‘directly but is given with 
fair accuracy by the vorticity equation, 
obtained by eliminating the pressure field 
from the equations of motion. When motion 
is on a large scale and the Richardson number 
large compared with unity we may probably 
ignore the contribution to vorticity change 
due to overturning in a vertical plane and write: 


diva v (f + curly v) + (y curly v) =o (2) 
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If at the same time we use (in place of the 
“actual” atmosphere) an incompressible model 
with the same static stability (but little varia- 
tion of mean density with height) then we have 
for the continuity equation 


Vz 
diva v + oe Sins (3) 
whence: 
Ov» d = \ 
(f+ curly v) ET (f+ curlav).. (4) 


As is well known, a much more accurate 
continuity equation, of similar form to (3) 
is obtained when instead of z we use for a 
vertical coordinate the mean pressure p and 


at the same time replace v; by v, = ah The 


dt 
only modification in (4) is the replacement of 
Ov ~ ) 
= by = and in fact the whole of the 


subsequent analysis may be carried out using a 
“pressure” vertical coordinate. Then instead 
of making our approximations in the conti- 
nuity equation we make them in the “thermal 
wind” equation. The results are similar in 
form and the only difference is in physical 
interpretation. In the present (spatial) model 
the means are to be interpreted as height- 
means with equal weighting for equal height‘ 
difference. In the alternative (pressure) model 
they are to be interpreted as pressure-means. 
It may be that something between these 
extremes gives best results in practice or alter- 
natively suitable weighting factors (relative 
to either interpretation) may be developed 
empirically. In the subsequent analysis it 
will be assumed that we are using a spatial 
vertical coordinate. 

Another approximation will now be made. 
In equation (4) we shall suppose that in most 
regions |curly v| may be neglected in com- 


parison with f so that: 
NW, ad 
fee = (f+ arlav)....... (5) 


This is certainly not true in the vicinity of a 
“jet” and the relative error may also be large 
near intense depressions or anticyclones. In 
general the errors are larger at high than at 
low levels. However, for motion on a large 
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scale we may expect to find the approximation 
not too bad in most regions. From the practical 
point of view the seriousness of the error 
depends on the extent to which it can be 
compensated for by the introduction of 
(readily computed) “weighting factors”. Here 
we may note that the error may be regarded 
as associated with incorrect “weighting” of 
dv: 
Oz 
that the appropriate weighting of other 
quantities is to some extent uncertain due to 
the incompatible requirements for simplicity 
of the thermal wind and continuity relation- 
ships. Moreover we shall subsequently need 
to adopt a rather crude and schematic distribu- 
tion for v; since our data are insufficient to 
determine its true variation. The total effect 
of all such errors cannot easily be estimated a 
priori. Hence the method will be to develop 
the simplest possible system of equations as a 
basis for experiment and subsequent critical 
re-examination. 

The stream functions Ÿ and ® are defined 
as mean values over a certain depth of atmos- 
phere, the boundaries of which will be taken 
to be +2, —2,. Thus — 2, corresponds 
to the surface of the earth and + z, the 
effective “top” of the atmosphere. In fact, of 
course, there is no top to the atmosphere but 
there is some synoptic and also theoretical 
evidence (based on perturbation theory — 
cf. Eapy 1949) that large scale disturbances 
are mainly confined to the troposphere and 
extend only a short distance into the strato- 
sphere. Hence the effective depth 22, will be 
supposed to correspond to a little more than 
the height of the tropopause but once again 
some modification may be desirable in con- 
junction with the use of weighting factors. 
ZQ is not a constant in fact but unless a simple 
approximate relation to the Y and ® fields 
can be found it may be necessary to use a 
mean value independent of the field of motion. 
This mean value certainly decreases slowly 
with latitude but the exact effect of this 
variation will not be included in the present 
analysis. The effective depth also varies 
because of unevenness of the earth’s surface 
and it is known that this is responsible for the 
existence of orographic stationary long waves. 
Since this particular variation of z, is well 
defined beforehand it is easy to include its 


in respect to height. Now we have seen 
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effects in the equations but for simplicity 
this modification will be omitted in the pres- 
ent treatment. Then if we suppose both 
the base and “top” of the effective part of 
the atmosphere as substantially flat we may 
take the boundary conditions to be vz = o at 
z=+2,. Integrating (5) with respect to z 
between these limits: 


+2 
0= [ZU + air) de ARE MG) 


By definition the velocity field in the model 
at any level z is: 


M aD, 
“us DEREN | 
EE DT] Er (7) 
en Oe nok 


approximately. This field of motion is strictly 
non-divergent but a good enough approxi- 
mation for all purposes except the direct cal- 
culation of divergence. In high latitudes Y 
and ® are very nearly constant multiples 
of the mean pressure and temperature fields 
respectively, consistent with the geostrophic 
approximation and the relatively small value 
of fi 

In low latitudes Y and ® are better deter- 
mined directly from the wind data. Some 
slight adjustments will in practice be needed 
to fit the data in middle latitudes and obtain 
the best Y and ® representation. 


Writing: 
OV ° OF 
Fe Oy TM MON 
Vy, = —- Vy ge Bau: (8) 
dy Ox 
and also: 


écurie NE Cu Curly ty = Va (9) 


we have on substituting in (6): 
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a) 0 2 
= ( r =) dy + a (& = = a) (10) 


Be di EVA DY) 
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where the contribution from if Ve = dz 


— 2 


has been neglected. Perturbation theory sug- 
gests that the contribution from this term is 
small, compared with that of the remaining 
terms, when the Richardson number is large, 
as may be verified by substituting the true 
values of the velocities corresponding to the 
development of a baroclinic wave (see EADY 
1949). We shall tentatively assume that the 
contribution from this term corresponding to 
vertical advection of vorticity may be neglect- 
ed, in more general conditions, over most 
of the region concerned. ‘The integrand in 
(10) is, so far as variation with z is concern- 
ed, a quadratic form and odd powers of z 
integrate to zero. Hence we obtain: 


d à à d 
tng + (+ YY Oy oa gree 


OE J v Pi Roar: 
3 HT "Oy 61 — 0 


and if we write: 


Die AR 
Dt eo Xo Ox I Vy, dy = at wir (12) 


corresponding to Langrangian differentiation 
following the mean motion, then 


Di. Par 
pill + So) = AC dy dy = (13) 


er) 


The term on the right hand side may be 
called the development term. It represents the 
change in absolute vorticity of a vertical 
column moving with the mean motion (but 
not of course consisting always of the same air 
particles) and expresses the difference between 
baroclinic motion and the motion of an 
“equivalent” barotropic model. We shall use 
the notation: 


0A 9B 
dy ox 


_2(48 
d(x, }) 
for the Jacobian of any two quantities A, B, 


with respect to x, y. Then for computational 
purposes (13) is conveniently written: 


d À 9B 
( — (A, B} (14) 


Ox dy 
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y, 2) _ df 


2 
a) dy ox } AZ vt + 


I 

E<— 10 7? ol EL SR 
az 

This gives one of the required equations, 


MALTE 
relating Tr © Y and ®. We now need an 


. 0® 
equation to determine et Now from the 
geostrophic relations: 

Wx gy. 
dz bag ‘i dy ; 
where @ is the logarithm of potential tem- 
perature, we obtain on comparison with (7): 


od gz Op, 00 
CORRE yak. (17) 


Avy  § & 
32 f Pr 


. (16) 


_ 820 Ip 
lan 


à ie > . 
If we ignore the slow variation of with 


f 


latitude and remember that p is a function 
of (x, y, z) whereas ®=® (x, y) by definition, 
then 


where g denotes a mean with respect to z. 
Since y varies also with z we shall write: 


> is 
@—=p+Bz— = D +Bz... (19) 

870 
where B = et the mean static stability. The 

dz 

relative variations of B are not usually very 
large and they will be ignored in subsequent 
calculations. From (19) it follows, if we 


ignore the slow variation of =, that: 
PA, 
SO 


If for the present we suppose the adiabatic 
approximation to be sufficiently accurate we 
have: 
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and the change in ® is apparently a function 
of z. However, our model is incapable of rep- 
resenting variation of this kind and all we 
require is an estimate of the mean change. 
Integrating between — 2, and + z, we get: 


+ 20 


which would be the required equation if we 
could express the mean value of v, in terms 
of ® and ¥. 

To do this we once again integrate equation 
(s) but instead of adding the contributions 
from each half of the atmosphere (z = — 2 to 
z=0 and z=o to z=+2,) we subtract 
them. Then if vz (o) is the Vertical velocity 
at the middle level: 


O + 2 


ne ff (medon) À + 
| (rs = vs.) = | 


erss)je en 


zo 


if once again the vertical advection of vorticity 
is neglected. On evaluation the even powers 
of z integrate to zero and we obtain: 


an df ie 2) 
pan — fe { 

= Vz (0) Vy, dy I (1 ax I Vy, dy 20 | 
2 à 


7 
| | = 
T ( ay | Vy, dy | 5) a, u (24) 


There is, of course, no necessary relation 
between the mean value of v, and the value of 
vz (0). However, with large scale disturbances 
we may expect to find that v, is rather simply 
distributed with, in most regions, the same sign 
for all z, a simple maximum near the middle 
level (z=0) and of course zero values at z= 
+29. The simplest distribution which satisfies 
these requirements is the parabolic one in 
which vz & (2,— 2) (2,+2). Then we have: 
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As a check we may note that the theoretical 
distribution for a growing baroclinic disturb- 
ance (EADY 1949) does not differ very greatly 
from this assumed distribution. The principal 
error arises from the fact that our model is 
not able to represent the (comparatively 
small) phase change of v, with height: in this 
example we know that the distribution of v, 
with z is not independent of x, y. However, 
this feature is associated with the fact that the 
disturbance is rapidly intensifying and it 
may be less marked in average conditions. 
We may note that an error of the same kind is 
involved in the assumption of constant ther- 
mal wind: in a growing disturbance there is a 
(comparatively small) phase change with 
height in the y-field. 


From (25) and (24) we obtain an expression 
for the mean vertical velocity. Substitution in 
(22) gives: 


NT N IR RE RE 
Dt = 3 f2 "29 [oni + (Mase tong )ee + 


| a aa 
j (0. OX r Vye dy 5) al ca (26) 
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which 


as we have seen is a pure number to be inter- 
preted as the square of the horizontal-vertical 
scale factor, then (26) may be written: 


if ap 
dy 0x 


If, as in equation (1), we write b = 


Viet Ho pl (27) 


bet Dr "Di 


where the notation of (14) has been used for 
the Jacobian. For computational purposes this 
result is conveniently rewritten : 

| DNA: | SE 


3. ~2\(a®\ dfo® 
(vs) CRETE 
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We now have a pair of equations (15) and 


. (28) 


wy 
(28) for determing —— and a) The right- 
ot ot = 


hand sides of these equations are known func- 


\ . Of 
tions so that to determine > we have a 


c 


FATPPE MDN 


~ 


00 


Poisson differential equation while os is 


determined by a Helmholtz equation. The 


sign of the constant a2 in the Helmholtz 
zo 


equation 1s such as to ensure a unique readily 
computed solution. To solve the equations 
we.need to apply suitable boundary conditions. 
If Y and ® are given over the whole of the 
earth’s surface this is a simple matter. The 


y 
condition is that Es and 2 
dt dt 


must not have singularities anywhere and this 
determines the functions uniquely. If the 
data are given over a hemisphere we may 
suppose that both functions vanish on the 
equator. With data over more limited regions 
the appropriate boundary conditions are less 
obvious and this question will be discussed in a 
subsequent paper. We may note that the 
equations (15) and (28) are easily adapted to 
computations over a spherical surface. Wheth- 
er it is necessary or convenient to do this is 
another question the discussion of which is 
postponed. 

In the above account we have neglected 
two features which in the long run must play 
an important part in determining atmospheric 
motion — surface frictional drag and influx 
and efflux of heat through radiation and con- 
vection from the earth’s surface. Formally it is 
a simple matter to include both these features 
in the model. For example if (to make a 
crude estimate) we suppose the surface stress 
proportional to and in the direction opposite 
to the geostrophic wind at z= — 2, then the 
torque acting on the column of air above will 
be proportional to the curl of the wind field 
at z= — 25, Le. to Ve, (W— ®). This. torque 
measures an additional rate of change of the 
vorticity of the column and we have only 
to add this term to the right-hand side of (15) 
to represent the effect of friction. The effects 
of heating and cooling are expressed by adding 
Q, the rate of change of » due to non-adia- 
batic heating or cooling, to the left-hand side 
of (21). Proceeding in the same way as before 


“boundary” 


we obtain an additional term, ae Q where Q 
z 


À 0 
is the mean value of Q averaged with respect to 
z, on the right-hand side of (28). Of course, 
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in order to apply this result we must have 
some simple means of estimating Q when 
Y and @ are given. 

Perhaps the best test of the formulae would 
be a prolonged series of experimental calcula- 
tions for comparison with observed behaviour. 
Alternatively we may test the model for types 
of motion where the 3-dimensional solutions 
are known. We have already made some 
comparisons with the results of perturbation 
theory. A more elaborate test may be made 
by seeing to what extent some of the quan- 
titative results of perturbation theory are 
reproduced by the model. Undisturbed hori- 
zontal baroclinic flow may be represented by a 
mean flow U= — (=) and a thermal wind 

C 
2® us 
De (52) where U and T are constants. 
0 
If now we superpose a small perturbation 
represented by stream functions W,, and @,, 


we obtain, on substituting in (15) and (28) - 


and picking out the first order terms, the 
perturbation equations: 


à 0 2 df d 
o= (us I l'AS 
ER 
+ — Ts. Ve D, docasée (29) 
and: 
= d 2 2 3 df d | 
O ( 2 +5\(v; b 10 dy 7. Pi I 


d 3 
re)” 1130) 


Corresponding to the assumptions made in 
the 3-dimensional theory (EADY 1949) we 
suppose that Vb z, = H and ay may be taken 
as effectively constant. Then the simultaneous 
equations (29) and (30) evidently have solu- 
tions of the form: 


tAx+uy+ ot) 
Le OD, = fe 


tAx+uy+ 9) 


(31) 


where E, F, À, u, are constants. Equations 


(29) and (30) will be satisfied if: 
(0,0) E+-TF=0 


[U, G + X)—U,] F+T(1—X) E=0 (32) 
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where we have written: 


= a. I df, 
el, a rn 
X= 3 = : 33) 


(242) He (B+ pe) be 


Eliminating E/F from equations (32): 
ee 2 Ur | 
(+x) (7) (2 + x) (6) De 


whence: 


427? : 
(2+ X) + +5 xy 


w r 


(35) 


(de 2 (1 + X) 


The boundary conditions in the horizontal 
directions require that both À and w should be 


real. Hence X and (= are necessarily real 


and positive. The diccueboices will be un- 
stable if, and only if, # is complex, i.e. if 


Ww 


= is complex. The condition for instability 


; 
is therefore: 


T? 
Pere nth AD ot) 
ù ef 


Now VX is a number proportional to the ratio 


N 5 PURE, 
of the “effective” wavelength ———— to the 


VAE + u 
“dynamic depth” 2H. There will be instability 
for some wavelengths if (rearranging (36)) 


Ore HE ee (37) 


Me Ta N3 74% 

for any value of (42 + u), i.e. for any positive 
X. The maximum value of the right-hand 
side is 1/3 (when X?= 2). Hence the condi- 
tion for instability of the initial flow is: 


| He df 
ae 


ya 
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The condition may also be written in the 
form: 


(sie ere 


where s is the slope of the (unperturbed) isen- 
tropic surfaces. This result is of considerable 
interest, though it has yet to be checked by 
three-dimensional calculations. The corre- 
sponding wavelength for the disturbances 
which first become unstable as T is increased is 


given by X=V2. Eliminating H? we obtain 
I df 2 T 


(+) dy V3 


ie: (ol ele; ie 


and this result may be compared with the 
stability criterion discovered by CHARNEY 
(1947). 

The simplicity of the above calculations 
illustrates the power of the 2 %-dimensional 
model. As a more precise test of its accuracy 
we shall consider the case when i is neglected. 
Then U,= 0 and in place of (34) we have: 


AR, I— xX 
eT pie 


In this case there always exist unstable waves, 
as is evident also from (38). The condition to 
be satisfied is that X > 1 or: 


I I ar I 7 
——— — + b “SS . 2 
era 


This compares with the true value: 


bee 


I I 
2 — ss x 
VA? + u? 1.1997 


(Eapy 1949 p. 39). The only error is in the 
numerical factor and this is clearly of a kind 
which could be eliminated by modifying the 
constants in the equations (15) and (28). 

It is easily verified that the disturbance of 
maximum growth rate (maximum imaginary 
part of #) corresponds to u=o and maximum 
Us 
x" 
gives the corresponding wavelength so that: 


.. (43) 


From (41) we find that X= 1+ V2 


Ee aE AWD py) 
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as compared with the true value: 


(44) 


I I 


À 08031 Vbzo. 


- (45) 


The corresponding value of # is purely imag- 
inary as in the accurate calculations. The 
numerical magnitude deduced from (44) and 


(41) is: 


= IF I: 
= *—— ~ 0.4142 —— (46 
ala one, 
as compared with the true value: 
if 
Ol = 0.3008 Te Mr den eee 47 
ja] =0.3008 (47) 


The value of r computed for the model is: 


0 | I 
AN Va Vi + va 


which compares favourably with the true 
value: 


20.3720 8 (48) 


3) == 0.285 


= ve a ea ave ele ee tre ete 


The above is probably a fairly severe test of 
the model. For comparison we may note that 
the assumption that temperature is horizon- 
tally advected (neglect of vertical motion) 
leads to values of |9| which are very much too 
large for the short wave-lengths — in fact 
there is no minimum wave-length for instability. 
For very long waves the effect of vertical 
motion is much less but on the other hand it is 
precisely in these conditions that an equivalent 
barotropic model can be constructed to repro- 
duce the behaviour of the mean flow with 
some degree of accuracy. 

The writer would like to express his thanks 
to Professor C.-G. Rossby and members of 
Institutet för Meteorologi, Stockholm, for 
lively and constructive criticism, during a 
recent visit, of the ideas presented here. 
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Abstract 


The logical design of an electronic digital computer is described. A simple computation 
problem: is proposed; this is used to illustrate the main steps leading to formulation of 
a problem for machine computation — the propram and the flow diagram. In conclusion, 
some general comments are made concerning areas of application of high-speed computers 


in meteorology. 


High-speed digital computers now are to 
some extent available to meteorologists for 
research. Doubtless there will be intensive 
developments in certain fields of theoretical 
and applied meteorology where heretofore 
progress has been inhibited by the prodigous 
amounts of computation required to yield a 
useful result. As this work advances, it is 
likely to take on some of the aspects of an 
engineering science, and the underlying prin- 
ciples of the use of computers for meteorolog- 
ical problems may tend to become obscured 
by a superstructure of special methods of 
computation and procedure. 

It is inevitable that high-speed computers 
will open new avenues of investigation in 
theoretical studies of the atmosphere, and 
although among meteorologists there is wide- 
spread latent interest in these prospects, it is 
probably too early to predict when to expect 
results of practical value — say, to the fore- 
caster. Nevertheless, it seems rather evident 
that to maintain a well-balanced approach to 


problems that offer some promise of yielding 
to computation methods, there should be 
effective communication between those who 
gain experience in the techniques of computa- 
tion, and those whose contacts are mainly 
with applied meteorological problems — for 
example, in the fields of forecasting or syn- 
optic meteorology. 

With this thought in mind, the writer has 
undertaken to sketch a broad view of some 
of the steps that enter into the task of domesti- 
cating a high-speed computer for meteorolog- 
ical purposes, and although the view is 
rom a very limited level of experience in 
these matters, it may nevertheless help to 
illuminate for the meteorologist this novel 
area of research and to implement his interests 
in a field that he may sense is intriguing but 
which otherwise may seem obscure. The 
results of this endeavor are reported in what 
follows; they are grouped under four headings: 
(1) acomputer, (2) a problem, (3) a computa- 
tion, (4) concluding remarks. 


HIGH-SPEED AUTOMATIC COMPUTERS AND THEIR USE IN METEOROLOGY 


The substance of this discussion has been 
taken from the reports by Goldstine and von 
Neumann, to which reference is made at 
the end of the paper. The reader should 
consult these reports for a more systematic 
and thorough treatment of the logical design 
of a high-speed computer and a more detailed 
exposition of methods of preparing problems 
for computation. 


1. A computer 


An electronic digital computer is an ex- 
tremely complex instrument containing many 
thousands of components and involving 
application of the most advanced techniques 
of electronic engineering. Fortunately, it is 
possible to comprehend the logical structure 
of the machine without understanding the 
design and construction whereby this logical 
structure is achieved, just as it is possible to 
know the logical structure of an automobile 
without understanding the construction of 
the various components or even the manner 
in which they accomplish the end result. 
Indeed, this is not surprising because in point 
of fact, the logical structure of the computer 
must be formulated more or less independently 
of the specific engineering techniques by 
means of which its physical counterparts 
ultimately are realized. 


The basic logical elements of an electronic 


digital computer are (i) the arithmetic organ, 
(ii) the memory, (iii) the control. In addition, 
there must be provision for transfering data 
into and out of the machine — namely, an 
input-output device. The significance of the 
primary logical elements can be illustrated by 
a simple example. Let a and b stand for two 
numbers which, at a certain stage in a partic- 
ular computation, must be combined arithme- 
tically to yield a derived number, F (a, b); 
for example, F (a, b) = a+b, or a—b, or 
axb, etc. Initially, a and b are stored in the 
memory; the control arranges their transfer 
to the arithmetic organ, instructs this organ to 
produce F (a, b), and returns this derived 
number to the memory (see figure ı). (Cer- 
tain operations involve only a single number; 
for example, F (a)=—a, or F (a) = |al, or 
F (a)= the first four digits of a, and so on. 

The arithmetic organ contains a small 
number of registers — analogous to the “dial” 
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Figure 1. The basic elements and their functions. 
registers of an ordinary desk calculator — 
each of which can hold one number (that is, 
a series of digits) at any instant. Not many 
registers are required if, as in most computers, 
the arithmetic organ performs only one 
operation at a time. When a number held in a 
register is no longer needed, the control can 
erase it, or send it to a specific location in the 
memory. If the specified memory location 
already holds a number, this number is erased 
and the number sent from the register takes 
its place. 

The usefulness of a high-speed computer 
depends to a large extent on the versatility 
of the control, and the capacity of the memory. 
Every computation of sufficient scope to 
justify the use of a high-speed computer will 
involve a great many intermediate numbers 
that do not represent the final result but are 
produced at some stage of the computation 
and set aside to be used at a later stage. Further- 
more, in some problems the initial data as 
well as the final computed result may in 
themselves involve a great many numbers. If 
the machine is to function as a high-speed 
computer, it is clearly essential to limit sharply 
the necessity for stopping the machine and 
supplying or removing numbers before the 
computation is finished; indeed, ideally, the 
machine should not be interrupted at all 
during the. course of the computation. To 
provide for a completely automatic sequence 
of computations, there must be a fully auto- 
matic control that can process the necessary 
instructions (orders) without interruption, and 
a large-capacity memory that can accommo- 
date simultaneously a great many numbers. 

One of the simplest types of memory to 
visualize is built around a specially designed 
electron tube, the so-called “electrostatic 
storage tube”, which in construction and 
operation is similar to a small television tube 
(iconoscope). At one end there is a network 
or mosaic of electrostatically charged “spots”. 
Each of these small charged areas can be 
represented schematically in the manner shown 
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Figure 2. Schematic illustration of the networks of 
charged points in five identical electrostatic-storage mem- 
ory tubes. The circles enclose corresponding points. 


in figure 2. Tube 1 is pictured as a network of 
16 points, each point representing a small 
area of electrostatic charge; tube 2 is identical, 
and so on. A number (series of digits) can be 
represented as equivalent to a group of 
charges located at corresponding points in the 
series of networks; for example, the circled 
point in tube 1 holds the first digit of the 
number, the corresponding point (circled) in 
tube 2 holds the second digit, and so on. The 
arrangement illustrated in figure 2 would in 
this way accommodate 16 five-digit numbers 
(16 points in each tube, and 5 tubes). In prac- 
tice as many as 40 tubes are used, each with as 
many as 1,024 points (a 32x32 network 
instead of the 4 x 4 pictured above); a memory 
of this size would have a capacity of 32 x 32= 
1,024 forty-digit numbers. 

In each tube is a device for producing a 
narrow beam of electrons that impinges on 
the network of changed points, and in fact 
can “scan” this network point by point and 
row by row. The required deflections of the 
electron beam are produced by a properly 
controlled magnetic field, which means that 
the scanning process can be performed with 
extreme rapidity. The tubes are synchronized 
so that at any instant their respective beams 
impinge on corresponding points of the 
charged networks; in this way the various 
digits in a given number can be “sensed” 
simultaneously and transferred to a register. 
By proper control of the intensity of the beam, 
it is possible to arrange for the scanning pro- 
cess to read, or to write, which means in 
effect that a number can be transferred from 
the memory to the register, or vice versa. 

The computations required to solve a given 
problem can be formulated as a list of orders 
that enumerate step by step, every single 
arithmetic manipulation that must be per- 
formed, including all necessary transfers of 
numbers to and from the memory. The 
orders are then translated into numbers 
according to a definite code, and as such are 
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themselves stored in the memory (thereby 
establishing the link between human control 
and machine control of the computation). 
In this sense the memory holds logical infor- 
mation as well as strictly numerical informa- 
tion; both types of information are expressed 
as series of digits). Since the orders are now 
located in the memory, they can be scanned 
and transferred, so to speak, to the control 
organ, which arranges for these orders to be 
performed, one after the other, as they are 
read. 

Virtually every problem of a scope that 
justifies the use of a high-speed computer 
will demand a highly versatile control organ. 
For example, a commonly recurring situation 
at various stages of the computation requires 
repetition (iteration) of the same sequence of 
operations and hence of the same sequence 
of orders. To illustrate: suppose we have, 
say, fifty pairs of numbers and that the same 
set of operations must be performed on 
each of these number pairs. Obviously, in 
writing out the orders we cannot afford, as 
a rule, to repeat the same sequence over and 
over again the desired number of times be- 
cause the memory space available for these 
orders probably would be quickly exhausted. 
It should be sufficient to list this sequence of 
orders just once; so the control should be 
able to “go back” to the beginning of the 
sequence and repeat as many times as required. 
In fact, we should generalize this requirement 
in the following way: normally, the control 
will take up the orders in sequence as they are 
entered in the memory; however, at any 
stage by means of an order that “transfers” 
the control, it will proceed not to the next 
order but to any other designated order, 
and there continue in sequence. 

There is another, vital feature of versatility 
of the control. Where a transfer of the control 
is required, we should in principle specify 
beforehand precisely the order to which the 
control should be transferred; or, if a sequence 
of orders must be iterated, we should specify 
beforehand exactly how many iterations must 
be performed. There are several difficulties 
here, connected with the fact that very fre- 
quently this type of information is not known 
until the computation has proceeded up toa 
certain point. If we are to avoid interrupting 
the computation, such decisions (to transfer 
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the control) must be made internally on the 
basis of specific orders provided for this 
purpose. Indeed, even when for example we 
know in advance exactly how many iterations 
of a given sequence of orders are needed, 
there would be no convenient way to ter- 
minate these iterations unless a decision can 
be made internally. All “decisions” of the 
type required can be reduced to an order 
providing for conditional transfer of control, 
depending on the sign (plus or minus) of a 
number held in one of the registers. Suppose, 
for example, that the first iteration of a sequence 
of orders produces a number a,, the next 
iteration a slightly smaller number ag, the 
third a still smaller number a3, and so on, 
and we want to iterate until a number, say a,, 
is produced that is smaller than a specific 
preassigned number b. This is accomplished as 
follows. After each iteration there is an order 
placing the difference a, —b in one of the reg- 


isters, followed by the conditional transfer : 


order. If this difference is positive (meaning 
a, is still greater than b) the control will 
transfer back to the iteration; if it is negative 
the control will proceed out of the iteration 
to whatever order may follow in sequence. 
This provision for conditional transfer of 
control is used repeatedly in most problems. 
It provides for extremely versatile and fully 
automatic control of the computation. 
Another typical application of the conditio- 
nal transfer order is this: suppose we know in 
advance exactly the number of iterations 
required for a given sequence of orders — 
say, 30 iterations. Before the first run through 
the sequence the number 1 is ordered into 
an empty location in the memory, and just 
prior to each repetition of the sequence this 
“counting index” or “iteration index” is 
advanced by adding 1 to the number already 
there, so that at the end of each sequence this 
particular memory location will hold a num- 
ber exactly equal to the number of times 
the sequence has been completed up to that 
stage. Further, at the end of each sequence 
the difference “29 minus the index” is ordered 
into a register and then the conditional transfer 
order follows. When this difference is positive 
(or zero), meaning that 30 iterations (the 
desired number) have not yet been completed, 
the control is transferred back to the beginning 
of the sequence; but after the 30th iteration 
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this difference is negative, so the conditional 
transfer order carries the control out of the 
iteration and on to the next part of the com- 
putation. 

We have now discussed the main features 
of the logical framework of a computer; they 
are (apart from an input-output device): 

(i) An arithmetic organ, consisting of a 
few registers each to hold a single number 
(series of digits) at any instant. The arith- 
metic organ is responsible for performing 
basic arithmetic operations such as addition, 
multiplication, division. 

(ii) A memory, capable of storing a great 
many numbers, with provisions for reading 
and altering its contents as the computation 
proceeds. 

(ii) A versatile control, which responds 
to a list of orders stored in the memory. 

A final comment regarding the speed of 
computation: in an electronic computer all 
internal operations are purely electronic and 
therefore extremely rapid. Indeed, the unit 
of time most convenient in this connection 
is the microsecond — one millionth of a second. 
A good approximation to the duration of 
most types of computation can be obtained 
by allowing about 500 microseconds for each 
multiplication (or division), and about 65 
microseconds for each order (including multi- 
plications and divisions). For example, in a 
typical meteorological calculation (the 24- 
hour barotropic forecast) carried out on the 
Princeton computer, 1,960,000 multiplications 
and a total of 42,300,000 orders are required. 
The approximate duration of this computation 
is therefore 500 x 1,960,000 +65 X 42,300,000 
microseconds, or about 62 minutes. This makes 
it quite clear why an extremely involved and 
lengthy computation — hopelessly imprac- 
tical to attempt by hand — may be quite 
practical for the high-speed computer. 


2. A problem 


The problem that will now be described 
has been chosen mainly to illustrate some 
typical steps leading to formulation of a com- 
putation for the high-speed computer. It is 
in fact a simple problem and requires only a 
rather small amount of computation, so it is 
not representative of the type of problems 
for which a high-speed computer is best 
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Figure 3. Polar stereographic projection; the dashed lines 
are latitude circles, the solid lines delineate the region over 
which the vorticity will be computed. 
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Figure 4. A network of grid points covers the rectangular 
region of figure 3. The grid shown here has 18 interior 
columns and 13 interior rows. 


suited — namely those that involve extra- 
ordinary amounts of calculation. 

Suppose we have an analyzed soo-mb map 
showing height contours. Problem: to calculate 
the geostrophic vorticity field that corresponds 
to the contour analysis. Let z denote contour 
height of the soo-mb surface, g the acceleration 
of gravity, and f the Coriolis parameter, so 
that the geostrophic wind components are 


Then the geostrophic vorticity is 


an 00, ON 
a en NR 


rz 
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d y? 
neglecting here the variation of f in the de- 
nominator of the expression for u. 


The contour analysis is drawn on a map 
that represents a plane projection of the 
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spherical earth; this projection will distort true 
distances. Let Ax be the true distance (on the 
earth) between two neighboring points and 
let Aa be the distance between these same 
points measured on the projection; then 
Aa= m Ax, where m is the scale factor for the 
projection and will usually depend upon the 
latitude where the measurement is made. 
To be specific, assume that a polar stereo- 
graphic projection is üsed; the scale factor for 
this projection is 


2 
re +sin®’ 
where ® is the latitude. Writing Ja=m Ax 
and Ab=mAy for distances measured on the 
projection, the equation for calculating the 
vorticity from a contour analysis drawn on 
the projection is, from (T) 


m? 92 z 9° 2 
g ( 5 =). (2) 


+ da? 
It will be recalled that the Coriolis parameter 
f is equal to 2@ sin ®, where « is the angular 
speed of the earth’s rotation. 


Figure 3 is a portion of a polar stereographic 
projection on which a rectangular boundary 
has been marked out, to delineate the region 
over which the geostrophic vorticity will be 
computed. (The contour analysis is not shown.) 
We now construct a network of grid points 
(figure 4) covering the rectangular region, 
such that the distance between grid points 
is the same in either diretion — say equal to 
d. In order to specify the location of a grid 
point in this array, let i designate the column 
in which the point is located (counting from 
left to right), and j the row (counting from 
bottom to top); thus, i=o designates the col- 
umn coinciding with the left-hand boundary, 
i=1 the first interior column, i = 2 the next, 
and so on; similarly, j=o designates the row 
coinciding with the bottom boundary, j = 1 
the row next above this, and so on. 


From the contour analysis, values of the 
height z must be interpolated (this can be 
done visually) for each of the grid points. 
The interpolated height at the grid point 
whose coordinates are i, j we denote by 2j. 
Consider a typical group of five neighboring 


HIGH-SPEED AUTOMATIC COMPUTERS AND THEIR USE IN METEOROLOGY 


grid points (figure 5), the central point with 
coordinates i, j; we wish to calculate the geo- 
strophic vorticity ¢; at this central point. 
Returning to formula (2) we see that the 
vorticity depends upon g (a constant, the 
same for all points), upon the scale factor mj 
at this point, the Coriolis parameter fi; at 
this point, and the second derivatives of the 
contour height in both directions at this 
point. 

To evaluate the latter two terms, consider 
the point half way between the central grid 
point (i, j) and the next grid point to the 
right (i+ 1,7); there, the first derivative dz/0a 
is approximately 

CRE, 
where d, as previously stated, is the distance 
between neighboring grid points. Likewise, 
half way between i, j and i— 1, j the first 
derivative dz/da is approximately 


(z..—z,__ .)/d 
1,7 EE, ‘ 


The second derivative 92?/9a? at the central 
point (i, j) is then approximately the differ- 
ence between these two first derivatives, 
divided by d (the distance between the two 
“half-way” points); this reduces to 


Cr 


el re 
ER 9 1,J 
The second derivative 022/0b? is approximated 
in the same way, by considering the points 
i, +1 and i, j — 1 above and below the cen- 
tral point; this gives 


(2, jar Tiji 2%, p/P. 


i, j= oF 
Adding the two second derivatives, as required 
in formula (2), we get 


= 2 
(a5 <i ie 47, ‚ld. 
The terms in parentheses consist simply of 
the sum of the height values at the four 
points surrounding the central point, minus 
four times the height at the central point; 
for brevity we will denote this group of 
terms by Sj. . 
The vorticity at i, j is now to be calculated 
from the formula 


ie 


era 9) 
ij 


173 


j+l 


itl 


Figure 5. A typical grid point i, 
J and its four neighbors. 


in which g and d are definite constants, mj; and 
fi; depend upon the latitude of grid point i, j 
(strictly, upon the sine of the latitude), and 
S;; depends upon the contour heights at i, j 
and the surrounding four grid points. 


There remains only to find the formulas by 
which mj; and fi; can be calculated for each 
i, j. Let r denote the linear distance between 
the pole and a point at latitude ® on the polar 
stereographic projection; the properties of 
this projection require that 


where a is the radius of the earth. If ip, jo 
denote the grid coordinates of the pole, then 
the square of the linear distance (on the pro- 
jection) between pole and grid point i, j is 


r= d? [Go +?) 


according to Pythagoras’ theorem. The last 
two formulas provide for calculation of sin ® 
at any grid point i, j: we need only the values 
of i and j and the constants jg, jo, d, and a. 

From the formulas already given, the reader 
can easily verify that the following computa- 
tion scheme will yield the geostrophic vor- 
ee 


Formulas: 
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Figure 6. Flow diagram for vorticity computation. The formulas for S;; and Ci in box III and Cy in 


box IV are listed in the program given previously. 


Data: 
zÿ= contour heights at all grid points. 
DORE D Pau Ont end | 1) 


Constants: 
g = acceleration of gravity, 
@ = angular speed of earth’s rotation, 
a = radius of earth, 
d = distance between grid points on projec- 
tion, 
ip, Jo = grid coordinates of pole, 
p == 


q = number of interior grid rows. 


number of interior grid columns, 


Roughly 10 multiplications are required to 
produce &;; for one grid point from these 
formulae. Suppose, further, that something 
like a total of 100 orders (including nonarith- 
metic operations, such as necessary transfers 
of numbers between registers and memory) 
are required to produce one £;. Allowing 
soo microseconds for each multiplication 
and 65 microseconds for each order, we get an 
estimate of 500x 10+65 x I00=11,500 mi- 
croseconds, or say roughly 0.012 seconds as 
the time required to compute the vorticity 
at one grid point. For a grid of 500 points 
the whole computation would be finished in 
about 500 x 0.012=6 seconds — which justi- 
fies our earlier remark, that this particular 
problem is not representative of the type for 
which the high-speed computer should be 


used, where the computation time may be 
of the order of minutes or even hours. 


3. A computation 


The computation scheme outlined at the 
end of the preceding section constitutes what 
is known as the program. In order to portray 
the precise sequence of logical steps by which 
the control will carry out this program, we 
construct a flow diagram. The flow diagram 
for this computation is presented in figure 6. 

Before interpreting figure 6, it is necessary 
to agree on a plan for enumerating the mem- 
ory locations, so that the position of each 
number held in the memory can be specified. 
It will be recalled that the various digits 
belonging to one number are distributed 
among the memory tubes; each digit is 
assigned to one tube and the digits of one 
number appear in corresponding points of 
the charge networks of the various tubes 
(figure 1). Suppose each tube accommodates 
1,024 points ees) — so the capacity of 
the memory is 1,024 numbers — and that these 
points are arrayed in a square network of 
32 rows and 32 columns. Each point in the 
array is a definite memory location. We 
enumerate these locations as follows: start 
in the lower left corner and call this point 0; 
proceed serially along the bottom row to 
31 in the lower right corner; continue with 
the first (left-hand) point in the second row 
as 32 and procced along this row to the end 
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point 63; likewise, 64 will be the first point 
in the third row, and so on; the very last 
point in the array (in the upper-right corner) 
will be number 1023, so that (since we started 
with o) there are 1,024 locations in all. 

Turning to figure 6, we note first that the 
flow diagram consists of a collection of “boxes” 
connected by branches marked with arrows. 
Each box (except those marked #) stands 
for a group of orders necessary to perform 
the operations indicated by notations within 
the box. The arrows on the branches connecting 
the Loxes indicate the sequence in which the 
control will take up these groups of orders. 
Starting at the place marked “begin”, the 
control proceeds to the first box, numbered I. 
The notation “1 to 1” in this box is intended 
to mean that the number 1 is ordered into 
memory location 1. The next box (marked +) 
does not stand for any operations or orders, 
but is intended to make an assertion as to how 
the number 1 introduced in the preceding 
box is henceforth to be interpreted; the nota- 
tion “I > j” signifies that we are to regard the 
number 1 (in memory location 1) as a value — 
in fact the first value — of the index j. This 
type of purely logical assertion is useful as 
an aid in reading and interpreting the flow 
diagram. 

Before proceeding to box II, we remind 
the reader that i and j identify the column 
and row location of a point in the grid (not 
to be confused with the charge network in 
the memory!) that was marked out on the map 
projection (figure 4). Interpolated contour 
heights z;; have been assigned to each of 
these grid points, and our objective is to 
compute the vorticity ¢, at each point. 
To fix ideas we will refer to a grid with 18 
interior columns and -13 interior rows (figure 
4), so that i runs from o (the left-hand bound- 
ary) to 19 (the right-hand boundary) and j 
runs from o (the bottom boundary) to 14 
(the top boundary). This grid has 18 x 13 
= 234 interior points; therefore, the results 
of the computation will consist of 234 values 
of ¢. The data required for the computation 
consist of the values of z at every grid point 
including points on the boundaries; therefore, 
20 x 15 =300 values of z comprise the data. 
(The four corner points are in fact not needed 
but it is convenient to include them in order 
that every row may have the same number 


Table 1. Storage table, 


Memory | Quantity Remarks 
location stored (see fig. 6) 
o = 
I 7 Formed in boxes Iand X 
2 i » » » II » VII 
3 Sij Computed in box III 
4 ij » ee dye ual 
Sage? a, 
20 o 
21 it 
22 (d/2a)? | Used in box III 
23 g/2 © d? DA) ON) SIN 
24 io ee LIT 
25 Jo DR GR EI 
26 17 Yat ae A 
27 12 yee ay VAT 
28— 31 _ 
32—265* Ë Computed in box IV 
266—287 == 
288—587t z Data 
588—607 == 
608— ? Orders 
2023| = 


* bi; is stored in location 13 + i + 18). 
ii aij » » » » 288 + i + 20/. 


of points.) In table 1 we have compiled a 
“storage table’, which is a record of the 
memory locations assigned to the numerous 
quantities involved in the computation. The 
234 values of € will be stored in locations 32 
to 265, while the 300 values of z are stored 
in locations 288—587. 

Proceeding to box II, we have a group of 
orders that send the number 1 to memory 
location 2; the following box is an assertion 
that this number is to be interpreted as a value 
of index i. In box III we compute two numbers 
Sÿ and C; needed in the formula for 6; 
(note that i and j are the nunibers in memory 
locations 2 and 1, respectively); Sj; is sent to 
position 3 and C;; to position 4 in the memory. 
We can now calculate ¢;; (box IV), which is 
sent to position 13+i+18j; for example, 
Cra (i=1, f=1) is sent to position 13 +1 +18 
= 32. The reader can easily verify that this 
scheme of storing the successively-computed 
values of ¢ in the memory will place the re- 
quired 234 numbers serially in memory posi- 
tions 32 to 265 inclusive. 


1 The last three paragraphs of this section give a more 
detailed account of the computations in boxes III and IV. 
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At this stage the computation of one of the 
é’s is completed — namely, the one corre- 
sponding to i=1, j= 1. The flow diagram is 
arranged in such a way that the ¢ ’s are taken 
up in the following order: starting with i=1, 
j=1 we proceed along the first interior row to 
the last interior point (i=18, j=1), then to 
the first interior point of the next row (i=1, 
j=2) and along this row to i=18, j=2, and 
so on; the last & computed will be for i=18, 
j=13. Now the orders involved in boxes III 
and IV can be used for every £, so the control 
should return to box III each time we pass 
on to the next grid point. In order to control 
this iteration we must alter the values of i 
and/or j (in memory locations 2 and 1) so 
they will always correspond to the grid point 
whose & is being computed. Further, every 
time we complete the sequence of orders in 
boxes III and IV we must determine whether 
at that stage we are at the end of one of the 
interior rows (i=18), because if we are, one 
of two steps must be taken: if we are at the 
end of the last row (j=13), the iteration must 
be stopped because the computation is finish- 
ed; if we are at the end of a row but not 
the last row, we must advance j and calculate 
for the next row, starting with i=r. 

The orders that provide for this decision 
are indicated in boxes V and VI. In box V, 
the number 17 — i is sent to one of the regis- 
ters, designated R,. Box VI represents the 
conditional transfer order described previously. 
If 17 — i is positive (or zero) then i < 18 and 
we are not at the end of a row; in this case the 
control is transferred to box VII (along the 
“positive” branch of box VI), where the index i 
is advanced so that we proceed to the next 
point in the same row. From box VII the con- 
trol is transferred unconditionally back to box 
II. When 17 — i first becomes negative we 
know that i = 18, the end of a row; we must 
then determine which of the two steps men- 
tioned above must be taken, that is, we must 
determine whether we are on the last row. 
This is accomplished in box VIII, which orders 
12 — j to register R,, and box IX, which is 
the conditional transfer order. If 12 — j is 
positive (or zero) then j < 13 and we are 
not on the last row; in this case the control 
is transferred to box X (along the “positive” 
branch of box IX), where the index j is ad- 
vanced, and then unconditionally to box II, 
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where i is returned to its initial value 1, the 
first step in starting along the next row. 
When 12 — j first becomes negative we know 
that j=13, the last row, and that all 2’s have 
been computed. (Recall that box IX is reached 
through the “negative” branch of box VI, 
so it is entered only when i=18, the end of a 
row.) The “negative” branch (j=13) of box 
IX proceeds to box XI, a halt order. 

It is evident that the flow diagram clearly 
portrays the sequence of operations by which 
the computation is performed, including the 
various points at which and to which the 
control must be transferred. The control 
cycles repeatedly around the loop (induction 
loop) formed by boxes III, IV, V, VI, VI, 
except that every 18th cycle carries it through 
the “larger” (less frequently traversed) loop 
formed by boxes III, IV, V, VI, VIII, IX, X, IL. 

The next step after completing the flow 
diagram is to write out the orders needed to 
perform the operations implicit in each box 
of the flow diagram. The complete list of 
orders required is the code. It will be recalled 
that the orders, properly translated into digits, 
will be stored in specific memory locations; 
referring to table 1, the reader will note that 
the orders for this computation are assigned to 
a portion of the memory beginning with 
location 608. The writer has not prepared the 
code corresponding to the flow diagram of 
figure 6, so the exact number of orders needed 
has not been determined. 

There is one further aspect of the computa- 
tion in boxes III and IV that should be men- 
tioned. In box III the quantity S;; is calculated, 
and this depends upon the values of z at i, j 
and at four points surrounding i, j (figure 5). 
The five values of z needed to compute Sj 
must be called out of their memory locations 
by means of appropriate orders that will 
specify the “addresses” of these z’s. Each 
time we pass through box III we are dealing 
with a different point (that is, a different pair 
of values i, j), which means that — as is 
obvious — a different set of five 2’s is required 
each time. Therefore, box III must include a 
provision for changing addresses in the orders 
that call out the z’s; this is certainly feasible 
because each time we enter box III we know 
what the current values of i, j are (they are 
stored in locations 2 and 1), and we also 
know — according to the previously stated 
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scheme for storing the z’s — that the address 
of 2); is 288+i+20/ (see table 1). This is an 
interesting and important feature of the flexibil- 
ity of the control: by means of appropriate 
orders within the code, other orders can be 
changed whenever such changes are needed to 
suit the requirements of the computation. 

In box IV a similar change of address must 
be included, in order to arrange for the storage 
of each newly-computed £;, in its previously 
designated location. According to table 1, 
Éi should be stored in location 13 +i+18 j; 
therefore, before giving the order to store {;,, 
the address in this order must be changed 
to correspond to the current values of i 
and j. 

Finally, we may note that after S;; and C; 
(computed in box III) are used to calculate 
&; in box IV, they will no longer be needed, 
SO it is not necessary to store these two quan- 
tities permanently. In fact, each time Sj; and 
C;; are computed and stored in locations 
3 and 4, the values previously computed (and 
stored there) for the preceding point are 
automatically “erased”’. 


4. Concluding remarks 


The central problem of applied meteorology 
is weather forecasting in all its phases. One is 
naturally led to speculate on the extent to 
which high-speed automatic computers can 
be brought to bear on this problem. In the 
first place, it is necessary to repeat a comment 
made earlier: a computer of the type we are 
considering excels in dealing with computa- 
tions that involve large quantities of data and 
(or), especially, prodigous amounts of arith- 
metic. The question is then, to what extent 
are advances in various phases of the problems 
of weather forecasting at present impeded by 
our inability to cope with difficulties of the 
type just mentioned? We will touch briefly 
only a few aspects of this question. 

As a specific example, consider the problem 
of forecasting radiational temperature changes 
in the free atmosphere. A great deal is known 
about the theory of such changes: they are 
controlled by the distribution of water vapor, 
carbon dioxide, ozone, and other principal 
absorbers of long- and short-wave radiation, 
and by the distribution of liquid water (clouds), 
as well as the reflectivity and temperature 
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of the underlying ground surface. Numerous 
graphical, tabular, and other approximate 
methods have been devised to account for 
most of these factors, but a more thorough 
treatment of the problem probably rests on 
detailed integration of the complete radia- 
tional transfer equations. The complexity of 
these equations is staggering, when one 
attempts to provide in detail for all of 
the factors that our present knowledge indi- 
cates are important, such as the distribution 
and intensity of lines in the spectra of the 
principal absorbers, and their variations with 
pressure and temperature, to mention only 
one element of the problem. A direct attack 
on the problem, attempting to account com- 
pletely and simultaneously for all known 
factors, probably would be fruitless. However, 
it seems likely that a great deal could be 
learned by using a high-speed computer to 
intégrate the transfer equations under condi- 
tions somewhat more complex than can reason- 
ably be treated by means of “hand” computa- 
tions. 

As another example, take the forecast of 
high-level winds (over relatively short periods, 
say 24 hours or less) — for example, for use 
in operation of modern aircraft. Here we are 
face to face with the problem of predicting 
the future state of motion of the atmosphere 
— a problem inherent in virtually all phases of 
weather forecasting. From the time of Helm- 
holtz, and particularly following the develop- 
ment of the Norwegian school of physical 
hydrodynamics under V. Bjerknes, meteorolo- 
gists have on the whole had implicit faith in 
the belief that the laws of hydrodynamics, dili- 
gently applied to the atmosphere, will advance 
our understanding of the motion of the atmos- 
phere and thereby improve our ability to pre- 
dict its future state. Certainly it is true that our 
knowledge of the behavior of the atmosphere 
is vastly improved as a result of the approach 
through theoretical hydrodynamics, but be- 
cause of the overwhelming complexity of 
the atmosphere, the advances as brought to 
bear on forecasting problems have been 
largely qualitative. It must be said that the 
meteorologist’s faith in the efficacy of hydro- 
dynamical methods has not been genuinely 
rewarded: modern forecasting techniques in 
most of the world’s weather services probably 
are preponderantly subjective, and for the 
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most part do not reflect appreciably in quanti- 
tative terms, the advances in theoretical “dyna- 
mic” meteorology. 

In these circumstances we are in fact en- 
countering the difficulties mentioned at the 
outset of this discussion: strict application of 
hydrodynamical methods entails an unearthly 
amount of computation, as well as vast 
quantities of data, and is far beyond the reach 
of “hand” computation. Perhaps it is not un- 
reasonable to advocate renewed faith in hydro- 
dynamical methods of forecasting the future 
state of the atmosphere, because the high- 
speed computer may provide a means — 
hitherto unavailable — for helping to develop 
such methods into an effective quantitative 
implement. 

At present the theoretical foundation of 
modern “dynamic” meteorology consists large- 
ly of classical hydrodynamics, amended let 
us say, as by Reynolds, in particular to include 
eddy transport phenomena. We cannot be 
certain that this framework will suffice ulti- 
mately to account for all of the complex 
phenomena observed in the atmosphere. To 
some extent it is reasonable to hold the view 
that modifications or reinterpretations of the 
classical framework will be suggested in the 
course of systematic scrutiny of the observa- 
tions now or eventually at our disposal. How- 
ever, no matter how thorough or inspired 


GEORGE W. PLATZMAN 


our interpretation of the observations may be, 
we shall not completely understand the limi- 
tations of the classical framework unless we 
learn a great deal about precisely what can 
be predicted within this framework; this of 
course means integration of the (nonlinear) 
hydrodynamical equations and use of high- 
speed computers. 

In this connection, we should mention the 
recent and very promising revival of experi- 
mental methods of approaching atmospheric 
circulation problems by means of hydrodyna- 
mical models studied in the laboratory (notably 
the work of Fultz and his coworkers at the 
University of Chicago). The high-speed 
computer may prove effective in studying 
these controlled experiments theoretically, 
and if so, this coordinated approach may 
establish a more rigorous basis for testing 
the theoretical framework of our ideas about 
the atmosphere. 
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Abstract 


A method is given by means of which quasi-permanent ‘‘velocity’’ fields: may be 
constructed in which absolute vorticity is conserved if the atmosphere is treated as 
barotropic. For vorticity displacements these fields may be considered as constant for time 
intervals of the order of magnitude 24 hrs. Some simple approximate integration formulae 
are given for the Poisson- and Helmholz-equations. A graphical procedure is introduced 
to derive the isolines of vorticity and to perform the integration work. 


i. Introduction 


In numerical weather prognosis one works 
with a certain system of partial differential 
equations. These equations may differ con- 
siderably in form in accordance with how the 
problem has been formulated _ physically. 
Apart from this essentially physical aspect of 
the forecasting problem there is one mathe- 
matical and one practical aspect. As to the 
first one a number of mathematical problems 
crise in connection with an attempt of solving 
the equations. Certainly one is forced to solve 
the equations numerically, i.e. to replace in 
some way differentials with finite differences. 
How this has to be done must to some exteut 
depend upon the physical assumptions waich 
are made. However, even for a given system 
of equations different ways of numerical 
procedures may be tried, and it is not altogether 
certain which one should be preferred from a 
mathematical point of view. The partial 
differential equations in question are namely so 
complicated that a pure mathematical discus- 
sion of how well a numerical solution approxi- 


mates the true one will have a limited value 
only. From a practical point of view, how- 
ever, this may not present any great difficulty. 
The relative goodness of different numerical 
approaches may namely be tested from a com- 
parison between the different solutions with 
respect to how well they approximate the 
actual development (integration) having taken 
place in the atmosphere. 

The practical aspect of the numerical weather 
forecasting problem is related to the time 
necessary for carrying out a numerical solution. 
For it is quite clear that a forecast, how good 
it else may be, is completely useless in daily 
weather forecasting if it cannot be completed 
ahead of the weather. From that reason one 
may prefer one numerical method for an- 
other even on the cost of accuracy if the 
computation time is reduced essentially. This 
will especially be the case if the additional 
errors introduced thereby are well within 
the limits of errors present tor other rea- 
sons. 
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In the following it is dealt with a particular 
numerical approach for solving the so-called 
barotropic vorticity equation. The chief 
achievements of this method may be summariz- 
ed as follows: 


1. For an area approximately as large as the 
one covered by the map in fig. 2 a forecast 
for 24 hrs. can be computed in the course 
of 2—3 hrs., reckoned from the conclusion of 
the analysis of the soo mb map. 


2. The work can be made entirely by one 
man and exclusively with the aid of standard 
means present at any weather forecasting 
section. 


Recalling the integration time for the same 
problem on high-speed electronic computers, 
especially if the necessary introductory readings 
of data from the map are included, one would 
not expect particularly high integration accuracy 
from such a method. However, from a total 
of about 45 forecasts, including the four 
Princeton-forecasts (CHARNEY, FJGRTOFT, VON 
NEUMANN, 1950) it has turned out as a very 
pleasant fact that the integration accuracy 
appearingly is of the same order of magnitude 
as the one obtained on the electronic computer 
used for the Princeton-forecasts. It will be a 
subject for further work to make a systematic 
investigation of the integration accuracy. Such 
an investigation will be much easified after the 
completion of more forecasts based upon 
different numerical approaches. In spite of the 
fact that no systematic investigation of the 
integration accuracy has yet been carried 
through the author will dare to write down as 
a third achievement of the particular method 
in question: 

3. In spite of the short integration time and 
in spite of the simple means by which the 
integration is carried through, the solutions 
obtained seem to be good approximations to 
the solutions of the barotropic vorticity 
equation. | 

In this article weight will be put on ex- 
plaining the particular features of the method 
in question which are responsible for the 
surprisingly short integration time. 

It is hardly worth while mentioning that it 
is the author’s hope that weather forecasting 
centers will start using the method. This 1s 
not because of an opinion that forecasts based 
upon a barotropic model always will give 
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useful results. Certainly it will be necessary 
later to include effects from the changes in the 
mass-field of the atmosphere. However, pro- 
vided the numerical solutions may be con- 
sidered as useful approximations to baro- 
tropic developments the meteorologists will 
themselves be able to encircle all the cases 
where barotropic forecasts fail and thereby 
gain valuable experience with regard to the 
importance of additional physical factors and 
how they possibly could be incorporated into 
the whole complex of computation and argu- 
ments building up a forecast. In this connection 
it should be mentioned that the numerical 
approach presented in this paper also will be 
valuable when applied to more extended 
physical models. This will be seen already in 
this paper from some arguments concerning 
the solution of the so-called Helmholz-equa- 
tion \7?« — a:x= h which appears in more 
extended models. Otherwise, this problem of 
using the particular numerical approach which 
is discussed in this article for less specialized 
models will be discussed more throughly in an 
article to appear later. 


2. The barotropic problem 


In the barotropic model it is assumed that 
the motion at levels close to the soo mb level 
is horizontal and non-divergent thereby con- 
serving the vertical component of absolute 
vorticity by the additional assumption of con- 
servative forces. To reveal most clearly the 
essential features of the particular numerical 
method which is discussed in this article it 
will at first be assumed that the motion is a 
plane one. The fundamental property of 
periodicity around the earth for any quantity 
a derived from some atmospheric variable may 
in the plane representation be expressed by 
writing the quantity as a Fourier series. Since 
we are always dealing with numerical ap- 
proaches to the solutions we will approximate 
the quantity in question by a finite Fourier 
series 


(IN) ani(~ + 7) 
m Am ®, 
(1) X— ZZ Xm,ne REC 
where [yj now indicates the shortest wave- 
length occuring in the series and c is the 
constant in the series. 


INTEGRATION OF THE BAROTROPIC VORTICITY EQUATION 


The law of individual conservation of 
vorticity may be written 


(2) T=—v vi, 


n symbolizing vorticity, v velocity and V7 
the del-operator in the plane of motion. 
Because of the assumption of non-divergence 
for the velocity field, velocity may be ex- 
pressed by a stream function y, v= — Vy x k, 
k denoting a unit vector perpendicular to the 
plane of motion. Vorticity is then given by 
Vy =. The problem is now to solve (2) 
under the condition of periodicity expressed 
in (1). To get numerical solutions one will 
have to replace \7?y by an approximation, 
the corresponding finite difference expression. 
One may write 


(3) a 


where y is defined from 


4) B= Tee thy) vl d, y) 4 
+v(uytd) tyra)! 


One may remark that even if \7?p is com- 
puted from the approximation (3) we do not 
necessarily require the solutions for gridpoints 
only. Even if a numerical solution cannot from 
evident reasons yield a solution everywhere we 
will to begin with in theory require a solution 


valid for any x and y. If now 7 = — + (y—y) 


1 ! : 4 
is substituted into (2) and the factor — 7, 
divided out of the equation one obtains as a 
new continuous problem to be solved under 


condition (1): 


an 
(s) DE im v- Vn 
where 
n=Y—Y. 


The numerical problem of solving (s) for 
a sufficiently short time At may be divided 
into the two independent steps: 


I8I 


a) To displace properly the initial isolines 
n = const over a time interval At to find 
the new position of these lines and by 
subtraction the increment An up to time At. 

b) To solve thereafter Ay — Ay = An to 
find the corresponding increment in the stream- 
function of the velocity field. 


These two steps will be discussed separately. 
At first, however, we shall in the next section 
discuss more closely the importance of the 
operation defined in (4) and some consequences 
which may be drawn from this operation. 


3. Some consequences of the operation ~_ 


Performing the operation L times on each 
term of the Fourier representation (1) one 
obtains 


el sales 
(6) = dee Amin € (in on) Le 
with 
I lc, 44% 4d x 
oe. = = (cos aes + cos = a 


From this last formula it may be concluded 
that for all components for which 4d/A,, 
and 4d/w, are not both equal to an even 
number one must have 


Bee 


The operation means therefore particularly 
a damping of all components with wave- 
lengths Am, w, > 2 d. The damping is complete 
for À, ©, = 4 d, and is from here on decreasing 
asymptotically to no damping for A,,,@, = 2 d, 
and to a reduced damping when 4, wn 
approach their highest values A, and w,, 
respectively. By assuming Ay, @, > 2 d one 
obtains therefore by repeated ~~ operations in 
the limit: 


re 


(8) ee c, provided 
L—>o 


Am, On > 2 d. 
In stata 


The expression for «— x will in a finite 
Fourier series be given by 


MU ed 
(9) ee = SD De ae 
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> Am. Un 


a ET 
Fig. 1. The curve illustrates the goodness of the 
approximation (3) for different wavelengths. 


where B,, , is given from 


AU Reid 
(10) Ban = (sin ey sh sin? an =) Lim, n - 
d2 
On the other hand—— 72% will as obtained 


from differentiation of (1) be given by 


A x Y 
=) ge sey 


The ratio R,, between corresponding coef- 
ficients in (9) and (r1) gives us then a measure 
for how well the approximation V?x = — 


= (x— à) is valid for the different compo- 


nents. We obtain R,,,, from 


PRIT ONE Saad oe 
Sint ee sind Sy 
de 4 On 4 


2 2 
d? 307] Am “aie dr], 


(12) 


Re = 


A curve is drawn in fig. 1 to illustrate the 
variation of Ry,„ for similar variations of 
An and w,. As would have been expected 
Rinn< 1, and is the closer: to zero the 
smaller are the wave-lenghts. It is worth 
while mentioning, however, that it is the 
relative variation of R,,., and not the value 
of Run itself that introduces errors in the 
vorticity equation 
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when 7 there is replaced by — À (y — y) to 
give 


2 (y — y) 


er v-V(y—#). 


The last equation may namely always be 
brought to give complete identity between 
the two equations for à single arbitrary com- 
ponent by multiplication with a suitable 
constant. The relative variation in R,,,1s seen 
to be slight in the range of wavelengths 
Am, On > 3 d, but some approximation to \/?a 
is obtained even down to Am, @, = 2d. 
In the following we shall fix the range of 
wavelengths in which we may expect useful 
solutions of (2) from the solutions of (5) to 
the range 


hp Oped. 


4. The problem of time-integration 


In Lagrangian terms the time-integration of 
an 
ot 
by repeated displacements of the isolines in 
velocity fields which are considered as constant 
for a sufficiently small time increment Ar. 
This Lagrangian formulation of the time- 
integration problem will make sure that the 
fundamental assumption of conservation of 
vorticity is realized during the computation in 
the sence that no isolines 7 = const then will 
be present at any later time which did not 
already exist in the beginning, nor will any 
isolines be lost.’ 

The assumption of constant velocity fields 
during the partial displacements of the isolines 
of vorticity introduces necessarily in the 
general case an error in the distribution of 7 
at a later time t. One may, however, keep this 
error within arbitrarily small limits by choosing 
a sufficiently small value for At. Decisive for 
how rapid a forecast can be computed for a 


= —v-\/n may in theory be accomplished 


L RAN: 
period T is the ratio —-. In the atmosphere 


Ae 
the actual velocity fields are usually changing 
so rapidly that one is forced to keep At less 


! In this connection it should be recalled that in 
the Princeton-forecasts referred to earlier the vorticity 
was not always very well conserved in this respect. 
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than 3—4 hrs. It is possible, however, to 
increase this time considerably if use is made 
of the following considerations: 
We write 

pews ale 


and correspondingly 


vev Ever 


y), c = const 


pe —»y Cc (y 
ey) 

Substituting v = v* + v** into the vorticity 

equation 


aCe Oe 
dt 


v-V (y—¥) 


one gets the equivalent equation 


em 

(14) 2 2 >= v7 (y —ÿ). 
ot 
The term — v**. \/»7 drops namely out of 
the equation because v** is parallel to the 
isolines 7 = y—y = const and consequently 
ineffective in changing the distribution of vor- 
ticity. One may therefore as well displace the 
equiscalar lines of vorticity in the fictitious 
v*-fields as in the actual velocity fields. This 
will now imply a considerable advantage if 
this fictitious velocity field v* is considerably 
less varying in time than the actual velocity 
field itself. It is quite clear that if the p*-field 
shall possibly possess this property, the constant 
c in (13) cannot be taken too small, nor can 
it be taken too large. If particularly c is put 
equal to 1 one obtains 


(15) 


Provided a useful displacement field exists it 
would imply a considerable advantage from 
a practical point of view if it could be derived 
from (15) instead of (13) because initially y 
already exists when 7-9, is computed, and 
because at later time-steps the solution of 
Av— Ay = An gives us immediately the 
change in the y-field. The question then arises 
what size should be taken for the distance d 
involved in the definition of y to make y 


vr 
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most useful as a displacement field. Obviouslv, 
since y—> y when d—>o only slight advan- 
tage could possibly be gained if d were taken 
too small. Recalling the consequences of the 
operation it is clear that in order to make 
y most suitable as a displacement field d 
should be taken equal to the order of magnitude 
of the wavelengths of the most rapidly 
changing components of the p-field. Applied 
to the atmosphere the value 


4d 2,500 km 


suggests itself as a natural value in this con- 
nection. In the forecasts referred to in the 
introduction d was put equal to the length 
of 12 degrees long. at 60 degrees N. It then 
turned out to be possible to take At = 24 hrs. 
Appearingly this only introduced negligible 
errors in most cases in the forecasted position 
of the equiscalar lines of vorticity. The map in 
fig. 2 illustrates a case of computed changes 
in An A(z—2), where z is the height 
of the soo mb surface. The isolines A (z— 7) 
= constant are drawn for every 40 m. The map 
in fig. 3 illustrates the corresponding observed 
changes. As it may be seen the observed and 
computed changes generally compare well. It 
should be remarked that the presented case is 
not one which is particularly favourable to 
the method of displacing described in this 
article because it was a case with relatively 
strong changes in the mean field itself. It 
should also be remarked that only relative 
vorticities were displaced, and also that dis- 
crepancies between observed and computed 
vorticity changes when they are of a disorderly 
character always have a tendency to disappear 
when the integration of Ay — Ay = /\n have 
transformed the vorticity changes into cor- 
responding changes in the velocity field. 


It should now be noted that having found a 
distance d involved in the definition of y 
which makes y most suitable as a displacement 
field, this distance may be too large for a 
sufficiently accurate determination of the vor- 
ticities from the approximation \/? y = — 


— + (y— y). Obviously if the distance d in 


this connection has to be taken much smaller, 
for instance equal to a dy < d, we have to 
return to the expression in (13) to get if possible 
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Fig. 2. November 25, 1951, 0300 GMT: The lines represent equiscalar lines for the computed 24-hr changes in z — 2, 


A ee , and are drawn for every 40 m. The arrows point in the directions — V A(z 


a useful displacement field. We shall abstain 
from a discussion of what value should then 
at best be taken for the constant c in (13). 
Recalling namely what was pointed out at the 
end of the preceding section the approximation 
(3) will only be applicable for components 
with wavelengths down to approximately 
3 d, which with the chosen value of 4 dx 2,400 
km means down to wavelengths of the order 
of magnitude 1,800 km. Phenomena on a 
considerably smaller scale than this cannot 
easily be observed in the free atmosphere 
with any appreciable accuracy: It might also 
be that such phenomena, at least for relative 
short periods, do not have any noteable influence 
at all for the development of the more large- 
scale components of the motion. In that case 
there would also be a physical reason for 
abstaining from an attempt to include compo- 
nents with wavelengths considerably less than 
approximately 1,800 km. 


—2) comp. x k. 


5. Thenumerical solution of the Helmholz- 
equation \7?2«— a « = —h. 


When the changes /\77 in vorticity are found 
from the displacement of the isolines 7 = const, 
one is left with the problem of solving 


p— A = An. 


This equation may be considered as a special 
case of 


(16) 


which, with M and h determined from 


a—Ma=h 


I 
(17) oo GA 
Lcd rente 
4 
and 
deh 
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Fig. 3. November 25, 1951, 0300 GMT: The lines represent equiscalar lines for the observed 24-hr changes 
in 2— 2, A (z—2),,,, and are drawn for every 40 m. The arrows point in the directions —V A (—2), Xk. 


is the finite difference equivalent to the Helm- 
holz-equation 


V2a—a@a =—h. 
The following identity exists: 


(18) «= («—Ma)+M(a—Ma)+...+ 


where 


(19) 


We are looking for a solution of (16) where 
the required quantity « is expressed as a 
finite Fourier series 


(IN) ani = +2) 
Ca 2, Xm,n € ig 2 


+c 


with In > 3 da 
The last restriction on /N results from the fact 
mentioned at the end of section 1 that the 


approximation 2% & — 5 («— a) does not 


apply to the components with wavelengths 
considerably smaller than 3 d We may now 
apply the result of (8) in section 1 in (19). 
We then obtain: 


(20) Ri—>0, a+zo; M<ı. 
L—>o 

and 

(21): Rr—>e, es Ons MSs 
L—+> oo 


We make further the assumption that the 


. 7. I 
relative variation of M = er as compar- 
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ed with the one of & is sufficiently small 
to secure the approximate validity of 
aoe TA I 


apes 


(22) 


if a? is not already a constant. In the case 
of a Poisson equation, a? = 0, and (22) is 
authomatically satisfied. The identity (18) may 
therefore by applying (22) be written 


(23) «= («— Ma) + M(x— Ma) + ...- 
\L 
ben Cr 


However, «— Ma is equal to the known 
quantity h in (16). We may therefore also 
write (23) as 

E93 

= Rey 

CG = h EE Mia ..... MO HER; 
Applying now (20), (21), the solution for « 
may be written as the infinite series 


=} : 

(25) a BANE rare OE ip ede 
+c 

where c=o, when a? + o 

and c = arbitrary, when a? =o. 


The arbitrary constant occurring in the case 
of a Poisson-equation is of course of no 
physical importance, as we, for instance, may 
understand when & is equal to the stream 
function. In the case of a Poisson-equation 
x — x = h it is immediately understood that 
the function h cannot contain any constant 
in the Fourier-development. This may, how- 
ever, be the case when we have a Helmholz- 
equation «—Ma-=h. The solution for « 
will in this case contain a constant determined 
from the solution (24). By the application of 
the Helmholz-equation in the more extended 
computational models of the atmosphere « is 
the temperature-change. From physical reasons, 
now, it may easily be understood why tem- 
perature in the mean may increase or decrease 
in a horizontal level extending all around 
earth. 
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The question about the rapidity of the con- 
vergence for each Fourier component of the 
infinite series (24) may be discussed easily by 
means of the expression (19) for the remainer 
Rr in (24). 


Using (6) one obtains 


Fe a 
(26) Ren (Rr)msn e 
Up) IE: zai (~ + 7) 


= > Xm,n Gra à M) e 


where now Ayn is given from 


It is easily understood that the convergence is 
most rapid for the component with 4,, ©, = 4d 
in which case A,,,= 0. From there on, how- 
ever, the rapidity of the convergence is de- 
creasing for Am,@, decreasing downwards to 
3 d which is the lower limit of the range of 
wavelengths we are interested in. The rapidity 
of the convergence is also decreasing steadily 
for Am, #, increasing upwards from 4d. As 
it may be concluded from (26) and (7) one 
must for instance take five terms in the in- 
finite series (25) to get a solution in the Poisson 
case which is correct to 25 % in the whole 
range 4d < Am, ©, < 8d. However, one can with 
fewer terms obtain solutions with approximate 
validity that approximate the true solutions 
more equally well over a wider range of 
wavelengths. The argument is as follows, 
taking for sake of illustration the Poisson case 
with M= 1: 


Let now 7 denote any of the operations 


= Sh 2 


2 Aug 
eer | defined from 


ee AG (x + du, y) + B(x—d,, y) + 


+B ley +4) + By — di) 


where 


din. Berge bel 
. = 0 
Let us particularly define 6 from 


B =B6. 
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We may then put down the following identity: 


= =, 5 a+tı 
(27) DA 4 (ae fe) + FL 
— 0 =0 
DH : | 
where 
—L+1 
— I 
(28) e =a 


Corresponding to the result (6), (7), in section 
3 one must have 


un tas 
it ey 
3 270 (— ern 
4 @ À 22 
(29) aa == ( | 40) Xm,n € om é 


g=1 
where 


(30) 


(q) I 4d. TT 
a, = er (co = 1 | 


For comparison the remainer Rı in (26) for 


the Poisson case with M = 1 is: 
UN) L+1 aai(~ ++) 
A @ 
Rı == >» CE) Xi, n € £ 


While Rı is zero only for An, On =4d, rı 
Peiseci CO" be- zero fOr A, On = Ad Ada, +, 
4dr..1t would therefore so far bean advantage 
to use (27) instead of (18) as a base for obtaining 
approximate solutions of the Poisson equation 
« — à =o. It will then, however, further be 
necessary to express at least approximately 
the general term in the identity (27), after the 
idee has been left out, in terms of the 
known function h. Denoting this general term 
with 


5) ©. 

P,= ‘a? —\a- q=0, I L 

and letting further Q, be defined from 
Sm) 

Or= er — an Fl eam a AB 


one obtains easily for the ratio between cor- 
responding Fourier components of P, and Q,: 
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(a+) 
| \ == 
\ PB; fn I Ar: 


Q, \ m,n I — A (a) 


m,n 


(31) = Be 


Hence by substitution for oe and ey from 
(30), one obtains 


ane 4dgtı 7 9 Adarı 7 
[ 
Am 4 On 4 
ae d 4d 7 
sin? —4 — + sin? —4 
Am On 


With the notation 


Kq >= 


Run will for similar variations of An, @n, say 
his On = AVE tas 


TU I Ess 

sine eee 

WEN 
k, = a 

PTE 

sin? 22 

À 4 

It is seen that 
4—> co 


It is also understood that k, will be but slightly 
varying for large values of 2. If one wants 
k, as a quasi-constant for 


1e, On > 4 dy, 
one must take 


HMS: approximately. 
In the following we shall take 


(32) 


We may then replace k, with the constant 2 
which is approximately the mean value of 


bia and ko 


a Bi 


Bun N 2, ns On = 4d, : 


Consequently we may write (31) as: 


dn) ani (= +2) 
“m on 
20 =P, = Sir Ve ne 


(SEP. 
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(34) A Ont Ade: 


We note from the definitions of Q, and P, 
that 
Sl 


(35) Q,= Pr 
Combining (33) and (35) one therefore gets: 
PAPE, Figs dE dd 


By repeated use of this relation and condition 
(34) one finally obtains in connection with (33): 


er 
(36) PAPA, 
provided 
da On = 4 d, 
Using 


| q 
> (r) 
{Poln.a = ( | | 4°) Palme 


T=1 


in connection with (36), we may write: 


q 
fp*\ =2'(| | 49 {pr 
ieee ( ee 


r=I 


or, when similar variations of A,,, @, are con- 


sidered: 


toe 4da,x 4d, x 4d, x 
VP ja 2" cos u rang M ‚08 7X 
x 1 Pet min: 

Consequently P; is zero TOR Any On dd), 


44, tag tA des and iswalso Small, for 4,40, 
having values between the values of 4d,, 
4d, ...,4d,, if q is not taken to large. This 

r 
Opa 
been chosen so small. 

Making use of this last fact one may by 
substitution from (36) into the identity (27) 
write this approximately as: 


is because the ratios 


[ED thes Glave 


—L 


2 ur 
(37) a=P, +2 Bele gO Py = ry 
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for all 
4 d, = Ans On 


if L is not taken too large. The relation (36) 
exists namely in parts of the range of wave- 
lengths which are considered, and in the 
remaining parts where it does not exist the 
terms in consideration are negligible anyway 
in comparison with other more dominating 
terms. The remainder rz given from (29), (30) 
will be small in the whole interval 


4d, = Anns On = 4dL+:- 


By disregarding the remainder in (37) above 
and by substituting 


ee h 


one should therefore get an approximate 
a 5 é ak 
solution « of the Poisson-equation «— x —has: 


—L 


(38) et 


provided 
4d, < A On < 4dıı 1° 


The accuracy of the solution may be examined 
for the different Fourier components in the 
following way: 


We have 


Bun D (a ine; ae == (1 Fuge AN.) Kın,n- 


If this is substituted into (38) one obtains for 
each Fourier component: 


* 
Xn, (1) CO. 
(39) = (: = Ann) > (: Sr 2 Ann .. + 
L (x) (L) 
ee tin ... An) : 


When this ratio is + 1 the solution is correct 
for the corresponding component. In the 
tabulation below this ratio is given as a function 
Of Am, @, = À measured in the unit 4 4, 
and for L assuming successively the values 
142,0 . 
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Table I. 
4 di 5 da, 6d, 7 dy 8 di MIROIR Ne | 87 
* ! 
En, nl m,n? UT M ce 1.00 1.212 1.00 .85 .82 ST .62 .54 .36 = 
* 
aa: Ve ra 1.00 .82 97 1.04 1.00 .92 .85 77 .67 — 
* 
RER VE  E 1.00 1.40 1.02 .85 .90 .89 .89 .88 .87 .76 


It is seen that as approximate solutions of 
Ber 


&—x = h we may use the following ones: 


(40) dr = h + ah, ANS AS On OG, 


id = 1.002 100, 


=o 
—2 —2 


+ ae ——I ae! 
M2) ch +22 Euh +8h , 
AAS hin On 164. 


If one wishes a solution with only two terms 
as in (40) which is better for wavelengths 
around 8 d, one could use a solution: 


(43) 


The goodness of this approximate solution is 
measured by 


eo = h + ah. 


x (1) (1) 
a= (:— An) k (: os i 3 ay, 
Im,n 


which is tabulated below: 


Table II. 
4dı|5d,|6d, | 7dı| 8d, | 94, |10d, 
KX] 
Sn, nl mn: I RES | 2:34 a s | 09 


To demonstrate the usefulness already of 
the approximation (43) a map for the 24 hr. 
observed changes (/\2)ws in the height of the 
soo mb surface, fig. 4, was used to compute 
corresponding “Laplacians” (Az— Az) = h 
which was then in turn substituted into 
formula (42) to give a set of integrated values 
PAZ)", fig..52 Assisseen the correspondence 
is surprisingly good in spite of the fact that 


the case was not a particularly favourable one 
to the integration method. 


It is worth while mentioning that in the 
general case of a Helmholz-equation « — Max = 
— h, the convergence for larger wavelengths 
may be speeded up for the series (24) by pro- 
ceeding in a similar way as was done for the 
Poisson case. But in the case of a Helmholz 
equation this is not as necessary as in the case 
with a Poisson equation because the conver- 
gence for larger wavelengths have a tendency 
to be rapid already because of the multiplica- 


tion with M**' in the remainer (26). 


6. Modifications due to the curvature of 
the earth. — Geostrophic approximation 


The curvature of the earth modifies the 
results of the preceding section physically and 
geometrically. Let now \/; denote the spherical 
del-operator and k a vertical unit vector. 
Relative velocity v, and relative vorticity 7 
will then be expressed by a stream function y 
as follows: 


v=—V;yxk 
and 

n= Ver. 
If 

f=-202 sin çp 


where 2 is the angular velocity of the earth 
and @ is the latitude, the vertical component of 
absolute voticity Na will be given by: 


Nabs = Ve LA +f. 


The conservation of absolute vorticity may 
now be expressed by: 


I Nabs 
(44) 3 


Zu ni V sNabs- 
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Fig. 4. November 25, 1951, 0300 GMT: The lines represent equiscalar lines for the observed 24-hr changes in 
ZEN ZUR, and are drawn for every 40 m. The arrows point in the directions — V a SS k. 


To find a useful displacement field one since the term v**- \7s 7a drops out because 


may proceed in exactly the same way as in v*™ is parallel to the isolines of absolute 
section 4. We then write vorticity and therefore ineffective in changing 
PR: re the vorticity distribution. Eq. (45) above 
D Un therefore expresses the conservation of absolute 
b vorticity in the field 
and correspondingly ; 
N° = — V5 {pate (VE y + f}\ xk. 


v = v* + NL 
We now introduce the geostrophic approxima- 


where fone 


_ gVszxk 


and 1 
F Hence relative vorticity will be given as 
TA ur D) 


p* = p—c(Viy +} cæconst 


y= 


Viv =F Viz + VE xfv-k 


or approximately as 


Substituting v = v* + v** into (44) one ob- 
tains 


; 4 IN abs 


(45) St 7 Pan (46) Vive Viz 
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30 I 30 20 30 
Fig. s. November 25, 1951, 0300 GMT: The lines represent equiscalar lines for the 24-hr changes in z, Nee 
found from the approximate integration formula Az,, = A CE NE NET) a 


The arrows point in the directions — V Az. 


Sur 
since — may be considered as a virtual constant 


under the curl-operation V,X v : k. With 
(46) the vorticity-equation (44) can be written 


rev (E Vie +f), 


Working in a conformal map projection with 
the magnification factor m, one obtains 


Viz ES AN trs 


Here V7 denotes the plane del-operator. Sub- 
stituting the finite difference approximation 


Che. (z— 2), d = const. 


in the last equation one obtains after division 


Sake 


int. 


with the factor — a the finite difference equi- 


valent of (47): 
m 9 (z— 2) v: ale 


ee ac ua gas 


It may be supposed that the relative variation 
Uh 
in Te 
variation in the quantity 2— 2. The last 
equation may therefore also be written with 


sufficient accuracy: 


is small compared with the relative 


RAC nm VU TRE € 
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N mag 
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Let us suppose that m is a function of latitude 
only. Then the last term in (48) may be 
written as 


m? 40? d? sin p cos y 
V . = 


f m= (p) 48 


V:® 


or, with the notation /(p) defined from 


(42? d sing cos p dp 
(4s) J@)= [Sate ag 


as 
m2 


se VsJ(¢)- 


Substituting now this for the last term in 
(48) one obtains after dividing out the factor 
m? 


0Ë 
(50) Aa Te Vet) 


(s1) §=z—2z—J(¢). 


This equation, then, expresses the individual 
conservation of the quantity £ = z—z— ](q) 
in the geostrophic wind field v. 


where 


To get an useful displacement field we now 
simply introduce in line with earlier con- 
siderations instead of v a fictitious velocity v* 
determined from: 


* 8 * 
win AT xk 
4) 


z z 2 — z 


Jo), € = const. 


Choosing the constant equal to unity one 
gets: 


z"= z+ J(y) 
and correspondingly 


with 
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Instead of (50) one obtains now the equivalent 
equation 
dé 


(53) ae ee eee 


This equation now expresses the conservation 
of the quantity & in the field v + cu. ca as 
obtained from (52) and (49) is seen to represent 
a velocity directed westwards and is in fact 
the velocity of propagation of a so-called 
Rossby-wave with wavelengths approximately 
equal to 4 d. Since d has been taken to be 
only equal to the length of 12 degr. long. at 
60 degr. N this wavespeed will be a small one 
and in most cases only imply a small correc- 
tion to the field v. The main contribution from 
the variability in the Coriolis parameter will 
result from the term + v- \7;J(~) when 
the smoothed fields v have noteable compo- 
nents in meridional direction. 


7. Practical operation in a weather 
forecasting section 


To accomplish a forecast in accordance 
with the principles outlined in this article it 
has from the following reasons turned out to 
be very convenient to make extended use of 
graphical additions and subtractions, rather 
than reading data in selected gridpoints only 
and furnish the necessary computational work 
with these data: Firstly, the analyzed map is in 
a sence utilized everywhere when graphical 
methods are used, whereas in a grid the 
gridpoints may have a very unfavourable 
position for a sufficiently accurate determina- 
tion of the required quantities. Secondly, the 
graphical method is much faster than the grid 
method. For instance the whole time for the 
computation of a map for the isolines z— z= 
const reckoned from the conclusion of the 
analysis of the s00 mb map amounts to only 
30—45 min, the area of the map being ap- 
proximately as large as the ones covered by 
the maps presented in this article. Thirdly, 
one gets in the graphical method the isolines 
for the required quantities directly. 

At a weather forecasting centre, the actual 
operational procedure then has to be as 
follows, or approximately as follows: 


I: Construction of the initial field £ = z — 
— 2 — J(y) = const. 
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II: Displacement of the isolines &= const in 
the displacement field v + cy, to find 
* the new position of the lines after 24 hrs. 


II: Solving the equation A z— Az=fé 
to find the changes in the height z of the 
soo mb map after 24 hrs. 


I. 


The soo mb. map, here denoted as the 
z-map is copied over on a blank map dis- 
placed in an arbitrary direction x the distance 
d — equal to the length of 12 degr. long. 
at 60 degr. N. Only every second line is 
copied. This gives us a z.:4-map for the 
isolines 2, +4=!/g32(x + d, y) = const. There- 
after is the z-map copied over once more 
on to the zy+y-map, but now displaced a 
distance —d in the x-direction. Again only 
every second isoline is copied. This gives us 
a map for the isolines zx_4= const. By 
graphical addition of the z,+-map and the 
Zx--a-map one obtains when only every 
second line is drawn, a map representing the 


field 


I 


A = - {arts — 2-4) —const 
4 


By repeating this for the direction perpen- 
dicular to the chosen x-direction one obtains 
similarly a representation for the isolines 


„ul 


I | é 
= <2ytd I <y—dy — Const. 
4 | } Y | 


Then by the graphical addition of the zl-map 


and the 2!-map above one obtains 


sw 


N 
N | 


Then finally from a graphical subtraction of 
the Z-field from the field z—7J(p) one 
obtains the initial distribution of € given as 
isolines 


E = z—z—J(¢) = const. 


It will obviously be most practical to have the 
isolines J(g) = const printed on blank maps 
once for all. 


. 


4—203045 


193 
H: 


The isolines & = const are displaced in the 
field v and then afterwards corrected for the 
relatively small westwards displacement due to 
the velocity cy. A wind-ruler of the ordinary 
hyperbolic type may be used. Instead of 
beginning with the displacement of the isolines 
& = const, one may at first draw a suitable 
line crossing all the isolines z = const and 
then compute the displacement of this line 
for Le ER ti 2/\t,..., untilthe whole 
map is sufficiently densely covered of such 
lines. The displacement of the isolines &= const 
will then afterwards be relatively easy. 


Ul. 


After having got the new position of the 
&-lines one obtains A Ë£ from the graphical 
subtraction 


Son var So = N £ . 
For integration of the equation 


Az—Az=AE 


one may use either of the approximate solu- 
tions (40), (43) or (41). Again, only graphical 
additions are involved in the computations. 


8. Final remarks 


In some future results obtained from the 
integration method in this article will be 
presented more systematically, also its exten- 
sion to more extended models where absolute 
vorticity is not conserved. When the 45 fore- 
casts referred to in the introduction were 
computed the integration formula (40) was 
used. It has later turned out that an use of the 
integration formula (41) instead, or even 
of (43), would have improved the integration 
results considerably in: many cases. Also in 
most of the forecasts the influence from the 
variation of the Coriolis parameter was com- 
pletely disregarded. It would therefore not be 
fair to the method to publish actual results at 
the present time. It will particularly be 
interesting to see how the forecasts compare 
with the ones obtained with electronic com- 
puters. 
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The time necessary for carrying out the 
three steps above amounts to between 2 and 
3 hrs. After considerable experience one should 
expect that it will be possible to bring the 
time down to 2 hrs., particularly if one could 
use carbonpaper for the copying of the maps 
involved in the integration work. Having 
brought the time down to 2 hrs one could of 
course take the time and use systematically 
the more accurate integration formula (41). 
One could also possibly forecast for 48 hrs. 
if the barotropic model applies to such time 
intervals. 
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Abstract 


Methods now in common use for making prognostic charts to a large extent are based on 
time derivatives, while the numerical methods start from space derivatives. That is the funda- 
mental difference between them. It is shown that the present methods give at times very good 
results, and that the degree of success varies with the weather situation, characterized e. g. by 
the average interdiurnal pressure change over the test area. Because of this variation it would be 
desirable that every test of a numerical method includes a comparison with results obtained at 
some forecasting center or, at least, with the results of pure extrapolation. 


Introduction 


At present, in forecasting the sea-level pressure 
pattern as well as the upper flow pattern, 
in particular methods of extrapolation are used, 
i.e. the forecaster knows the previous history 
of the flow pattern, and he studies the changes 
or tendencies, the time derivatives. He may 
then introduce, in addition, certain correc- 
tions to the extrapolated pattern, argued from 
space derivatives, such as advection, diver- 
gence, vorticity etc., but extrapolation will 
always be the first step. It may be the 
areas of 24-hour changes which are moved 
either in accordance with their previous mo- 
vement or along the track of the preceding 
changes, or it may be the troughs and ridges 
that are moved in a similar way. By extra- 
polation he can also make ridges build up or 
weaken and make troughs deepen or fill. 
Such a graphical extrapolation is rather easy 
to carry out, and it is not very subjective. 


In general, however, extrapolation is not 
enough. The forecaster also must try to find 
signs of real evolutions, and then the methods 
grow more and more subjective. He has to 
look for areas of confluence or other kinds of 
advection of warm and cold air, and he has 
to estimate the effects to be expected from in- 
teractions between patterns over different parts 
of the map. He also has to apply rules of ex- 
perience, general rules, rules valid for a cer- 
tain season as well as special rules, just applicable 
to a prevailing type of weather. 

But we should keep in mind that all these 
more subjective steps are just refinements, 
and that pure extrapolation already can give 
surprisingly good results. 

In contrast to this the numerical methods, 
actualized by the advent of high-speed elec- 
tronic computing machines, are all based on 
space derivatives, and the time derivatives 
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‘and the previous events do not enter at all. 
But at a meeting on numerical forecasting held 
in Stockholm in May 1952, there was briefly 
discussed if the time derivatives could not be 
entered in some way in the numerical predic- 
tion, €. g. to improve the analysis of the 
vorticity pattern in areas with sparse observa- 
tions. Also the possibility of taking the previ- 
ous events into account was mentioned, e. g. 
by introducing corrections in the forecasts 
for “past errors”, as involving effects of 
neglected factors. 

The development of the flow pattern is in 
present theories fully determined by the initial 
data, including all sorts of space derivatives. 
As however the time derivatives have proved 
invaluable for the non-automatic forecaster, 
it may be necessary to try to introduce them 
also in the numerical process. This ought to 
be done also if we do not want to replace the 
present methods entirely by numerical ones, 
but restrict ourselves to let the two kinds of 
methods supplement each other. 


Present routines in forecasting 


The goal for the meteorological prediction 
is in general to forecast the weather, i.e. 
surface wind, temperature, cloudiness and pre- 
cipitation. It is true that forecast upper flow 
patterns may be used directly in preparing 
winds for aviation and ballistic purposes, but 
although this side of the problem is impor- 
tant, it will be neglected in the following. 

To obtain a weather forecast we can either 
begin with an upper-air prediction and pass 
on to the sea-level and from there to the 
weather. Such a routine was recommended 
e.g. by RIEHL and collaborators (1951). Or 
we can begin with a sea-level prediction, as 
proposed by SCHERHAG (1948). In the latter 
case we can then pass to the weather almost 
directly, but our forecast will certainly im- 
prove if we construct an upper-air chart 
from the sea-level chart before making the 
weather forecast. 

The advantage of the first way is that an 
acrological chart gives a much simpler picture 
of the weather situation than does the sea- 
level chart. The disadvantage is that upper-air 
charts must be translated into sea-level isobars 
and fronts before forecasting the weather. 

The advantage of the sea-level way to the 
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weather forecast, on the other hand, is partly 
the more direct relation between a sea-level 
chart and the weather, partly the existence 
of 3-hour tendencies on that chart. The 
disadvantage is the fact that sea-level prog- 
nostic charts can hardly be constructed directly 
for a longer time interval than 24 hours. 

The choice between these two alternate 
ways depends on which one gives the best 
forecast in the actual case. Generally a com- 
promise is the best way. 

At the Acrological Division of the Swedish 
meteorological and hydrological institute the 
second method, in a form drawn up by Sm, 
is used for 24-hour forecasts, i.e. prognostic 
charts are prepared for sea-level pressure, rela- 
tive topography and soo-mb topography in 
the said order. For 48-hour forecasts, on the 
other hand, the first method is used. Then the 
already prepared 24-hour s00-mb prognosis 
can serve as a hint. 


Testing prognostic charts 


Forecasts made by numerical methods have 
been tested by computing the coefficient for 
the correlation between computed and ob- 
served changes. Such tests have been perfor- 
med by Bouin & CHARNEY (1951), SAWYER & 
BusxBy (1951) and Starr MEMBERS, Univer- 
sity of Stockholm (1952). It can be questioned 
if the correlation coefficient is the most 
appropriate number to characterize the degree 
of success in forecasting. Other numbers 
could be and have also been chosen for that 
purpose, such as the mean error, the ratio 
between mean error and mean change etc. 
However, we shall accept here the use of 
correlation coefficients, in spite of the fact 
that such a coefficient is not sufficient to 
state if a certain numerical forecast was suc- 
cessful. It is a common suggestion made by 
experienced forecasters that you must know 
in addition the corresponding coefficient for 
the forecast chart made for the same day by 
means of the old methods. This is necessary 
since the degree of difficulty in forecasting 
varies considerably from one situation to 
another. 

When the Swedish routine is used, the 
correlation figure for the first prognosis, 
the 24-hour sea-level chart, will in general be 
higher than the figure for the secondary 24- 
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hour soo-mb prognosis. That will be shown 
in the following (see Table I). When the 
Opposite routine is used, e.g. by numerical 
forecasting, we should then expect that the 
correlation figure for the soo-mb level will 
in general exceed the one for the sea-level 
chart with about the same amount. There- 
fore, if only soo-mb forecasts are compared 
by correlation methods, a considerably higher 
coefficient for the numerical forecast is re- 
quired than for the routine one before we 
can say that the two methods will give just as 
good sea-level charts, which are in fact the 
most important ones for the forecaster. 

As to the last step, the weather fore- 
casting, we do not know exactly to what a 
degree the success in forecasting depends on 
the qualities of the prognostic chart. Hence 
it is not quite sure that improved prognostic 
charts, made in an automatic way, will im- 
prove automatically the weather forecasts. But 
to make the prognostic surface charts we need 
at present a good deal of reasoning, and it is 
quite sure that this reasoning is very impor- 
tant for the success in weather forecasting. 


Correlation coefficients obtained 
by present methods 


For certain periods during the last two years 
coefficients have been calculated at the Swe- 
dish meteorological and hydrological insti- 
tute for the correlation between expected and 
observed 24-hour changes in sea-level pressure. 
In all, 228 prognoses have been tested in that 
manner, and the obtained mean coefficient 
is 0.71. The success varied within wide 
limits. Thus the mean for the 30 cases in 
April 1951 was as high as 0.82. The variation 
may depend on varying ability of the fore- 
casters. The variation due to the weather 
situation is probably more important. This 
effect has been studied in the following way 
(see fig. 1). In each one of the 228 cases, 
the coefficient of correlation, r, has been 
plotted against [A p|, i.e. the 24-hour change 
in sea-level pressure, averaged over the test 
area. This area covered northwestern Europe, 
and the number of points in the test grid was 
26. The grid is shown in fig. 2. 

We see from fig. 1 that forecasts of high 
quality are generally obtained when the ave- 
rage pressure change is large. When on the 
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Fig. 1. Coefficients for the correlation between expected 

and observed pressure changes plotted as a function of 

the average diurnal pressure change over the test area. 

Crosses stand for cases studied in Table I. The curve 
represents average conditions. 


other hand the average pressure change is 
small, the success is highly variable, high corre- 
Jations occur but also very low ones. In situa- 
tions where the changes exceed in average 
to mb per 24 hours, the mean correlation is 
0.85. The figures corresponding to the inter- 
vals 8—10, 6—8, 4—6 and 2—4 mb/24 hours 
are 0.83, 0.73, 0.69 and 0.55, respectively. 
The curve in fig. 1 represents this average 
relation between the correlation and the “in- 
tensity” of the weather situation, given in 
terms of the average diurnal pressure change 
over the test area. 

The general experiences represented by 
fig. 1 are probably well-known, viz. it is 
easier to make a good prediction of the 
changes to occur when large changes are 
going on than in a flat and diffuse situation. 
But it is necessary to remind of the important 
fact that even if a prognostic chart gives a 
low correlation figure in a situation with 
small changes, the chart itself often looks 
fairly correct and is of good help in weather 
forecasting. 

In 8 cases of the 228 ones, correlation 
coefficients were computed also for the 
changes in relative topography between 
1,000 and 500 mb and for the changes in the 
5oo-mb topography. These cases are denoted 
by crosses in fig. 1. The results are given in 
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Fig. 2. Sea-level isobars (solid lines) and 24-hour isallobars (dashed lines) for the actual and prognostic chart of 
Nov. 13 and for the observed chart of Nov. 14, 1950. Contours labelled in mb. Correlation coefficient equal to 
0.92. Test points inserted. 


Table I. We sce from the table that it seems 
to be easier to predict the sea-level pressure 
changes with accuracy than to predict the 
changes in the thickness pattern. We also see 
that there is in general, as mentioned above, a 
decrease in the agreement between expected 
and observed changes when the relative topo- 
graphy is added to the sea-level chart. 


Table 1. Correlation between predicted and observed 
changes in sea-level pressure, relative topography 
and 500 mb topography. 


Date Sea-level | 209-7590! oo mb 
mb 

22—23 Oct. 1950 0.93 0.45 0.56 
23224 09 1950 0.84 0.76 0.84 
2A 25023 1950 0.68 0.81 0.86 
13—14 Nov. 1950 0.92 0.67 0.82 
I4—I5, » 1950 0.88 0.18 0.92 
10—11 Apr. 1951 0.90 0.73 0.75 
TT 127), 771651 0.88 0.79 0.85 
I2—13 » 1951 0.98 0.84 0.78 
Mean | 0.88 | 0.65 | 0.80 


Bold-face types are used for the best figure 
in each column. From the forecast for 12—13 
April we see that the best sea-level forecast 
combined with the most successful relative 
topography forecast resulted in a soo-mb 
forecast which was one of the least successful 
of the forecasts studied here, and from the 
case of 14—15 November is seen that the best 
soo-mb forecast was obtained in spite of a 
failure in predicting the changes in the thick- 


ness pattern. This is peculiar, but we have not 
cases enough to see if it is just an accident. 
Anyhow, it demonstrates the complexity of 
the forecast problem but also, to some degree, 
the imperfection of the correlation method. 
In the case of 14—15 November the sea-level 
prognosis had its largest errors in some peri- 
pheric points in which the predicted changes 
in relative topography were also in error, 
while the agreement in the other points was 
good in the sea-level prognosis but moderate 
in the other one. 

For the 48-hour changes in the soo-mb level 
the correlation coefficient has been computed 
in two cases, only, viz. 13—15 November and 
14—16 November 1950. In both cases the 
figure was 0.79,11. e. a figure but slightly lower 
than the one for the 24-hour changes. 

In Fig. 2, 3 and 4 is given an example of 
combined forecasts for sea-level pressure, 
relative topography and soo-mb topography. 
These forecasts were made on November 13, 
1950, and the correlation coefficients are close 
to the mean values for the eight cases given in 
Table I. From fig. 2 we see that both the mo- 
vement and the deformation of the principal 
fall-area are well predicted, but when the sea- 
level charts themselves are compared, we 
find a remarkable discrepancy as to the inten- 
sity of the low west of Norway, in spite of 
the high correlation 0.92. Regarding the rela- 
tive topography we see from fig. 3 that the 
small changes were correctly sketched as to 
their main features. In fig. 4, giving the soo-mb 
charts, the movement of the pattern of 24- 


COMPARISON BETWEEN DIFFERENT METHODS OF FORECASTING 


Bias a. 


1,000—$00 mb relative topography (solid lines) and 24-hour changes (dashed lines) for the actual and 


prognostic chart of Nov. 13 and for the observed chart of Nov. 14, 1950. Contours labelled in tens of metres. 
Correlation coefficient equal to 0.67. 


Fig. 4. soo mb topography (solid lines) and 24-hour changes 
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(dashed lines) for the actual and prognostic 


chart of Nov. 13 and for the observed chart of Nov. 14, 1950. Contours labelled in tens of metres. Correlation 
coefficient equal to 0.82. 


hour changes turns out to be well estimated, 
giving an almost correct trough position. Over 
the southwestern part of the chart, however, 
there is a significiant discrepancy between the 
expected and observed upper flow pattern. 
Nevertheless the correlation coefficient is as 
high as 0.82. 

Fig. 5 demonstrates an extraordinary good 
prognostic chart, tested in a somewhat different 
way. It is the 48-hour soo-mb forecast made 
December 17, 1951. It is tested in 54 points 
covering an area from Hudson Bay to Russia. 
The correlation coefficient was 0.92. If we 
leave out the nine westernmost points, the 
correlation goes up to 0.96. It is worth men- 
tioning that in this case the situation was such 
that the forecast was made by pure extra- 
polation. Thus pure extrapolation can give a 


high correlation even for 48-hour changes 
over such a wide area, but the extremely 
good result is of course an exception. Neverthe- 
less it demonstrates the well-known fact that 
the difficulty in forecasting meteorological 
charts varies considerably from one day to 
another. In this very case, numerical predic- 
tion must give just as high a correlation to 
be considered acceptable for practical use. In 
other cases, however, simple extrapolation 
gives perhaps not more than 0.20 correlation. 
Then the numerical prediction would be very 
successful even if the correlation is about 0.50. 
Therefore it is desirable that every test of a 
numerical forecast includes a comparison with 
result of pure extrapolation. Only such a 
comparison clearly points out the practical 
qualities of the numerical methods. 
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Fig. 5. A comparison between expected and observed 48-hour changes in $00-mb topography 17—19 Dec. 1951. 
Contours labelled in tens of metres. Correlation coefficient 0.92 computed from test points inserted. 
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Abstract 


Elementary analyses of the mean patterns of free convection from a line source and a 
point source are presented without regard to the specific means by which the gravitational 
action is produced. The derived functional relationships are then verified and completed 
through use of velocity and temperature measurements above sources of heat, the general- 
ized form of the results permitting characteristics of the mean flow to be determined 
over a considerable range of the primary variables. These results should enable meteoro- 
logists to evaluate the role of the basic convective process in the more complex movements 


of the atmosphere. 


Introductory Remarks 


Atmospheric disturbances are generally so 
complex in nature that the relative impor- 
tance of the various factors which they in- 
volve can be appraised only by investigating 
the effect of each factor independently. In 
conducting an analysis of this nature, one 
often notes a close similarity between partic- 
ular aspects of meterological and other flow 
occurrences, indicating that experience gained 
in related fields can be adapted to the problem 
in question. A case in point is the phenomenon 
of large-scale thermal updrafts in the atmos- 
phere, one aspect of which is the comparatively 
simple process of gravitational or “free” con- 
vection from a boundary source, recently 
studied experimentally at the Iowa Institute 
of Hydraulic Research. The writers believe 
that the results of this investigation will per- 
mit meteorologists to evaluate the role of the 
basic convective mechanism in atmospheric 
phenomena. 


Free convection due to a point source of 
heat is very simply illustrated by the plume 
of smoke which rises from a cigarette in 
otherwise stagnant air. Because of the buoyancy 
of the heated air in the immediate vicinity of 
the burning end, a continuous upward current 
is induced, with a corresponding radial inflow 
for reasons of continuity. While the initial 
steadiness of the smoke filament indicates 
purely laminar motion for a considerable 
distance above the heat source, the flow there- 
after becomes unstable and the filament breaks 
up into the eddying clouds normally associated 
with turbulent motion. After the onset of 
turbulence, the effect of the convection upon 
the surrounding fluid becomes far more pro- 
nounced, the molecular shear and heat transfer 
of the initial laminar motion becoming 
dwarfed in scale by the macroscopic mixing 
action. Whereas the magnitude of the veloc- 
ity along the vertical axis is then rapidly di- 
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Fig. 1. Definition sketch. 


minished, the angle of radial diffusion (and 
hence the rate of radial entrainment) is accord- 
ingly increased. 

Although free convection of this nature is 
usually produced by a source of heat, it is 
the gravitational rather than the thermal 
aspects of the flow which are fundamental 
to the phenomenon. Indeed, essentially the 
same convection pattern can be produced in a 
liquid instead of a gas, and without any tem- 
perature change whatever. For example, the 
atmosphere may be replaced by a large body 
of water, and liquid having a different specific 
gravity may be introduced at a constant rate 
at a point on either horizontal boundary — a 
lighter liquid at the bottom or a heavier one 
at the top. In each event a similar pattern of 
convection will result, provided only that the 
local changes in density are relatively small in 
magnitude. 

In the following pages an approximate anal- 
ysis is made of the mean pattern of free con- 
vection from both a line source and a point 
source, without regard to the specific means of 
producing the gravitational action. Experi- 
mental evidence is then introduced to verify 
and complete the derived functional relation- 
ships, and generalized diagrams are presented 
to permit evaluation of the ‘primary flow 
characteristics over a considerable range of 
the parameters involved. 


Theoretical Analysis of the Mean Flow 


If it is assumed that the scale of the initial 
laminar zone is small compared with that of 
the subsequent zone of turbulent convection, 
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notation for either two-dimensional flow 
over a line source or axially symmetric flow 
over a point source may be reduced to that 
shown schematically in Fig. 1. Herein the 
origin O represents the location of either source. 
The local mean component of velocity in 
the vertical direction x is denoted by u, and 
that in the lateral direction y or radial direc- 
tion r is denoted by v. 

In addition to the mean velocity of convec- 
tion above the source, there will also be a 
mean local change Ay in the weight density 
of the fluid. This entails, to be sure, a com- 
parable change in the mass density @. How- 
ever, the following analysis is based upon the 
assumption that, whereas the unit buoyant 
force —Ay is sufficiently great to produce 
vertical acceleration, the corresponding varia- 
tion in the mass density of the fluid under- 
going acceleration is sufficiently small, in 
comparison with the density itself, to be 
neglected. 

The foregoing assumption as to constancy 
of the density makes it possible to proceed 
in close accordance with the elementary anal- 
ysis of flow in boundary layers, jets, and 
wakes. Three additional approximations are 
thus involved: the pressure intensity is assumed 
to be hydrostatically distributed throughout the 
field of motion; transverse forces are ignored in 
comparison with those in the vertical direc- 
tion; and turbulent mixing in the vertical 
direction is ignored in comparison with that in 
the horizontal. The fundamental equations 
of motion thus reduce to the equation of 
continuity and simplified equations of ver- 
tical acceleration and diffusion. The develop- 
ment of the relationships describing convec- 
tion above a line source will be described in 
some detail, but for purposes of brevity only 
the results will be indicated for the case of 
axial symmetry. 

Under the foregoing assumptions, the 
equation of vertical acceleration in two- 
dimensional flow reduces to 


Ou By ou IT 
Ox 


ou 


in which the intensity of vertical shear + is 
proportional to the mean product of the tur- 
bulent velocity components u’ and v’: 


T=—ouv 
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The equation of continuity is simply 


nt 0 () 


Finally, the equation of diffusion of the 
quantity Ay has the form 


x 0 (Ay) a 


ne 


g 2 (Ay) ow 
dy x dy 


(3) 


in which the rate of lateral diffusion w due 
to turbulence is proportional to the mean 
product of the eddy scale / and the lateral 
velocity of fluctuation: 


er À / 
7 2 (dr) 


w= i 

Since 4 = 0 at y= ©, v =0 at y =0, 
and t = 0 at y = Oo and y = ©, integration 
of Eq. (1), through use of Eq. (2), over a 
horizontal plane from the axis to infinity 
results in the following form of the momen- 
tum relationship: 


oo [ee] 


d 
fer dy =— | yay (4) 


O 


This states that the vertical gradient of the 
momentum flux is equal to the buoyancy of a 
horizontal stratum of unit thickness. As Ay=o 
ao y= CO and af = © at y =o and at 

= ©, a similar integration of Eq. (3) leads to 
the following expression for the constancy o 
the incremental weight flux past all successive 
planes: 


co 


ff wardy = (5) 


Oo 


Finally, after multiplication of each term of 
Eq. (1) by u, integration of the resulting ex- 
pression by parts, with the help of Eq. (1), 
yields the following form of the energy rela- 
tionship: 


oo oo [ee] 
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Because the quantity v® is negligible in com- 
parison with uw, this relationship states that 
the vertical gradient of the flux of kinetic 
energy is equal to the rate at which work is 
done by the buoyant force less the rate at 
which work is done by the turbulent shear. 
The last term at the same time represents 
the rate at which energy is lost to the mean 
motion through the generation of turbulence. 

Although no indication is given as to the 
form of the functional relationships for u, 
Ay, and t which will satisfy these three 
necessary conditions, the customary hypoth- 
esis of dynamic similarity of the mean (and 
turbulent) motion. at all elevations permits 
each of the relationships to be considered of 
the same form for every value of x. It is thus 
finally assumed that 


—_ =f (n) 


u 


A 2 
—<g9(n) ou. 2/2 = h(n) 


max Y max S 


Herein 7 = y/o, the quantity o representing 
some linear characteristic of the velocity pro- 
file u = u(y), such as the value of y at which 
the ratio u/umax has some arbitrary magni- 
tude. 

Inasmuch as the unknown functions f, g, 
and h are independent of elevation, any power 
or combination thereof must also be constant 
with x. It is hence convenient to formulate a 
series of pertinent integrals for substitution 
into the three primary relationships which the 
functions must satisfy : 


ne [ far = [Pan = [Pan 
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Introduction of these terms into Eqs. (4), 
(s), and (6) permits the latter to be rewritten as 


d 
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Simultaneous solution of these relationships 
will yield the following significant results: 
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Evidently, no matter what forms the distribu- 
tion curves may take, the convection zone 
will expand linearly with elevation; the maxi- 
mum velocity (and hence the velocity along 
any line of constant y/x) will be independent 
of elevation; and the maximum incremental 
weight density (and hence that along any 
line of constant y/x) will vary inversely with 
elevation. 

With this information it is possible to reach 
corresponding conclusions as to the variation of 
the volume flux Q, the momentum flux M, 
the kinetic-energy flux E, the unit buoyant 
force F, and the flux W of the incremental 
weight above a source of length L. Thus, per 
unit length of source, 


sols i à ae 2d mx 
T 2 u dy = x a] dy = x 
12 hose F r 
——2]|- mx — = —2 | Aydy~ x 
T. f& dy ~~ x L fay x 
O oO 
7 =2 | udydy = 20, 
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Herefrom it is seen that the flux of volume, 
the flux of momentum, and the flux of kine- 
tic energy must increase continuously with 
elevation, for the force producing the con- 
vective motion is the same at all successive 
levels. As the flux of the incremental weight is 
also the same at every level, the quantity W 
must indicate the output of the source. 

The corresponding analysis of the axially 
symmetric pattern of convection above a 
point source proceeds in a closely comparable 
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manner. The resulting integral equations of 
momentum, diffusion, and energy have the 
forms 


fer = — f Ay rd (7) 
d = 
7 [waved =o (8) 
d how i if ou 
Fr era — [uayra— fr Sr (9) 


Through the assumption of dynamic similarity 
of the mean flow at all elevations it is found 
that 
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The following expressions are then obtained: 
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Owing to the differences between these 
expressions and their counterparts for the two- 
dimensional case, basically different relation- 
ships are indicated for the flux of the several 
characteristic quantities past successive planes 
above a point source. Although the convection 
pattern again expands linearly with height 
above the source, the maximum velocity 
(and hence that for any constant value of 
r/x) now decreases continuously, and the in- 
cremental weight density decreases even more 
rapidly. On the other hand, the buoyant force 
on a stratum of unit thickness becomes steadily 
greater with increasing elevation, and Q, M, 
and E all increase at different rates. Once 
again, however, the flux of the incremental 
weight is the same at all elevations, so that W 
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must be equivalent to the output of the 
source. 


Experimental Procedure and Results 


In any experimental program to determine 
explicit relationships for either two-dimen- 
sional or axially symmetric convection from a 
boundary source, there are evidently four 
variables which can be controlled independ- 
ently — two coordinate positions, the mass 
density of the fluid, and the strength of the 
source. These four quantities govern the mag- 
nitudes of such dependent variables as veloc- 
ity, weight density, shear, turbulence, and 
their various combinations. While all of these 
factors would be of interest in a complete 
investigation of the phenomenon, attention 
is focused in the present study upon the major 
characteristics of the mean motion: the veloc- 
ity component # and the change Ay in weight 
density. The resulting functional relationships 
between the dependent and independent var- 
iables for the two types of source can be 
indicated schematically as follows: 


F W/L 
Pure We) 


Inasmuch as five variables involving three 
fundamental dimensions appear in each of 
these four functions, it is to be expected that 
the dimensionless form of each will contain 
only two terms. Thus, for two-dimensional 
conditions, the relationships may be expressed 
as 


Lan) 
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For conditions of axial symmetry, the counter- 
parts of these expressions are 


In each instance the grouping of terms is 
seen to be in accordance with the general 
relationships obtained analytically. 
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Experiments with heated air to determine 
the forms of these functions were conducted 
in closed rooms, particular care being required 
in the study of axially symmetric conditions to 
eliminate all sources of disturbance. In this 
instance the room was approximately circular 
in plan and had a diameter of 25 feet and an 
overall height of ır feet. The heat source was 
placed at the midpoint of a centrally located 
platform 8 feet in diameter and a few inches 
in height. For the two-dimensional study the 
flow was confined between two parallel 
walls 4 feet high, 8 feet long, and 4 feet apart, 
and the source was placed across the midsection 
of a low platform extending the length of the 
walls. The heat sources in both cases con- 
sisted of recessed gas burners yielding low, 
blue flames approaching as nearly as practi- 
cable the desired point and line concentrations. 

Measurements of the temperature distribu- 
tion were made by means of a copper-constan- 
tan thermocouple, with its cold junction 
placed well away from the heat source, in 
combination with a potentiometer reading to . 
0.002 millivolt. Velocity indications were 
obtained with a specially constructed vane 
anemometer 11/, inch in diameter, the jewel 
bearings of which permitted the indication of 
velocities as low as 0.2 foot per second. These 
instruments were mounted alternately on 
remotely controlled traversing mechanisms. 

The rates of heat output from both the line 
and the point source were evaluated from the 
measured distributions of velocity and tem- 
perature according to the following thermal 
counterparts of the equations for W/L and W: 


jl 
CR 2,0 fuATdy 
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(10) 
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for the closely approximate condition that 
Ayly = —AT/T, the corresponding value 
of W in either case was computed from the 
resulting conversion equation, 


Te te 


Supplementary tests performed during the 
axially symmetric study were directed toward 
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Fig. 2. Distribution functions for two-dimensional convection. 


determination of the level at which, under 
ordinary circumstances, the initial laminar 
flow could be expected to become unstable. 
Because of the low rates of heat input required 
to yield readily measurable conditions of 
transition, it was necessary to use as the limit- 
ing heat source a burning cigarette (from a 
package carefully calibrated for rate of heat 
release), the smoke serving as a convenient 
tracing agent. The critical Reynolds number 
determined in this manner was of the form 


Xe 


R tamis al 
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in which » is the kinematic viscosity. As can 
be determined herefrom, the critical height 
x. was small in every case for which detailed 
measurements were made. 

All velocity and temperature traverses 
were carried entirely across the convection 
zone, the plotted curves being essentially 
symmetrical about the maximum but often 
indicating a slight deviation of the flow from 
the vertical. In analyzing the results, the refer- 
ence axis was adjusted to agree with the 
axis of symmetry. Upon reduction of data 
to the pertinent dimensionless forms, com- 
posite plots were made of velocity and weight 
distributions for both the two-dimensional 


and the axially symmetric cases, to verify 
the assumed similarity of flow at different 
levels and with different strengths of source 
and to determine the explicit distribution func- 
tions. For ease in comparison as well as econ- 
omy of space, these plots are shown in half 
section: in Fig. 2 for two dimensions and in 
Fig. 3 for axial symmetry. 

At once apparent is the appreciable scatter of 
points which seems to typify experimental stu- 
dies of this nature — in part because of the 
difficulty of precise measurement and in part 
because of the sensitivity of such flow to slight 
disturbances. The scatter is not appreciably 
systematic, however, nor is it difficult to 
construct mean curves to indicate the distri- 
bution functions. The curves shown are nor- 
mal-probability functions which best fit the 
data and at the same time satisfy Eqs. (4) and 
(s), and Eqs. (7) and (8), respectively. At least 
as a first approximation, there seems little 
doubt that the plotted curves describe the 
basic free-convection phenomenon. To be 
remarked, however, is the fact that use of the 
Gaussian function — not to mention the sim- 
plified analysis itself — loses significance as 
y/x or r/x becomes large. 

If the probability curves plotted in Fig. 2 
are assumed to represent the true distributions 
of u and Ay, the corresponding expressions 
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Fig. 3. Distribution functions for axially symmetric convection. 
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can be used conveniently to evaluate other 
characteristics of the mean pattern of motion. 
Thus, since y = fudy, the stream function is 
expressible implicitly as 


Ÿ 
SST Le 2 
y= 1.80 ( ) fe »3)% (15) 
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in which values of the integral are obtainable 
from probability tables. The lateral velocity 
of entrainment just beyond the diffusion zone, 
otherwise available from the condition that 
v = — ÿy/ox, is instead computed from the 
expression 2v, = —d (Q/L)/dx as 


a "/s 
v,= + 0.28 ( 4 (16) 


Specific relationships for Q, M, E, and F per 
unit length of source are obtained by integra- 
tion as follows: 
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For conditions of axial symmetry these 
expressions take the following forms: 


was 1 
WIRT (—*) exp (—96 5] (21) 


sb 2/5 2 
Ay=-—-11.0 Bes exp (es 5) (22) 


— Wx5\ ls r? 
y=0.024( : ef: exp( 065) | (23) 


208 
400 


HUNTER ROUSE, C. S. YIH, AND H. W. HUMPHREYS 


400 


300 


250 


se Kane 
Cw/ Ly 9%)% 


200 


150 


100 


30 


100 


y 
WEY Gl 


Fig. 4. Convection pattern over a line source. 


M = 0.36 (—W)'! 9's xs (26) 
E=0.57 (—W) x (27) 
F = 0.49 (— W)' op": x's (28) 


So far as the mean flow is concerned, it 
remains only to indicate in generalized coor- 
dinates the pattern of stream lines and lines 
of constant incremental weight for the two 
cases. This, however, requires three rather 
than two dimensionless groups of variables 
in each instance — the additional variable 
which is needed being a characteristic length to 
serve as a reference coordinate. Such a length 
is found in the quantity (W/Lyg';)' for two 
dimensions, and in the corresponding quantity 
(W/yg')"s for axial symmetry. Upon intro- 
duction of these quantities, the dimensionless 
coordinate relationships for y and Ay be- 
come, for the line source, 


Y = | x n = Y | 
WILy °° L(W/Lyg')'!? (W/Lyg"ÿh 


Ay | a y | 
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and, for the point source, 
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Why "wg (Wet 


Ay _ \ | d X_ j LA | 
» PL(Wpe Te (Wiyg)" 


The corresponding plots of stream lines and 
lines of constant incremental weight are shown 
in half section in Figs. 4 and 5. 
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As can be seen at once from these undis- 
torted diagrams, the zones of pronounced ver- 
tical motion and change in weight density are 
concentrated about the area of symmetry and 
expand only slowly with elevation. Because 
of the fact that the measurements were restrict- 
ed to these zones and that the analysis itself 
is not exact, the outlying portions of the 
stream-line patterns are purely schematic. 
Aside from this qualitative indication of the 
pattern of inflow, the diagrams are considered 
to represent with good approximation the 
basic characteristics of convection at any scale. 
They may readily be interpolated or extra- 
polated, as all curves in any diagram are 
geometrically similar by virtue of the dynamic 
similarity of the convection mechanism at 
successive elevations. 
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Conclusions 


From the foregoing analytical and experi- 
mental studies of gravitational convection 
from both à line source and a point source, the 
following conclusions may be drawn: 


(1) The motivating force in gravitational 
convection is the positive or negative 
buoyancy of the fluid, and the primary 
characteristics of the convection pattern can 
hence be analyzed without regard to the 
means whereby the buoyant effect is pro- 
duced. 

(2) Through use of an approximate method 
of analysis similar to that for boundary 
layers, wakes, and jets, it is possible to ex- 
press as proportionalities the mean charac- 
teristics of gravitational convection for both 
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two-dimensional and axially symmetric 
turbulent motion. 

(3) Laboratory measurements of the mean 
convection patterns above boundary sources 
of heat have verified the approximate 
analysis and provided explicit relationships 
for the primary flow characteristics. 

(4) Although both the analysis and the ex- 
periments represent a considerable simpli- 
fication of the actual state of motion, the 
correlated results should permit evaluation 


1 Shortly before proof of the foregoing paper was 
received, there came to the authors’ attention a German 
wartime publication on the same subject, the analysis 
proceeding from assumptions as to the secondary rather 
than the primary characteristics of the motion: SCHMIDT, 


W., 1941, Turbulente Ausbreitung eines Stromes er- 
hitzter Luft, Z. angew. Math. Mech. 21, 265—278, 
351 — 363. 
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of the role played by the primary convec- 


tion process in more involved natural 
phenomena. 
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Abstract 


During a cruise across the Gulf Stream in October 1950 measurements of surface velocity 
were made both with the Loran-space-dead method and the electromagnétic method. 
A short account of the results is given with special reference to the velocity profile and 
the horizontal shear across the Gulf Stream. 


The transverse surface velocity profile 


In October, 1950, a short cruise was made 
across the Gulf Stream south of Cape Cod to 
observe the transverse surface velocity profile. 
Current velocities were observed by means of 
dynamic sections (Atlantic stations 4853 to 
4865), the Loran space dead reckoning method, 
and the electromagnetic method (von Arx, 
1950). A preliminary reconnaissance of the 
area revealed an anticyclonic bend which was 
then crossed repeatedly at nearly right angles 
to the flow. Some of the results, shown in 
Figure 1, disagree among themselves but 
suggest the narrow width and sharply peaked 
form of profile which seems possibly to be 
characteristic. 

Disagreement among observations stems 
partly from the average error of Loran, which 
at a range of several hundred kilometers from 


t Contribution No. 621 from the Woods Hole Oceano- 
graphic Institution. This work was supported by the 
Office of Naval Research under Contraet No. N6onr- 
27701 (NR-083-004) with the Woods Hole Oceano- 
graphic Institution. 


the transmitters amounts to about + 1 kilo- 
meter. This amount of uncertainty in the posi- 
tion of each hydrographic station and of the 
drift of the ship while on station in the strong- 
est part of the current, leads to a relative error 
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Fig. 1. Comparison of simultaneous measurements of 
surface velocity across the width of the Gulf Stream 
near the 69° W meridian, made by navigational, dynamic, 
and electromagnetic means. The area. under a curve 
drawn through the eastward phase of each set of observa- 
tions is 6.2, 6.0, and 5.4108 cm?/sec., respectively. 
The dynamic observations were made and computed 
by Mr. L. V. Worthington. All observations made in 
October 1950 from R. V. ATLANTIS. 
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Fig. 2. Results of simultaneous electromagnetic measurements of the velocity at the surface and 100 meter 

levels. Open symbols denote points at which the direction of flow was measured. Closed symbols and the 

trend of the lines joining them denote the current speed at intermediate points assaming that the direction re- 

mained the same or changed uniformly between fixes of direction. The abscissa was taken at right angles to 

the direction of flow at fixes 128 and 129. The ordinates shown are 23 per cent less than the average of navigational 
measurements of current speed. 


in the distance between stations 5 km apart 
perhaps as great as 20 or 30 percent. Similarly 
the errors of Loran and dead reckoning for 
navigational current measurements made On 
an hourly basis, may approach so cm/sec and 
become proportionally larger as the time 
interval is decreased. Measurements of sutface 
current by the electromagnetic method are 
reduced in magnitude when the current struc- 
ture extends downward through more than a 
negligible fraction of the total depth of water 
(STOMMEL 1948; Markus and STERN 1952). In 
that the current structure of the Gulf Stream 
may extend a kilometer or more below the 
surface, and the loss of signal strength is 
roughly proportional to the ratio of the depth 
of quiet water below the current to the total 
depth, decrements as great as 40 percent are 


sometimes encountered in the Gulf Stream 
off the New England coast. 

In the face of these uncertainties it has been 
necessary to counteract the limitations of one 
method with the advantages of another. In 
fine weather the evidence of navigation is a 
useful index of surface motion, but the method 
lacks resolving power. The electromagnetic 
method will RER direction of hie and 
detailed changes in one component of motion, 
but requires “calibration in magnitude. Both 
methods can be used underway. 

Figure 2 shows the uncorrected results of 
electromagnetic sections made simultaneously 
at the Saisie and 100 meter levels. Navigatio- 
nal measurements made on the same traverse 
indicate that the electromagnetic measure- 
ments of surface current speed averaged 23 
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Fig. 3. Comparison of the intensity of cyclonic shear (between the Gulf Stream front and the current maximum) 
with the curvature of the front. Curves A and B in the upper frame represent the smoothed results of electro- 
magnetic and navigational measurements of shear, respectively. The chart of the position of the Gulf Stream 
front in the lower frame was made from measurements of the greatest change of average temperature in the 
upper 200 meter layer. ‘‘Edgar’’ is a deep cyclonic bend that was observed to detach from the main current. 
Within this feature navigational measurements of shear were 1.7 times greater than electromagnetic measurements. 


ercent too small. Nevertheless a sharply 
peaked distribution of velocity across the cur- 
rent seems evident. The similar pattern of 
motion at the 100 meter level was displaced 
to the right of the surface profile along a 
seaward slope of approximately 1/80. The 
horizontal normal shear on either side of the 
current maximum was asymmetrical. 


Horizontal shear 


Data obtained during the multiple ship 
survey of the Gulf Stream of 1950 (FUGLISTER 


and WORTHINGTON, 1951) made possible a 
study of the change of the intensity of hori- 
zontal shear along the meandering current. 
The intensity of shear was first evaluated from a 
comparison of Loran and dead reckoning navi- 
gation, and then corrected for the local obliquity 
of the traverse to the direction of flow as 
indicated by electromagnetic measurements. 
In order that the curvature of flow might be 
neglected, the shear was evaluated and correct- 
ed for the obliquity of traverse to the direction 
of flow in increments of one or two hour’s run. 


214 


It was found that the average intensity of hori- 
zontal shear normal to the flow on the right 
of the current maximum was equal to or less 
than the Coriolis parameter, while the average 
values of cyclonic shear on the left were equal 
to or greater than this. Such asymmetry of 
shear has been suggested from theory by 
Haurwitz and PANOFSKY (1950). 

These data from the multiple ship survey 
also show that the intensity of cyclonic shear 
on the left of the current maximum varies 
systematically with the sense of curvature of 
flow. Greater than average intensities seemed 
to be associated with cyclonic curvatures, 
while less than average intensities seemed to 
occur where the sense of curvature was anti- 
cyclonic. Due to the unsuitable distribution of 
observations, it was not possible to test for a 
similar correspondence of curvature with the 
intensity of anticyclonic shear on the right of 
the current maximum, or for the persistance 
of a seaward counter current. 

In Figure 3 the upper frame contains the 
measurements of horizontal cyclonic shear, 
‘and the lower frame a schematic chart of the 
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curvatures of the Gulf Stream front during the 
period when these observations were made. 
It can be seen from a comparison of the curves 
in the upper and lower frames that the inten- 
sity of shear does not seem to depend on the 
length of the radius of curvature of flow. 
The points through which curve A was drawn 
represent electromagnetic measurements of 
shear, and curve B represents the smoothed 
result of the navigational measurements. There 
is a systematic increase in the proportion of 
disagreement between the electromagnetic and 
navigational estimates of shear in the vicinity 
of cyclonic curvatures which suggests that the 
current extends to greater depths in the vicinity 
of cyclonic curvatures. Such an interpretation 
is consistent with an increase in cyclonic shear 
since vertical stretching and horizontal con- 
vergence of the flow in cyclonic curvatures 
would result in increased cyclonic vorticity 
on the left of the current maximum. There is 
some evidence that the current maximum 
itself shifts toward the left as well as that 
surface current speeds are a little higher in 
cyclonic curvatures of flow. 
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Introduction 


The purpose of this paper is to present a 
survey of published data on chemical analysis 
of atmospheric precipitation; to discuss dif- 
ferent cycles proposed for inorganic com- 
pounds in atmospheric precipitation and, 
finally, to assess the importance of these com- 
pounds from different points of view. In- 
vestigations of rainwater with a view to de- 
termine atmospheric pollution in industrial 
areas will be mentioned incidentally but no 
results will be given in detail as they concern 
only small areas and abnormal conditions. The 
literature on this special subject is, furthermore, 
rather voluminous and should sell deserve a 
separate review. 


A. Historic review 


The historic review presented below is 
largely built on existing reviews [see e. g. 
Miter (213, 214), RUSSELL & RICHARD (254) and 
Way (298):] and will therefore not contribute 
<n new. The history of the research on 
chemical compounds in rain and snow is the 
history of the research on nitrogenous com- 
pounds in atmospheric precipitation. For a 


1 References are listed in a bibliography to be published 
at the end of the second part of this study. 


long time ammonia and nitrate nitrogen were 
the only compounds of interest while anal- 
ysis for compounds like sulfate, alkali and 
alkaline earth metals have been carried out 
practically only for the last three decades. 
Chloride takes an intermediate position. 

The first discovery of nitrate in rainwater 
was made by the German chemist MARGGRAF 
in the winter 1749—50 and in 1780—go it was 
rediscovered by a Swede, TORBERN BERGMAN. 
In the beginning of the 1800’s DE SAUSSURE 
in France showed the presence of ammonia 
in the atmosphere. Thus the two compounds 
of greatest interest were known to be present 
in the atmosphere at a time when the revolu- 
tion in agricultural chemistry took place. 

This revolution was mainly due to two other 
very important discoveries. In 1799 a Dutch, 
Jouan INGENHOUSS, published his studies abate 
the assimilation of CO, from air by plants 
and five years later DE SAUSSURE in France 
showed that inorganic constituents of plants 
were taken up by the plant roots from the 
soil. These two discoveries founded our present 
concepts of plant nutrition but it took, how- 
ever, about half a century before they were 
accepted in the scientific world. 

An important man in conveying these ideas 
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was no doubt Justus von Ltesic although his 
role to a great extent was that of a propa- 
gandist. The actual outline of modern agri- 
cultural chemistry was made by another 
German, CARL SPRENGEL. At about the same 
time also BoussINGAULT in France, although 
leaning towards the old humus theory of plant 
nutrition, made a valuable contribution through 
some classical field experiments. Among other 
things he observed the difference between 
legumes and non-legumes with respect to 
nitrogen; the legumes were able to use atmos- 
pheric nitrogen for their nutrition. Despite of 
many years work he failed to find the reason 
for this difference. 

BoussINGAULT was also the first to make 
analyses of rainwater in order to estimate the 
yearly contribution of combined nitrogen from 
rainwater. 

In 1840 Liesic published his well-known “Die 
organische Chemie in ihre Anwending auf 
Agrikultur und Physiologie”. In this work he 
strongly criticizes earlier concepts of plant 
nutrition and stated the modern concepts. In 
a later edition he changed his views on the 
uptake of nitrogen by plants. As he had found 
ammonia in plant sap he concluded that this 
compound, which was then known to be 
present in the atmosphere, was assimilated by 
the plants directly from the.atmosphere in the 
same way as carbon dioxide. This was, of 
course, a rather obvious conclusion at that 
time. It may be of interest to cite him in a 
letter to “Revue Scientifique et Industrielle” 
given in English translation in the Farmer’s 
Magazine (182). He states — “‘if the soil be suit- 
able, if it contains a sufficient quantity of 
alkalies, phosphates and sulfates, nothing will 
be wanting; the plants will derive their 
ammonia from the air as they do carbonic 
acid. We know well that they are endowed 
with the faculty of assimilating these two 
aliments; and I really cannot sce why we 
should search for their presence in the manures 
we use” .— “Ammonia added with the manure 
may be useful but is certainly not necessary.” 

His quoted data on the amount of combined 
nitrogen which was brought down to the soil 
yearly with precipitation was very high, 
27 kg/ha (1 ha = 10% m?), and it was mainly 
this figure which initiated the many long- 
termed research projects on combined nitrogen 
in rainwater which were carried out during 
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the later part of the 19. century. Among these 
the work carried out at Rothamsted is probably 
best known. (See Lawes & GILBERT (171), 
Lawes, GILBERT & WARINGTON (173, 174), 
MILER (213), Russert & RICHARD (254), Wa- 
RINGTON (296), Way (208, 299, 300). During 
the discussion which took place Liesic’s theory 
about ammonia in plant nutrition was turned 
down on the basis of data from field experi- 
ments at Rothamsted [Lawes & GILBERT (171)], 
and his estimate of combined nitrogen in 
rainwater was shown to be too high as the 
yearly data from Rothamsted showed. 

It is very probable that this controversy 
between LIEBIG and the Rothamsted scientists, 
combined with the increasing interest in the 
new concepts of plant nutrition, initiated many 
of the investigations which were carried out in 
Europe during this period. As will be shown 
below, similar investigations in U.S.A. were 
started much later. 

Most of the inorganic compounds present 
in rainwater were known before 1850, but 
naturally, they did not arouse the same interests 
as nitrogen. 


B. Nitrogen compounds in rainwater 
I. Introduction 


Most of the investigations carried out on 
rainwater concern nitrogen compounds, most 
frequently ammonia and nitrate, although a 
few cases also included nitrite and organically 
bound nitrogen, generally called albuminoid 
ammonia. Most of these investigations were 
carried out in the period 1870—1920. Fig. 1 
shows the frequency of such investigations 
during successive ro-year periods, the fre- 
quency being given as the total number of 
years of investigation within each period. It 
will be seen from this figure that a maximum 
occurs in 1881—1890. The distribution for 
Europe is, further, approximatly normal while 
for other countries it is asymmetrical, starting 
with a relatively great number of years of 
investigations in the 1881—1890 period. In the 
period 1911—1920 there is a minimum in the 
number of investigations in Europe, probably 
connected with the occurrence of the first 
world war. 

The contribution. from Rothamsted is seen 
to take a central part in the distribution. 
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2. Nitrogen as ammonia and nitrate 


a. Sampling and analytical errors 


Before discussing the results some remarks 
on the sampling, the analytical methods and 
their errors will be made. 

There are obviously two groups of errors 
which can be introduced; one group consists 
of errors due to the sampling technique and 
the storage of rainwater samples; the second 
group consists of errors in the analytical pro- 
cedure used. 

The sampling errors can be classified into 
inevitable errors and non-inevitable errors. 
Dust carried from the soil to the atmosphere 
(cyclic dust) and precipitated together with 
rain is an inevitable source of error as the dust 
particles may carry with them absorbed 
ammonia and nitrate or organic nitrogen com- 
pounds which are liable to give off their 
ammonia to the rainwater during the analytical 
procedure. The magnitude of this source may 
be considerable but is difficult to estimate. 
Existing data on yearly income of nitrogen 
compounds with the atmospheric precipita- 
tion may represent the total income but not 
the net amount, which from an agricultural 

oint of view is of greater interest. 

Cyclic dust seems to be the only inevitable 
error. Sampling errors which can be avoided 
include 


a. Bird droppings. 


b. Losses of ammonia from the sampling 
device due to aeration. 


c. Transformation of the nitrogenous com- 
pounds in the sampling device by reactions 
with the sample and the material of the 
device or by microbiological activities 
during storage. 


Bird droppings can produce errors in two 
ways, firstly by giving off nitrogen compounds 
as ammonia or as nitrate during the storage 
or the analytical procedure, and secondly by 
giving off substances which interfere with the 
analytical method used. In the second case it 
is, however, an error of the analytical method. 

Bird droppings certainly can increase the 
content of nitrogenous compounds in rain- 
water. JOHNSON (154) compared protected and 


unprotected gauges and found that the bird 


droppings in the unprotected gauge raised 
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Fig. 1. Number of years of investigation (station-years) 
on nitrogen compounds in precipitation for successive 
10-year periods. 


excluding Rothamsted 


the ammonia content considerably. The nitrate 
content was not materially affected. 

Harrison & WILLIAMS (133) found definitely 
higher nitrogen contents when using copper 
funnels for sampling rainwater than when 
using glass funnels. They explained this as due 
to a greater contamination by bird droppings 
in the copper funnels. The birds found the 
copper funnels suitable as a resting place while 
they probably got scared by the light reflec- 
tions in the glass funnels. 

The error produced by bird droppings will, 
of course, depend upon the population density 
of birds at the place concerned and further 
upon the habits of the different species present. 
At seaside places the errors caused by birds 
may be very great. In the second paper by 
MILLER (214) concerning investigations of rain- 
water from the west coast of Scotland several 
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samples were contaminated, apparently by 
bird droppings, to such a degree that the anal- 
yses and to be excluded. In inland places the 
errors due to bird droppings may not be great. 
A simple and effective device for preventing 
birds from using sampling gauges as resting 
places has recently been described. 

Losses of ammonia by aeration can be pre- 
vented by a suitable design of the collecting 
gauge. If the rainwater reacts with the material 
of the gauge the pH of the liquid will generally 


rise and thus enhance the ammonia losses. 


Reactions between the sampling gauge and 
the rainwater have been studied to some extent. 
CROWTHER & STUART (82 cf. 81) found that 
copper containers were likely to neutralize 
rainwater, especially acid rainwater and in the 
above mentioned Swedish investigation it has 
been shown that galvanized iron in sampling 
gauges is liable to cause large errors due to 
reactions between Zn and water. As Zn(OH): 
is formed, the pH will increase to such values 
that losses of ammonia can occur; further the 
hydrogen formed in this reaction is very active 
and reduces nitrate almost quantitatively to 
ammonia even when present in such low con- 
centration as in rainwater. 


Microbiological activities may occur during 
storage of rainwater. HANSEN (130) found that 
toluol was effective in preventing microbio- 
logical activity during storage. In many cases 
lead-lined gauges have been used. The small 
amounts of lead which go into solution are 
then sufficient to prevent microbiological 
activities during storage. This type of con- 
tainer was used in Rothamsted. 


The analytical errors are due entirely to 
unsuitable methods for the analysis of ammonia 
and nitrate in rainwater. The most common 
methods for these compounds are colorimetric, 
using Nessler’s reagent for the determination 
of ammonia, and phenoldisulphonic acid for 
the nitrate. The Nessler reagent is rather 
sensitive to the organic compounds always 
present in rainwater. If the determination is 
made directly in the rainwater, by adding the 
reagent, large errors can be produced, the 
magnitude and sign of which will depend on 
the special type of reagent which can be made 
up in a number of ways. ROMELL (251) has 


! EGNER, ERIKSSON and EMANUELSSON, 1949, Ann. 
Agr. Coll. Sweden 16, 593—602. 
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recently shown that water contaminated by 
pollen ‘interferes with the Nessler reagent in 
the direct determination to a very high degree 
probably due to formaldehyde given off by 
the pollen grains. Many investigations have 
been carried out using this straight forward 
procedure and their results may be misleading. 
It is, further, probable that the errors are 
systematic and positive so that the given data 
are too high. | 

An indirect method using the Nessler rea- 
gent was used at Rothamsted. The ammonia 
in the rainwater was distilled off and the 
ammonia determined in the distilled portion 
by the reagent. This method, which has also 
been used in other places, gives quite reliable 
values. The ammonia in the distilled portion 
can also be determined by titration provided 
special care in taken during the distillation. 
This method has been used by Hansen (130, 
131) in Denmark and must be regarded as 
highly reliable. 

The direct determination of nitrate by the 
phenol-disulphonic acid has been widely used 
but is not satisfactory as losses of nitrate may 
occur during the procedure (130). In older 
investigations of this type the nitrate figures 
may therefore be too low. In Rothamsted 
WARINGTON (296) introduced a new method of 
determining nitrate by reducing it to ammonia 
by a copper-zinc couple, followed by distilla- 
tion and estimation of the ammonia by the 
Nessler reagent. Later Devarda’s alloy was 
used for the reduction. This process seems 
to be quite satisfactory and has been used also 
in other places. The technique of using ion 
exchange resins in the sampling of atmospheric 
precipitation was described in the paper by 
Egnér and coll. cited in footnote on this page. 

It is, of course, difficult to divide published 
rainwater nitrogen data into a reliable and an 
unreliable group based on the analytical 
methods used. In many cases no statement 
concerning the method is given, in other 
cases the data seem quite normal despite less 
satisfactory methods. There are so many 
factors involved (including the sampling 
factors) that rejection of data on the basis of 
unsatisfactory analytical methods cannot be 
done. The best way is to point out the most 
reliable investigations, basing the discussion 
mainly on these, leaving the rest as “ probable” 
or in single cases as “quite improbable”. 
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in precipitation plotted against the yearly precipitation. 
Filled circles: European places. 
Open » : non-European places. 


b. Results 


A survey of data up to 1900 has been given 
by Mitıer (213). The investigations published 
since then do not change the general impression 
although they are quite numerous. For the sake 
of completeness existing data have been com- 
piled; the averages for different periods and 
places being summarized in tables 1—3. 

The results in tables 1—3 have been arranged 
diagrammatically in figs. 2 and 3 to show the 
distribution of yearly amounts as a function 
of the rainfall. In these figures a division 
between European and non-European places 
has been made. 

The total nitrogen as ammonia and nitrate 
shows a very wide variation (fig. 2). There is 
a fairly marked minimum level of around 1.9 
kg N/ha and year and it may be concluded 
that this is probably the smallest amount 
which is brought down by precipitation (except 
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in deserts). The two exceptionally low values 
marked a) and b) in fig. 2 belong to coastal 
stations at Vifelstadir, Iceland, and Butt of 
Lewis in the Hebrides. It appears to be 
possible to draw another boundary for the 
distribution (the sloping line) which can be 
interpreted as the maximum nitrogen amount 
for a given yearly rainfall outside Europe. It 
is seen from fig. 2 that this maximum amount 
is greater in Europe than in other places. 
The places which lie above the sloping line 
are with one exception from Europe, namely 
from Russia, Germany (five from 1860—70 
and one from 1928—35) and Denmark. The 
exceptionally high value c) comes from Buenos 
Aires. In the Danish investigation (HANSEN 
129) the analyses were very accurate so there 
can be no analytical error involved, further 
there cannot be any great immediate cont- 
amination from industry. 

The nitrate-nitrogen distribution is shown in 
fig. 3 and it will be seen that the minimum 
amount is not so distinct as in the case of 
total nitrogen. The sloping line of the maxi- 
mum amount is, however, much more distinct 
than in fig. 2. On the other hand there is no 
marked difference between European and non- 
European places although there is a slight 
indication. The exceptionally high value a) 
is from an old German investigation (PROSKAU 
1864). 

It should be mentioned that the division of 
the data into a European and a non-European 
group has significance only from the point of 
density of population and industrial activities. 
Another possible division is tropical and non- 
tropical places; this is, however, of but little 
value as the investigations carried out in tropical 
places are rather few and incomplete. This can 
be seen in fig. 4, where the amount of total 
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Fig. 3. The same as fig. 2 for nitrate nitrogen. 
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Table 1. Ammonia and nitrate nitrogen in precipitation in Europe 
a POS De Ae BETRETEN EEE RAP PORN Des Bee A, RES en a ee en — 


Precip. mg N/l kg N/ha 
Ref. Place Period in as as as as as as 
mm | HN | NO, | tot. N| HN | NO, | tot. N 
| 
Austria 
103 Reldkirche rer er. 1876 — 2.79| 1.76| 4.55 = — — 
194 Scie Nichaelwrr eek 1885— 1887 TFILA ET 67.1, 1.84. 073.03 7.44| 20.47 
Belgium y : 
2251230 GÉMDIOU er re 1889—1893 693 1.14 935 1.49 7.92| 2.39], 7oSE 
146 Bulgare EA. — — — — | 17.86 .96 12.82 
Denmark 
291 Copenhagen sente. 1880/81—1884/85 S7 TAN T-OE AS) -2:36\\|\) 02.00) | 22.621070202 
131 ASKOV TEE Re eee 1922/23—1926/27 740 AE BS 1.06 5.25, 2,0302 7588 
131 BANDE RE EE Eee 1925/26 600 | 1.07 EST) 7,04 6.44] 3.41] 9,85 
131 Spangsbjeror nen a: 1925/26 730 .88 .60| 1.48 6.42 | 4.39] 10.81 
Ta Hlornumere. ARMES. 1925/26 495 74: .$4 1.28 3.66 | 2.65] 6.31 
Finland 
112 Stork TO ee — == — == 1.8 = = 
France : 
176 IMethro ya ters dees 1877 760 AT — _- QUE — — 
179—181 Montsourise rer. 2er 1876— 1900 545 ARTS .66 2.79 | 11.60 3.60] 15.20 
Germany 
213 Liebwerd, Bohemia....... 1877/78 620| 1.3 .6I 1.91 8.05 | 3.78| 11.83 
213 Pecek ee Ato. he. 1883/84, 1885/86 490 1.26 .$O 1.76 6.17 2.45| 8.62 
53 Ida-Marienhütter 20, . ve 186$—1870 580 — — 1.88 — — NITRO 
314 CUI HCH Bee 1864/65—1865/66 369 .48 .I6 .64 BT SSO - 234 
314 InstenDUr Sep ee 1864/65— 1865/66 642 .6S .39 1.04 41.15 2.49| 6.64 
314 DAME seen eee ton oe 1865 427 1:42 .30 1.72 6.07 1.28] 7.35 
314 Recehwaldesn atts arene 1864/65—1865/66 568] 2.08 62 | 2.70") 17:82 779.21 184 
314 Proskauer s ARE 1864/65 445 3.21 7073 4.94 | 14.29 7.70| 21.99 
255 \Wehenstephanıns. "7 1928—1935 750 — _ 1.97 — — | 14.80 
Great Britain 
254 Rottiamstedke MER wae oor 1888—1916 730 AI .20 .6I 2.98 1.48] 4.46 
81 Garforth, Leeds en. 1906—1909 685] 1.04 5321110 2236 7-09 | 2.16) 9.25 
214 andale Ard courses ei 1907 2,250 14 .06 .20 372 1.41] 4.53 
214 Barrahead, Berneray...... I910—1913 895 .IS .I4 .29 130.132 2.54 
214 Shilay, Mon. Islands...... 1910—1913 1,230 T2 .0$ 17 1.41 .66| 2.07 
214 Butt of Lewis, Stornoway 1908—1913 1,060 .04 .03 .07 .40 .34 .74 
Holland 
144 Groningetinme: MN 1910—1912 700 172 .24 .96 5.05 1.68| 6.73 
177 Hilversum (aver. of 6 places) 1932—1937 -- — — = 3.8 7 4.5 
Iceland 
214 Nitelstaditze ze ee ge I9II—I9I2 949 .09 .03 Re) 86 .28] 1.16 
Italy 
24 FIOTENS ER in RER ere er 1869—1875 970 1.00 457 1.57 9.70 5.55] 15.25 
24 Wall omibnosa ats ee. 1872—1875 1,520 .62 25 .87 9.40 | 3.80] 13.20 
232,233 SCADdLCC ee Ne eee ke 1888— 1890 740 OI AT .88 4.55 2.00] 6.55 
2 Catania (outside town) ... 1889/89 466 38 .I6 .49 1.52 TSI 0 
Norway 
252 IHRE ROE REN RS 1864/65 670 37 = sp air = = 
$I HTONADIEMEE Free 1928/29 889 eng .I0 .26 I.17 Ar) DE 207 À 
51 Voll(7kmeast of Trondhjem) 1928/29 758 Id .0$ .19 1.06 .38l 1.44 
Russia 
239 GOUCSS LR dence acc ee ER 1903 275 .49 .I8 .67 1:35 .§0] 1.85 
292 INOVOStl Enr ae ee m 1905 ee 
164 ER trate sauna Cr 1945/46—1946/47| = 400 6 a 7 235 4 2.0 
30I, 302 PISTES Winks cay een RE 1900—1903 445 85 .06 OI 3.88 27 US 
Sweden : 
106 HART tees taste 1909 827 45 .18 .63 272 TA CINE SES 
Switzerland 2 
1202712 1) er ae 1870/71 890 — .46 — — 4.05 — 
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Table 2. Ammonia and nitrate nitrogen in precipitation in Asia, Australia, New Zealand and Africa 
PR OP ER SE Re 


Precip. mg N/l kg N/ha 
Ref. Pema Ge Period in 
ae as as as as as as 
Asia 
East Indies 
200 Pasoeroean,Hjavatres...a. 1891/92 920 .I4 09 23 127 79| 2.06 
250 Sumatra elite PE. 192$—1940 2,045 — — 15 — — | 30.4 
80 AE ES RME PTS RE 1941 — == = = = el = 
India 
243 Allahabad seen dence de 1932 — .50 88| 1.38 _ — — 
213 (CERN ES RER En 1891 1,170 17 12 29 2.001 1.340 3.34 
FAT TS CAWDDORS RS ee 1904/05 1,260 22 06 28 213 .85| 3.58 
1370775 Pébr Dunes En. Li 1904/05 2,220 — — 18 — _— 4.02 
3 NASSEN ate 1888— 1893 995 = — 19 — — 1.91 
86 Sylheesafsın ah hand, 1931/32 3,950 ie} .09 22 5.06| 3.74| 8.80 
17, 60 Hanon Tonkın....s..s4s: 1902—1909 1,650 A PROVEN Pee TG OS SUIROT.S 
18 Geulon LANCE ARR 1898/09 2,090 .20 .07 27 | 4.09| 1.431 5.52 
162 » Reredeniva +. 1941/42—1942/43| 2,850 27 18 45 7.82 5.22| 13.04 
Japan 
157 TRS) Ce are ac 1883/84 1,337 pe) 07 20 1.67 97| 2.64 
216 NACRE EU SEE PE 1938 — .58 
216 EOD ema sed: 1938 — 28 — — 
216 Famamatın) «tba cso oe 1938 — 19 — — — — — 
Australia 
56 Brisbane, Queensland. .... 1908 1,150 22) 19 .4I 2.48, 2.7141 4.62 
s6 Roma, EN CORRE 1908 660 .42 27 -63 2.77. 1737 E40: 
New Zealand 
135 Lincoln, Canterbury...... 1884—1887 755 .08 17 25 0. :127| 81.83 
Africa 
156 Bloemfontein, South Africa] 1910/11—1911/12 550 .87 HO. ai 4,70 1m 1267 6:12 
156 Cedera » » 1912 680 .78 AL .92 5.27. .97| 6.24 
156 , Durban » » IQII—IY9I2 910 .47 8105 .62 4.23 1.39) 5.62 
156 Grahamstown, » » 1911/12 590 .18 .I4 32 1.05 82| 1.87 
156 Kokstad » » 1912 670 .19 zu .30 12$ 751. 2.00 
149 Pretoria » » 1904/05 620 1.19 .20 120 7.40 1.21| 8.61 
44 U LERIe UN nr ee ne 1895 (1,500) | (.43)| (.40)} (.83)| 6.45] 6.00| 12.45 
IOI De ot PET os 1938 1,002 3.38 — — 3.80 — — 
222 REUNION IE IRA 1886— 1887 (1,000) — .69 — = (6.9) — 


nitrogen and nitrate nitrogen is plotted against 
the latitude of the place in question. It will 
be seen that the tropical zone is so poorly 
represented that it is difficult to make a fair 
comparison between this zone and the tem- 
perate ones. Fig. 4 is interesting, however, as 
it shows a very high frequence of large total 
nitrogen values in the northern hemisphere. 
The nitrate nitrogen seems to be fairly evenly 
distributed throughout all the zones. 
Returning to figs. 2 and 3 the sloping border 
lines suggest that the variation in maximum 
nitrogen amount tends to be proportional to 
the amount of precipitation, or, expressed in 


a different way, that the maximum concentra- 
tion of nitrogen is independent of the annual 
amount of rainfall. The significant difference 
between European and non-European places 
is that the maximum nitrogen concentration 
in European rainfall is able to reach somewhat 
higher values than in non-European. For 
nitrate nitrogen there is no such pronounced 
difference (fig. 3). 

As the nitrate nitrogen does not show any 
pronounced difference in the distribution 
between European and non-European places 
the difference in fig. 2 must be due to the 
ammonia nitrogen; the ammonia nitrogen 
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Table 3. Ammonia and nitrate nitrogen in precipitation in America 
D RL ne ee ee Se eee 


Ref. ; Place Period 


Precip. mg N/l kg N/ha 
in as as as as as as 
mm | HN | NO, | tot. N| HN | NO, | tot. N 


ge" > 


North America 


271—275 Ottawa Ganadas zn 1908— 1924 860 .SI .28 .80 4.43 2.39] 6.82 
310 Alfred, N.Y. U.S.A....| 1923/24—1924/25 820] 1.67 .03 1701970 28| 13.98 
310 Brockport» » 1923/24—1925/26 800 .36 cok EES 2.91 2010807 

75 Geneva » » 1919—1928 900 .99 12 TE 8.91 1.07 32 
Tee, BAS Ithaca » » 1915/16—1925/26 750: 1.07 AD) Lebo 8.00 91| 8.91 
110 Goodwell, Oklahoma 
WES As Re RS 1930 390 .28 12 .40 eae .48] 1.59 
115 Kentucky, six places U.S.A. 1922/23 1,090 1.19 TA 1.932 13:00, 8:021027.02 
104 Manhattan, Kansas U.S.A. | 1886/87—1889/90 740 39 .16 aS 2.87| 1.16] 4.03 
TS DIS, 

LE LOO: 

TOE 5234: Mont Vernon, Iowa U.S.A.| 1912/13—1931/32 740 35 LS .$0 2.63 T:08 | Pagar 

246, 269, 

(8 years) 

290, 305, 

306, 313, 

102 Salt Lake City, Utah U.S.A. 1891—1893 5.65 .40| 6.05 

289 INHSSOULL AUS Are ee 1894—189$ 1,120 .24 .07 37 2.63 821 WAS 
West Indies 

133 Barbados. ra 1886— 1889 1,620 07 .20 27 FAIR 4.36 

61 dad ceteris Mee 190I— 1912 1,480 .16 .OI PL 2.37 Toll 2.47, 
South ‘America 

134 Georgetown, Br. Guiyana. 1890—1909 2,520 .0$ .08 LEA 113 2.06| 3.19 

196, 222 Caracas. Venezuela... 1889/90, 1895 (1,000) 1.55 ss 2.10. (ESS) 1) (SH REA) 
268 Montevideo, Uruguay .... 1912 1,504 24 27 SE 3:68 F4. 03 TT 

58 Buenos Aires, Argentine .. 1906 750| 2.96 .18 Sikh.) 222 T.33 253 


content varies more in Europe than in other 
places. 

The proportions of ammonia to nitrate 
nitrogen is shown in fig. $, arranged in the 
same way as in fig. 4. It has been concluded by 
MILLER (213), that trophical rains contain a 
higher percentage nitrate nitrogen than non- 
tropical. As seen in fig. 5 this seems to be true 
when comparing zones II and III and is prob- 
ably due to the earlier mentioned greater 
concentration of ammonia in European rain- 
fall. Comparing zones III and IV it seems, 
however, that there is no reason to believe 
the proportions to be really different. As 
a matter of fact there are only two places in 
the tropical zone that show a low ammonia 
percentage. 

In this connection it should be mentioned 
that some investigations show an abnormally 
high ammonia percentage, especielly those 
carried out in the state of New York, U.S.A. 
There is probably some error involved, either 
in the sampling or the analysis which have 


led to a reduction of nitrate to ammonia, 
or in the analysis, in which contamination 
present has affected the ammonia figures. The 
second explanation is more probable. It has 
been pointed out earlier by ROMELL (251), that 
the ammonia figures from these places are 
abnormally high and he suggested interferences 
in the direct determination of ammonia by 
Nessler’s reagent, caused by substances given 
off by pollen grains in the rainwater samples. 

It is difficult and probably not justified to 
work out an average of nitrogen brought to 
the soil by rainwater in different parts of the 
world as the variation between different places 
probably is more significant than an average. 


c. Monthly variation 


The amount of combined nitrogen brought 
to the soil with the rain each month will vary 
not only with the rainfall but also with the 
general climatic conditions and the habits and 
occupation of the population. For some of the 
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Fig. 4. The average yearly amounts of total nitrogen and 
nitrate nitrogen for different places plotted against their 


latitude. 
Filled circles: total nitrogen. 
Qpen Fa nitrate *) 


places recorded in tables 1—3 monthly averages 
are given and a few representative investiga- 
tions of these are shown in table 4. 

In Gembloux, a place where atmospheric 
pollution is fairly high, the monthly amount 
of nitrogen is fairly well correlated with the 
monthly amount of rainfall. The season has, 
however, also a direct influence on the amounts. 
This can already be seen from the figures in 
the table and is clearly shown by statistical 
analysis. Assuming that the relationship be- 
tween amount and precipitation is linear and 
the seasonal influence is a simple periodic one, 
the relation arrived at reads 


N = 0.274 cos 30 (T + 0.3) + 0.0143 p — 
BG FOG, eee Meme le. CRE. (A) 


where N is the amount of ammonia nitrogen 
in kg/ha/month, T is the number of the 
month starting with January equal to zero and 

is the monthly rainfall in mm. The argument 
of the cosine function is expressed in degrees. 
Thus, there is a higher ammonia content in 
the precipitation in winter than in summer 
which may be due to an increased coal com- 
bustion during the winter. 


28 


In Rothamsted there is a maximum in the 
amount of ammonia nitrogen in August and 
a minimum in February. The influence of the 
rainfall is obscured by the fact that it also 
varies periodically showing a maximum in 
October. The ammonia concentration during 
different months (weighed averages, not given 
in the table) shows, however, two maxima, 
one in the spring and one in the late summer 
as a generally somewhat higher level during 
the warm season than during the cold. 

A significant feature for Ottawa is that the 
rainfall varies very little from month to month. 
The variation in the amounts of ammonia 
nitrogen can therefore be ascribed mainly to 
seasonal influences. An analysis of the data 
gives a relation which reads 


N = 0.150 sin 30 (T— 4.1) + 0.099 sin 60 
RR EEE 


and is graphically illustrated in fig. 6. The 
seasonal variation has thus been resolved into 
two components, the first of which follows 
the expected yearly temperature of the air at 
such a place. The second one has two periods 
within the year with maxima in spring and 
autumn just as indicated by the figures on the 


oO 
ammonia content in Rothamsted. Both of the 
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Fig. 5. The percentage ammonia nitrogen of total for 
different places, plotted against their latitude. 
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Table 4. Monthly distribution of amounts of nitrogen compounds in precipitation in a few 
representative places 


Ref. Place and item Jan. | Febr. [March] April 
235 | Gembloux 
H;N—N, kg/ha...| .54 05 .79 .57 
NO;—N, kg/ha...| .19 .03 LS .18 
rainfalls, goer «7 44 18 $I 46 
254 | Rothamsted 
H,;N—N, kg/ha...| .21 .19 23 PR 
NO;—N, kg/ha...| .10 | .10 | .ır | .ıı 
rainfall mmer an. 55 $1 61 43 
273, | Ottawa 
H,;N—N, kg/ha...| .21 ATS 22 .43 
NO;—N, kg/ha...| .14 | .17 | .17 | .26 
rainfall mme... 81 64 72 65 
133 Barbados 
H3N—N, kg/ha...| .06 .0$ .0$ .05 
NO;—N, kg/ha... 12 .10 .10 10 
Tinta) eer eye 76 52 42 30 
133 | British Guiyana 
H3N—N, kg/ha...| .10 ‚16 9 1-19 
NO;—N, kg/ha...| .33 PASS || || re 
raintalls malen... 3301| 22S97 1 er S4an 205 
162 | Peredeniya : 
H,;N—N, kg/ha...| .57 AT .75 .26 
NO;—N, kg/ha...| .16 34. 1-09) | 229) 
FAP 64 45 146 | 279 


May | June | July | Aug. | Sept. | Oct. | Nov. Dec. 
553 .81 .8I .97 .42 .89 .76 .6I 
.32 .26 .29 43 217 23 MES 211 
59 93 | 104 85 45 62 52 42 
.26 27 .28 .31 .24 .26 .25 DS 
13 .13 =e 325 .12 .I8 .I2 .13 
SI 64 61 72 48 82 67 78 
.43 .40 .33 33 SE .6I .35 25 
.19 22 .28 SE 07 .20 .I4 #73 
74 72 70 74 68 80 69 71 
.06 08 17 .23 “12 .10 IT .06 
21 .18 .30 .40 21 2 .50 3 
93 114 877.1 2857.|7203 I4I 208 117 
.20 TE Sa: .19 .OI .07 ‘17 .09 
27 34 32 ES .08 .06 28 ty .19 
312 328 278 173 64 SI 142 380 
.46 .36 LME 371 v7 SSAVE23 .70 
.88 .58 LS .30 43 533 .54 CES 
603 296 112 183 203 349 359 226 


rainfall, mm 


sine functions in relation (B) are statistically 
significant. It would be of interest to see if 
this distribution of ammonia nitrogen might 
be related to the occurrence of forest fires. 

In the tropical places Barbados, British 
Guinea and Peredeniya the variation in rainfall 
is very great and obscures the direct influence 
of the season. A correlation between the 
amounts and the rainfall is indicated by the 
data. 

The nitrate nitrogen does, in general, show 
less influence of the season than ammonia but 
is fairly well correlated with the rainfall. From 
the Ottawa figures it can be seen, however, 
that there is a certain seasonal influence. A 
relation of type (B) for nitrate nitrogen reads 


N = 0.070 sin 30 (T — 2.6) + 0.022 sin 60 
(T + 0.3) + 0.199 


both sine functions being statistically signifi- 
cant although the contribution from the second 
one is rather small. The maximum due to the 


first component occurs between June and July; 
it differs from the ammonia relation in that 
respect by one and a half month. The maxima 
of the second component occurs in February 
and August. Its significance from meteoro- 
logical point of view is doubtful. 


d. Yearly variation and trends 


The variation of the yearly amounts of nit- 
rogen compounds at a certain place is in general 
fairly large. This may be due to a number 
of factors such as yearly rainfall, prevailing 
winds and the habits and occupation of the 
population dwelling in or around the place 
in question. These factors vary periodically 
but some may have a very long period and 
give rise to what is called yearly trends. 
A closer study of these trends necessitates 
long-termed series of observation, especially 
when the yearly variation is large. 

At two places, namely at Rothamsted and 
in Ottawa, the yearly amounts of ammonia 
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and nitrate nitrogen in precipitation have been 
estimated for relatively long periods. In 
Ottawa (from 1908 to 1924) the yearlyamounts 
have a tendency to increase with time. This 
increase is especially noticeable for the last 
four years. The average amount of ammonia 
nitrogen for the first ten years is 3.88 kg N/ha 
while the average for the whole period 
(17 years) is a high as 4.40 kg N/ha. The 
increase is probably due to increasing human 
activities. For nitrate nitrogen no trend at all 
is detectable. 


At Rothamsted (from 1888 to 1916) the 
amounts of ammonia does not show any 
clear trend but the nitrate nitrogen, as was 
also pointed out by Russet and RICHARD (254) 
increases slowley but steadily with time. 
The reason for this increase may be difficult 
to find. 


In the period 1908 to 1916 the two series 
overlap. There is, however, no obvious correla- 
tion neither between the amounts of ammonia 
nor between the amounts of nitrate at the two 
places. This is also interesting as it shows 
that the factors reponsible for the yearly 
variation may be fairly local in extension. 


3. Albuminoid nitrogen 


This name has been assigned to the ammonia 
released from organic matter in the rainwater 
when treating it with oxidizing agents, for 
instance KMnO,, or by Kjeldahl digestion. 
There are a few data available for albuminoid 
nitrogen in rainwater. In Mont Vernon, 
Iowa, (ref. 15, 114, 115, 160, 161, 234, 246, 
269, 290, 305, 306, and 313) it is estimated 
to be 2.7 kg/ha and year as an average for 
9 years and at an average rainfall of 755 mm. 
In Ottawa, Canada, (271—275) the yearly 
amount was 0.98 kg/ha as an average for 17 
years and a rainfall of 860 mm. At Dehra Dun 
in India (86), the amount for 1931—32 was 
6.05 kg/ha at 2,950 mm rainfall. In Rothamsted 
Miter (213) has estimated the amount to 1.5 
kg/ha/year, the estimate being based on 69 
analyses. In New Zealand Gray (125) found 
0,5 kg/ha. 

From these scattered data it is seen that the 
contribution may be considerable in some 
places. Some of this nitrogen is apparently 
derived from cyclic dust and, thus, does not 
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Fig. 6. The monthly amounts of ammonia nitrogen 
of Ottawa precipitation. 

Dashed lines: the two components of the full drawn 
line mentioned in the text. 


contribute anything to the nitrogen economy 


of the soil. 


4. Nitrite nitrogen 


The nitrite nitrogen, when reported, is in 
general very low. In some investigations it is 
automatically included in the nitrate fraction 
as the determination of nitrate is made by 
reducing it to ammonia by a metal alloy. 
Such a reduction obviously also includes the 
nitrite. 

FAILER and WILLARD (104) tested daily samples 
of rain in Kansas qualitatively for nitrite. Of 
171 samples 75 per cent gave a clear positive 
reaction, 8 per cent no reaction and the re- 
maining 17 per cent a doubtful reaction for 
nitrite. The frequency of nitrite present was 
higher during summer than during winter. 

The investigations carried out in Mont 
Vernon, Iowa, show variations in the daily 
samples from a trace to 0.051 mg nitrite 
nitrogen/l or roughly about 10 per cent of 
the nitrate nitrogen. The same was also esti- 
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mated at Sylhet, India (86) where the total 
nitrite nitrogen for a year amounted to 0.46 
kg/ha, most of it being delivered during 
April—June. Analysis of rainfalls near Bombay 
over a half year period gave nitrite nitrogen 
contents ranging from traces to 0.002 mg N/L. 
In Japan (216) the average contents for three 
places were 0.0027, 0.0040 and 0.0067 mg N/l 
as nitrite. An old French investigation (69) show 
comparatively high values, but from the more 
recent ones quoted here it seems likely that the 
amounts are relatively low. Some investiga- 
tions carried out in Australia will be discussed 
later. 

Most of the data quoted here were obtained 
on samples from single rainfalls, analysed very 
soon after the precipitation. Hence the figures 
should be quite reliable. 


5. Nitrogen compounds in fog, dew, frost and 
condensed water from the air 


The compositions of these forms of atmos- 
pheric precipitation have been investigated 
to some extent. 

a. Fog. In 1853 Bousstncautt (50) found 
more ammonia and nitrate in fog than in rain. 
The amount of ammonia ranged from 3—9 
mg N/l. This was later confirmed by Levy (180) 
who, however, found a low nitrate content 
in fog. In Belgium PETERMAN and GRAFTIAU 
(236) made three analyses which gave 2.60, 
5.64 and 5.00 mg total-N/l and in Bologna 
in Italy Casoxt (62) found values ranging from 
14—58 mg ammonia nitrogen/l. From Russia 
it is reported (244) that fog, dew and frost 
contained two and a half times as much am- 
monia as did rain and two and a half times 
less nitrate. 

b. Dew. Bousstncautt (50) found as much 
ammonia in dew as in fog but much less 
nitrate. In Cawnpore in India (137, 175) the 
analysis of dew gave 0.85—2.65 mg ammonia 
nitrogen and. 0; $n. 12 mg nitrate nitrogen/] 
of dew. For the six month period in 1904—05 
the amount of ammonia and nitrate nitrogen 
delivered to the soil by dew amounted to 0.06 
kg/ha of each. Also Levy (180) reported that 
dew was rich in ammonia. 

Recently QUITMANN and CAUER (241) com- 
pared the composition of condensation. water 
near the ground with that of rain on the same 
occasion. They got (two cases) 
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H3N mg/l | NO, mg/l 


Condensation water .| 10.0, 10.0 0.10. 0.05 
RAIDE Ans en 5 160 0.35, 0.10 
which shows a considerable difference in 


composition between these two forms of 
precipitation. - 

c. Frost. The Belgium (236) investigation 
gives an average of 7.52 mg total nitrogen 
per litre for five analyses; in Bologna (62) 
the range of concentration of ammonia in 
frost was 22.4 to 44.7 mg N/l. Some recent 
analyses from Germany of MROSE (218) show 
lower values with a variation ranging from 
o—4.8 mg ammonia nitrogen and 0—3.7 mg 
nitrate notrogen/l. 

From all these restilts it can be concluded 
that those forms of precipitate which are 
derived from near the soil surface as a rule 
are much richer in ammonia than ordinary 
rainwater. 

In this connection some observations on the 
rate of rainfall and size of raindrops on the 
ammonia content of rain should be mentioned. 
As early as 1870 Rosine (252) in Norway 
observed that a very fine rain (drizzle) had a 
higher ammonia content than ordinary rain. 
Later also HANSEN (131) pointed out that a 
long continous rain falling at a low rate had a 
higher ammonia content than showers. 


6. Origin of ammonia nitrogen in rain 


a. The equilibrium in the air—sea water system 


In connection with his fundamental plant 
nutrition experiments BOUSSINGAULT suggested 
that the sea was the source for ammonia in 
air. His theory was strongly supported by 
SCHLOSING who in a series of papers (257— 
265) discussed this possibility on the basis of 
his own experimental data. In an introductory 
paper he outlined the theory as follows: As 
all plants take up nitrogen only in combined 
form (nitrogen fixing bacteria were not known 
that time) and some loss of combined nitrogen 
occur during the decay processes of plants and 
animals, the total reserve of combined nitrogen 
in air must decrease unless there are some 
processes by which it is brought to the air or 
converted to combined form. It is known that 
electric discharges in the air oxidize nitrogen 
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to nitrate which is carried down in rainwater. 
This amount can, however, not account for 
the nitrogen taken up by the plants and lost in 
drain water; there must exist a big reservoir 
which can deliver combined nitrogen to the 
air. This reservoir is the sea which contains 
relatively large amounts of ammonia. The 
ammonia cycle is then this: The nitrate in 
drainage water is carried to the sea where it 
becomes reduced to ammonia which diffuses 
into the air and is either absorbed from the 
air by the soil or even the plants or is washed 
down by rain. The losses which occur (free 
nitrogen) during this cycle are compensated 
by formation of nitrate in electrical discharges. 

SCHLÖSING tried to verify his theories experi- 
mentally. He determined the equilibrium 
constant in the air and at different temperatures. 
Further he analysed air for ammonia and from 
his results he draw the conclusion that his 
theory was correct. He also studied the absorp- 
tion of ammonia by the soil and estimated 
that about 40 kg ammonia nitrogen was 
absorbed from air per year and hectare. Further 
he found that rainfall did not always make 
the ammonia content of air to decrease. Some- 
time it even increased after a rainfall all de- 
pending on how the equilibrium conditions 
in the place where the rain was formed 
differed from the equilibrium conditions near 
the surface. 

His theory was, however, not accepted 
everywhere. An investigation from Spain (10) 
indicated that the ammonia content in rain 
was in general lower when the wind came 
from the sea and the same was observed also 
by Posnyakov (239) in Russia. In Great 
Britain MILIER (214) started an investigation 
on rainwater collected at a number of seaside 
places and isles on the Scottish west coast and 
found in general low ammonia contents (see 
table 1), Similar results have been obtained 
later by Braapu (51), HANSEN (131), and 
LeercANG (177). Unless the higher inland 
figures are due to a greater portion of cyclic 
dust present it is obvious that these observa- 
tions do not tally with ScHLôsiNG’s theory 
about the sea as a source for ammonia in air. 

The actual ammonia distribution in the sea 
is fairly well-known. The concentration at the 
surface in the temperate and tropical zones has 
been estimated to 0.03 mg ammonia/l. (WAT- 
TENBERG (297) see also Bucx (57). It decreases 
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downwards to around 0.01 mg at 100 m 
depth and to around 0.005 mg at 200 m depth. 
Already this gradient could suggest that there 
is a steady absorption of ammonia from the 
air by the sea. But the gradient is normally 
explained in a different way. The nitrate 
shows an opposite gradient. It is practically 
zero at the surface but increases downwards 
to reach a comparatively constant value below 
200 m depth. This is due to the content of 
plankton which is greatest at the surface but 
which is limited there by the nitrate supply, 
in this way keeping a relatively high ammonia 
concentration at the surface. At greater depths 
the limiting growth factor is light; the nitrate 
level is therefore greater and that of the am- 
monia less. This organic life thus acts as a 
catalyst for the photo-chemical reduction of 
nitrate to ammonia, the equilibrium being 
determined by the light intensity. In polar 
regions and elsewhere where vertical con- 
vection takes place no such gradient is observed 
with the well-known result of an intense 
plankton production. 

From this point of view the concentration 
gradient of ammonia in the sea does not 
necessarily indicate a transport of ammonia 
from the air to the sea. 

Calculations by Bucu (57) show on the other 
hand that at equilibrium with the ammonia 
content in sea water the concentration of 
ammonia in air should be only 0.4 ug/m3 
air. Bucu, however, points outt hat reliable 
figures on the ammonia pressure at equilibrium 
have been obtained only at relatively high 
ammonia concentrations. An investigation at 
lower concentration would be desirable. 

The actual ammonia concentration in air 
has been determined in a few cases, mostly 
long time ago. The most reliable results are 
summarized in table 5. Of these investigations 
that of Levy from Montsouris covers a whole 
year with daily analyses. It is seen that the 
monthly variation is small and there was no 
clear trend between different months. In 
Budapest the values were higher in summer 
than in winter. The results from Buenos Aires 
are comparatively high probably due to 
atmospheric contamination. The figures from 
Pic du Midi (at 2,900 m altitude) are probably 
the most realiable ones existing for the am- 
monia content in uncontaminated air. The 
average of these is 13.5 ug/m3. 
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Table 5. Ammonia in air 


m 


Ref Author and year Place | ug H3N/mÿ Remarks 
179 Levy 1880 Montsouris, France 21—29 
113 Fodor 1881 Budapest, Hungary 39 average 
219 Miintz et al, 1882 Pic du Midi, summit, France 7.2—30.3 
58 Buenos Aires 1906 Buenos Aires, Argentine _ 75—87 monthly 
averages 


The ammonia concentration in air thus 
appears to be about 30 times higher than the 
equilibrium concentration calculated by Buch, 
and the sea as a source for ammonia in air 
seems to be out of question, unless it is brought 
to the air by some other means than pure 
diffusion. One such possibility is the spray 
which is known to carry chlorides and other 
salts from the sea up into the air. If the amounts 
of different salts brought into the air by spray 
were in the same proportion as in the sea, 
however, the amount of ammonia in spray 
should be very low. If the spray has a different 
composition from the sea water an enrichment 
of ammonia in the spray is possible. 

As a conclusion it can be stated that all 
evidence so far is against SCHLOsING’s theory 
about the sea as the great source for ammonia 
in air. 


b. The equilibrium in the air—soil system 


In 1872 BRETSCHNEIDER (54) published an 
account of experiments where he tried to 
estimate the amount of ammonia absorbed by 
the soil. These are probably the first experi- 
ments of this kind. He relates an earlier dis- 
covery by another German, SCHONBEIN, who 
found that H,NNO, was formed on eva- 
poration of water. The conclusions drawn by 
this discovery was that H,NNO, was syn- 
thesized in the plants during their transpiration. 
BRETSCHNEIDER calculated, however, that a 
field of red clover would receive only four 
per cent of its nitrogen in this way. (This 
was also before the discovery of the nitrogen 
fixing ability of legumes). He found, further, 
that water which had been in contact with air 
contained an excess of ammonia over nitrate 
and concluded that the excess must have been 
absorbed from the atmosphere. He then started 
to investigate the absorption ability of “syn- 
thetic” soil and found that a moist mixture 


of 95 per cent quartz sand and $ per cent 
ulmin (a humic acid prepared from sugar) 
absorbed 46 kg ammonia nitrogen/ha/year. 
Pure quartz sand absorbed only 1.6 kg and 
a water surface 0.8 g. His conclusion was 
that if SCHONBEIN’s theory was correct the 
amount of nitrogen fixed in this way was of 
much less importance than the absorption of 
ammonia from the air. Very soon it was found, 
however, that SCHONBEIN’s theory was in- 
correct, his results were due to contamination 
of the air in the laboratory where he made his 
discovery. But BReTSCHEIDER’s results were 
really important. Similar experiments usin 

acids instead of soil carried out in Rostock (138) 
in 1879—1891 gave 30.6 kg ammonia nitrogen/ 
/ha/year and in Tokyo (157) 11.6kg. SCHLÔSING 
found, as mentioned before, that the soil 
absorbed 40 kg. 

In 1911 Hatt & MILLER (129) in Rothamsted 
published results of a critical study on the . 
absorption of ammonia from air by acids. 
They were unable to show that soil could 
absorb ammonia from the atmosphere in 
direct experiments. By using two dishes filled 
with acid, mounted at different heights in the 
same place, they found, however, that the 
lower contained generally less ammonia than 
the upper one which they took as an indica- 
tion of a steady absorption of ammonia by 
the soil. To this interpretation it can be ob- 
jected that the difference in ammonia contents 
could also be due to differences in the volume 
of air brought in contact with the dishes. The 
turbulence of the air is generally increasing 
with the height above surface. The amount 
they found absorbed in the dishes were, how- 
ever, quite low. The results from 1908—09 
and 1909—10 for four different fields at 
Rothamsted are shown below, the amounts 
are all given in kg per hectare. The dishes 
were placed at 5 and 105 cm height respec- 
tively. 
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I II 


Il IV 


Field es S| Ne a 
Upper | Lower Upper | Lower Upper | Lower Upper | Lower 
I908—1909 ..... | 1.55 0.88 1.25 1.45 Tal 7 1.29 0.79 0.66 
I909—I9QI0 ..... 77 1.02 1.28 1.98 1.41 20T 0.95 0.88 


The fields II and II are manured with 
anımonium sulfate in the spring. This manure 
brought about a greater absorption of am- 
monia in the lower together with an ro-fold 
increase which persisted for some time. This 
explains why the “lower” figures are higher 
for these fields. Thus a reverse process may 
take place under certain conditions. 

The figures in the table are, however, quite 
low and they drew the conclusion that the 
amounts of ammonia absorbed from air by 
the soil is of no practical importance. They 
had, however, the dishes covered with fine 
meshed brass nets to prevent dust particles 
from entering the dishes. This arrangement, of 
course, does not allow a free contact with air 
which they also showed by comparison with 
uncovered dishes. Uncovered dishes absorbed 
about four times as much ammonia. Their 
results are thus not comparable with those 
related before. Obviously considerable amounts 
of ammonia can be absorbed from the atmos- 
phere by soils under suitable conditions. 

That a reverse process may take place has 
also been observed by EHRENBERG (96) and 
SCHNEIDEWIND (266). They showed that lime- 
rich soils, when receiving ammonia in fer- 
tilizers, may lose ammonia to the atmosphere. 

Any moist soil, with a pH lower than 7 will, 
however, absorb ammonia from air, the 
amount depending upon the ammonia content 
of the air, and the rate of vertical mixing, i. e. 
the volume of air brought into contact with 
unit soil surface. In this connection some recent 
results given by INGHAM (148) indicate that 
plant leaves can adsorb ammonia physically, 
this ammonia then being washed down to the 


soil by rain. 


c. The source of ammonia in air 


In their summary on the Rothamsted rain- 
water investigations Russet, & RICHARD (254) 
in 1919 listed the probable sources of ammonia 
in air as the sea, the soil and fossil deposits in 
which ammonia is released by combustion. 


They thought that the soil was probably the 
most important source, this conclusion being 
based mainly on the work of Hatt & MILLER 
(loc. cit.) and the fact that, at least at inland 
places, the ammonia content in rain is higher 
during summer than during winter. The work 
by Hırr & Miıtrer has been discussed in the 
preceding section and if the sea, with a pH 
of about 8, absorbs ammonia from the air 
most soils will certainly do the same. The fact 
that the ammonia contents are higher during 
summer than during winter at inland places 
may be due to increased convection in the air 
during summer, the excess over winter merely 
reflecting the contribution of cyclic dust 
during the hot season. Thus the sea and the 
soil cannot act as a real source for ammonia 
tal Abe 

The third source, namely ammonia in 
fossils used as fuel is considerable. Already 
RussELL & RICHARD calculated that in Great 
Britain the amount of ammonia nitrogen 
released by combustion of coal was about 5 
times higher than that actually received from 
rain water but they did not think it could 
be the only source. It is well known that the 
ammonia in rainwater in industrial centres 
is high but it generally decreases rapidly going 
out from the centre. Due to the general 
turbulence and movement of air masses the 
ammonia delivered from local spots will be 
rapidly spread out to cover extensive areas. 
To this source should be added the role of 
forest fires, as yet unexplored. As a matter of 
fact these sources of ammonia are the only 
ones possible as far as we know. Maybe in 
future some physico-chemical processes of 
ammonia formation in nature can be dis- 
covered, e.g. photo-chemical reduction of 
nitrate or even nitrogen to ammonia. So far, 
however, no processes of that kind seem to 
be known. 

The higher European values in fig. 2 com- 
pared to non-European are most probably 
due to a higher coal consumption in Europe 
than in other countries. 
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To summarize, it can be stated that the only 
net source of ammonia in air seems to be the 
ammonia released on combustion of fuels. 
Ammonia in cyclic dust adsorbed on fine 
particles from the ground does not increase 
the ammonia pressure of the air and does not 
contribute anything new to the soil when 
washed down by rain. 


7. The origin of nitrate and nitrite 


The formation of nitric oxide from nitrogen 
and oxygen at high temperatures has been 
known for a long time. It was therefore quite 
natural to assume that nitrate in rain was 
derived entirely from electric discharges in the 
atmosphere especially as many earlier in- 
vestigations reported increasing nitrate con- 
tents in rainwater during thunderstorms. Later 
there was, however, some doubt about this 
point; POSNJAKOV (239) for instance, failed to 
find any correlation between nitrate content 
and frequency of lightning; the same was 
also pointed out by Ram (243). Others (110, 
111, 164) found that both ammonia and 
nitrate increased during thunderstorms. In one 
place with exceptionally high nitrate content, 
namely Hanoi, Tonkin (17) there was a 
positive correlation but these values are so 
extremely high that they may be classified as 
“very improbable”. 

HUTCHINGTON (145) has, on basis of the 
analyses from Oklahoma (110, 111) estimated 
the amount of nitrate which could be derived 
from electric discharges in this place and found 
around 14 per cent could be accounted for. 
The rest must be derived from some other 
source. This low estimate is compatible with 
the fact that in tropical places with a very high 
frequency of electric discharges the total 
amount of nitrate nitrogen is not exceedingly 
high and not at all in proportion to the 
frequency of lightning. 

Some workers from India, namely Rao, 
GAPOLO & Duar (244), Duar & Ram (89) and 
RaM (243) have suggested another possible 
source of nitrate in the atmosphere. On ex- 
amining the relative amounts of ammonia and 
nitrate in rainwater in temperate and tropical 
places it was found (see Mier [213]) that 
the relative amount of nitrate nitrogen is 
higher in tropical places although the total 
nitrogen is of the same order. This they 
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explained to be due to a photo-chemical 
oxidation of ammonia to nitrate in tropical 
places, the intense ultraviolet radiation deliv- 
ering the necessary activation energy. Accord- 
ing to Ram (243) the ammonia molecule should 
be oxidized by O, possibly through the 
formation of ozone, O, from O, by the action 
of ultraviolet radiation, the ozone being a 
strong oxidant. 

Several objections can be raised against this 
theory. Firstly it is doubtful whether the 
nitrate fraction is higher in tropical places than 
in temperate ones. This has been discussed 
earlier in this paper (see p. 222). Secondly the 
amount of ultraviolet radiation received near 
the ground is highest at the poles, and decreases 
towards the equator. 

There seems to be, however, a positive 
correlation between nitrate in rain and ozone 
content in air. This has at least been shown 
qualitively in 1876 by EUGLING (103). This 
could, of course, be taken as an indication. that 
the nitric oxide formation in air is facilitated 
by the presense of ozone. But on the other 
hand the appearance of ozone at lower alti- 
tudes indicates great vertical convections in 
the upper atmosphere as suggested by REGENER 
(245). If nitrate is formed at high altitudes 
such great convections will bring fresh reserves 


of nitrate into the lower atmosphere as well, . 


where it could easily influence the nitrate 
content of rainwater. So it does not seem 
necessary to assume that ozone itself takes 
an active part of the nitrate formation. 

In the atmosphere there seems to be a consid- 
erable concentration of nitrous oxide. It was 
first discovered by ADEL (2) in 1939 from its 
absorption spectra in the infra red. Later it 
was estimated by the masspectrographic 
method on samples of air taken near the 
ground to 5-107 per cent by volyme. 
(SLOBOD and Krocıı [276] in 1950). On the 
basis of these and some other results ADEL 
(3, 4) proposes a nitrogen cycle in which 
nitrous Rs takes part. He suggests that N,O 
is formed in the soil by decomposition of 
H,NNO,. As N,O is a rather inert gas it 
escapes into the soil air from where it diffuses 
into the atmosphere. When it reaches high 
altitudes it is decomposed by the action of 
ultraviolett light (A < 2,000 À) into N,, O, 
and NO. The nitric oxide thus formed is, 
however, also photo-chemically decomposed 
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into N; and O,. The nitrogen formed diffuses 
into the lower atmosphere where it is recom- 
bined with ©, to NO under influence of 
electric discharges. The nitric oxide formed 
in this process is brought to the soil by pre- 
cipitation where it combines with ammonia 
to H,NNO, whereby the cycle is closed. 

Several objections can be raised against this 
theory. The decomposition of H,NNO, into 
N,O and water is well known. For the pure, 
dry salt, however, it does not occur below 
180° C. The reaction might, of course, still 
occur in the soils at ordinary temperatures 
as the result of some now unknown micro- 
biological action, especially as the energy 
released during the decomposition is con- 
siderable. The presence of H,NNO, in the 
soil is, however, also unlikely. The ammonia 
content in soils is generally very low, much 
lower than the nitrate as it is rapidly oxidized 
to nitrate as soon it is released on decay of 
organic matter. 

It is more probable that the N,O of the 
atmosphere is formed from nitrogen com- 
pounds in fossils in their combustion. 

The decomposition of N,O under the 
influence of ultraviolet radiation is known to 
occur, at least in presence of oxygen, when 
NO is formed. If the decomposition stopped 
at that stage the NO thus formed could well 
account for the nitric oxides present in the 
atmosphere. 

That nitric oxide exists to fairly high alti- 
tudes can be inferred from an Australian in- 
vestigation in 1914, ANDERSON (11,12), MASON 
(202, 203) et al. They made daily analyses of 
nitrate and nitrite in rainwater near Mel- 
bourne and found that the content and 
amount of these compounds was dependent 
on the origin of the air masses. They grouped 
the values according to the weather type into 
the following groups 


1. antarctic low pressures 

2. antarctic low pressures influenced by trop- 
ical low pressures 

3. tropical low pressures. 


It was found that each individual rainfall 
of types 2 and 3 carried down an amount of 
nitrate and nitrite which was independent of 
the amount of rainfall. The tropical low 
pressure rains gave the greatest amount of 
nitric nitrogen. 
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Further analyses of the air for nitric nitrogen 
confırmed their assumption that the amount 
of nitric nitrogen in each rainfall was to a 
great extent dependent upon the vertical 
extension of the air masses involved. From 
their air analyses they calculated the height 
required to give an amount corresponding to 
a certain type of rainfall. The calculated height 
was of the same order of magnitude as that 
estimated by meteorologists. They also found 
that the nitric nitrogen content of air was much 
higher on the forward side of a low pressure 
system than on the rear side. 

From their results it appears that nitric 
nitrogen should be fairly evenly distributed 
throughout the troposphere. 


8. The geochemistry of nitrogen 


Geochemical data on nitrogen gives an 
interesting picture of its distribution in nature. 
Some data are given by HUTCHINGTON (loc. 
cit.). The nitrogen content of igneous rocks 
is surprisingly constant, roughly 0.005 per cent 
(by weight). Most of it is present as ammonia 
while nitrate seems to be absent. Sedimentary 
rocks contain as an average 0.051 per cent, 
0.046 per cent thus being biochemically fixed. 
The total amount present is thus 


Igneous rocks 37,000 g N per cm? 
Sedimentary rocks 87 a DD 
In the atmosphere TSS dh Digan 


Thus it is seen that roughly 10 per cent of 
the nitrogen content in the atmosphere has 
been fixed by the biosphere during the time 
since life appeared. Part of this nitrogen is 
being released at the present time through 
combustion of coal but this amount is prob- 
ably small. The fixation process is, of course, 
still going on. 


9. The biological importance of combined 
nitrogen in air and precipitation 


For most plant species a supply of combined 
nitrogen is essential for development, and of 
course, they do not care about the origin of 
this. The rate of delivery is the only relevant 
question for the plants. 

It may be of some advantage to classify all 
the combined nitrogen brought to the soil 
in two groups, namely cyclic and non-cyclic, 
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the term cyclic taken in a short time sense. 
To the cyclic group belongs all nitrogen 
precipitated from the atmosphere as cyclic 
dust and nothing more. To the non-cyclic 
group belongs all other forms of combined 
nitrogen, even that released on combustion 
of coal. 

The cyclic dust does not contribute anything 
new to the total land surface, on the contrary 
some of it may be carried out over the sea 
and precipitated there. This fraction represents 
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a loss. The biological importance of cyclic dust 
lies in the fact that it will redistribute non- 
cyclic nitrogen and even out differences in the 
soil income of non-cyclic nitrogen. 

As cyclic dust will be represented in rain- 
water analyses, these may well represent the 
total income from atmospheric precipitation 
but are not a reliable figure for the net income. 
The possibility of a direct absorption of am- 
monia from the air by the soil must, however, 
also be taken into consideration. 


(To be continued in the next issue.) 


Electron-Microscope Study of Cloud and Fog Nuclei 


By S. OGIWARA and T. OKITA 
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(Manuscript received $ March 1952) 


Abstract 


Droplets of clouds on a mountain and of fog in an urban area were captured and the form, 
nature and size of their nuclei were studied by means of an electron-microscope and by a cham- 
ber of constant humidity. These nuclei have similar form and nature to the hygroscopic particles 
in haze and to the artificially produced combustion particles. No sea-salt nuclei were found in 
our observations, therefore sea-spray appears to be an insignificant source of condensation 
nuclei. It was found that both the cloud and the fog nuclei originated in combustion products 
which were the mixture of hygroscopic and non-hygroscopic substances, and that the greater 
part of the nuclei did not contain pure sulfuric acid. 


I. Introduction 


The nature and origin of cloud and fog 
nuclei have been studied by many investiga- 
tors and SIMPSON (1941) concluded that “sea- 
salt, sulfuric acid and nitrous acid drops were 
important nuclei in the atmosphere”. How- 
ever, these nuclei are so small that their size 
and nature cannot be studied by an op- 
tical microscope. Recently Kurorwa (1951) 
observed the sea fog nuclei by an electron 


microscope and found that many of them. 


were solid particles insoluble in water and 
that only a small fraction of them were sea- 
salt particles. No sea-salt nuclei were found 
on Mt. Niseko (1,300 m) by the same author 
(Kurorwa and TADANO 1951). However, he 
did not refer to the origin of the nuclei. 
From May to July, 1951 we observed the 
cloud nuclei on Mt. Zao and the fog nuclei in 
. Sendai (fig. 1). Observations and experiments 
were made mainly by an electron microscope 
and by a chamber of constant humidity. To 
study the nature and origin of nuclei we 


compared them with haze particles, sea-salt 
particles and artificially produced combustion 
particles. 


II. Observation of cloud and fog nuclei 


1. Capture of cloud and fog droplets 


Cloud and droplets were captured on a 
thin collodion film coated on the sample 
holder (Trager) of the electron microscope. 
As the droplets thus captured on the Trager 
evaporated quite rapidly the following tech- 
nique was employed: 

A small metal case was placed on the stage 
of an optical microscope as shown in fig. 2 (a) 
and the Trager was put into it. Drawing the 
air by a syringe which was connected to the 
case, the droplets in the air were captured 
upon the collodion film and their locations on 
the film were easily detected by the micro- 
scope. 


~ 
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140° ° optical microscope. The humidity in the room 
he was 60 to 70 % and the magnification of the 


de microscope was 900. Some of the nuclei were 


too small to be detected by the optical micro- 
scope. 

Then the nuclei were put into a small 
JAPAN SEA chamber of constant humidity shown in fig. 

2 (b). The chamber was a tapered glass vessel 
whose diameter was about 4 cm, and on its 
bottom a sheet of filter paper soaked with 
distilled water was placed. In the chamber 
almost all nuclei increased their diameter by 
absorbing the water vapour and the incre- 
ment continued for 3 or 4 minutes. The final 
diameter was about 2 to 3 times as large as 
the initial one. On the small nuclei, which 
could not be detected at the room humidity, 
water droplets were also formed. Two nuclei, 
however, retained their original size in the 
chamber. 

The temperature of the filter paper was 
equal to room temperature and the experi- 
ment was made in an underground dark room 
to avoid temperature fluctuations. Therefore 
the humidity in the chamber will be consi- 
dered as just 100 %. As the nuclei which grew 
in the chamber are able to form cloud and fog 
droplets in the atmosphere with a humidity 


Fig. 1. 


The total number of droplets captured 
were 60 on Mt. Zao and 15 in Sendai, respec- 
tively. The time and station of observation 
and the weather conditions are shown in 


Table r. 


2. Nature of nuclei. 


The size of the nuclei which remained on 
the collodion film after the cloud and fog 
droplets evaporated was measured by the 


Table 1. Observation of cloud and fog. 


Date ana Station of Observation Try Twet | R.H. ae; Wand ECS 
m (26) CC) (92) (m) Direction Form 
11 May | 03 30 | Tohoku University....... 8.8 8.8 | 100 — = Fog 
5 June | 07 oo » > Mine 17.4 1740, 98 250 — » 
15 » 04 00 » a a Sirota 27.540173 98 250 — » 
17» 05 00 » RTS 15.4 | 15.4 | 100 250 — » 
19 July | 05 50 | Mt. Zao (Hütte 1,400 m).| 13.0 | 13.0 | 100 200 WSW Cu 


23 » \ 08 oo | Mt. Zao (Mt. Katta 1,750 mi = » » 
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of 100 % or below we shall refer to them as 
“hygroscopic nuclei”. The two nuclei referred 
to above are incapable of acting as condensa- 
tion nuclei except in the supersaturated atmos- 
phere, shall be referred to as “non-hygrosco- 
pic nuclei”. The above results show that cloud 
and fog nuclei are formed in the atmosphere 
with a humidity of too % or below. There 
was one cloud sample in which no nuclei 
were found either by the optical or by electron 
microscope. 

After the above experiment was finished, 
the form and size of the nuclei were observed 
by the electron microscope. Miscellaneous 
forms of the nuclei were observed and typical 
forms of them are shown in the photographs 
on page 236 and 237. Those from (1) to 
(13) are hygroscopic nuclei and (14) is a non- 
hygroscopic one. 


(1) ... Nearly circular thin filmy substance, 
in about half of this type of nuclei 1 to 10 
small dark spots are also seen. These nuclei 
range in diameter from 0.3 to I u. 

(3) and (6)... A large well defined dark 
particle and thin filmy substance around it. 
The former is of polygonal, circular or 
ragged form and its diameter is 0.3 to 1.2 u. 
(8)... Special form as seen in the photo- 
graph. 

(9) . .. Spotted thin substances, in the centre 
of which a dark particle is seen. 

(to) sana (2) ee. Resemble (3) etc... As 
will be described later, however, they do not 
readily lose their hygroscopic nature by 
electron-microphotographing. 

(14)... Non-hygroscopic nucleus. Presum- 
ably it is a clump of carbon particles of 
oil as shown in (15). It should be noted that 
most nuclei belong to the former two types. 


After the electron microphotographs of the 
specimens were taken, the hygroscopic nuclei 
were again put into the chamber of constant 
humidity. Then it was found unexpectedly 
that the moisture in the chamber did not 
condense on the greater part of the nuclei, 
that is to say, these nuclei had then lost 
their hygroscopic nature. It is supposed that 
by electron bombardment the hygroscopic 
substances contained in the nuclei were evap- 
orated during photographing. Examples of 
these nuclei are shown in photos (1), (3), (6), 
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(8) and (9). On the other hand three nuclei 
did not readily lose their hygroscopic nature, 
two of which are shown in (10) and (12). It 
should be noted that whether the nuclei lost their 
hygroscopic nature or not did not primarily 
depend upon their size, since, even some small 
nuclei did not necessarily lose their hygroscopic 
nature, while some nuclei with large diameters 
lost it by photographing. These results may 
be explained by assuming that generally several 
kinds of hygroscopic substances are contained 
in the nuclei and the rate of evaporation of 
these substances by electron bombardment 
differs from one kind to another. 

To see whether the residues of the nuclei 
which lost their hygroscopic nature were 
soluble in water the following experiment was 
made. The residues were placed in a sealed 
glass vessel lined with moistened filter-paper 
as shown in fig. 2 (c). When the temperature 
of the filter paper was higher than the room 
temperature the moisture in the vessel con- 
densed on the residues and water droplets 
were seen by the optical microscope. After 
the residues were wetted with water like this 
and then dried, their photographs were again 
taken by the electronmicroscope. 

Comparing the image of the nuclei before 
and after wetting, it was found that most 
nuclei retained their original form as shown 
in (1) and (2), (3) and (4), (6) and (7). In 
three nuclei which retained their hygroscopic 
nature, however, the filmy substance changed 
form considerably, as shown in (10) and (11), 
(12) and (13). 

From the experimental study described 
above it may be concluded that cloud and fog 
nuclei are composed of filmy substances and 
non-hygroscopic solid particles. A part of 
the substances with thin filmy appearance are 
also non-hygroscopic and insoluble in water, 
while others are hygroscopic and soluble in 
water. 


III, Observation of haze particles 


1. Capture of particles 


For comparison with cloud and fog nuclei, 
haze particles were captured in Sendai. Near 
the ground there were large particles exceed- 
ing several microns in diameter. However, 
most particles captured on the roof, about 
10 m above the ground, were comparable in 


230 S:OGIWARA: ANDITJOKITAY 17 


ue 
\ 
à: \ 4 
« if. 
(1) Hygroscopic nucleus, (2) Nucleus (1) after rewetting. (3) Hygroscopic nucleus, Mt. Zaë. 


Mt. Zao. 


\ . j 
\ = 
m— | 
N 7 J > 4 
(5) Nucleus (4) after reaction (6) Hygroscopic nucleus, Sendai. 
with HF. 


| | | | a. | 
| we 
/ 


Lo Pe Fe > 
(7) Nucleus (6) after rewetting. (8) Hygroscopic nucleus, (9) Hygroscopic the 
Mt. Zaë. Sendai 


pr 

ed - Ne de + 

(10) Hygroscopic nucleus, Mt. Zaë. (12) Hygroscopic nucleus, 
Mt. Zao, 


(11) Nucleus (10) after rewetting 
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(16) Non-hygroscopic particles 


in haze. 
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(19) Hygroscopic particle in haze. 


(22) Smoke particle of coal 
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(14) Non-hygroscopic nucleus, (15) Smoke particles of light oil. 
Sendai. 


(17) Hygroscopic particle 
in haze. 


(18) Hygroscopic particle in haze 


(20) Hygroscopic particle (21) Spike particle of 
in haze Japanese cypress. 
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(24) Smoke particle of lignite. 
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(23) Smoke particle of grass. 
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Table 2. Observation of haze (1951). 
EEE 2225 a een ee 


Number of Particles 


Re Visi- Cloud 


Date ee Pe INoachyers Tary 2 bility | Wind 
h m woe u CORE Cs) (m) Amount | Form 

Se eee ee ee ee eS 
19 June |o5 30 I 13 

… In 5 x ne 40 18.0 | 14.3 67 — | Calm 9 as, ci, CC 
20) 06 30 3 14 79:88 0 69 = » 8 ac, as, SC, ci 
PR i 05 40 9 16 17.50 lad 72 2,500 » T CS ei en 

24 » .|o5 40 14 18 15.5 14.0 85 BEKO) |) —% o — 
2) 06 00 35 35 17.4 15.8 80 2,000 | NW 5 CSC 

4 Sep. | 06 40 41 154 1023110060 79 4,000 | WNW 8 ac, as, cu 

9 » 07 00 10 37 20:90 27.8.7 81 15000) == Io sc 

231001. 000 I 47 6.3 84 1,200 | NNW 9 ES (acs SC 


size with the cloud and fog nuclei. To capture 
these small particles we used a collector shown 
in fig. 2 (d), in which the particles passed 
through a metal tube M at high speed and 
struck the collodion film, to which they then 
adhered. 

The time and the weather conditions are 
shown in Table 2 together with the number 
of haze particles captured by the collector. 
It will be seen from the table that the haze is 
composed of hygroscopic and non-hygrosco- 
pic particles in various proportions. Experi- 
mental procedure to decide whether the par- 
ticles were hygroscopic was quite the same 
as that for cloud and fog nuclei. 


2. Nature of haze particles 


Electron-microphotographs of hygroscopic 
haze particles are shown in (17) to (20). As 
will be seen by comparing (17) with (1), 
(18) with (6), and (19) with (8) they are 
similar in form to the hygroscopic nuclei of 
clouds and fog. Furthermore these particles 
lost their hygroscopic nature by photographing. 
Thus hygroscopic particles in haze are similar 
both in form and in nature to cloud and fog 
nuclei. 

Two hygroscopic particles evaporated com- 
pletely as a result of the electron bombard- 
ment and left nothing on the collodion film. 
These particles probably consisted entirely 
of hygroscopic substances. 

One example of non-hygroscopic particles 
in haze is shown in photo (16). There are 
three particles one of which is similar in form 
to the smoke particles of fuel (?) oil and 
others have circular form. Since these par- 
ticles remained unchanged under exposure to 


hydrogen fluoride vapour it is presumed that 
they consisted of carbon. 


IV. Origin of cloud and fog nuclei 


1. Comparison with sea-salt particles 


Is has been considered that a large number 
of sea-salt particles are supplied from the sea- 
spray and that they are the most important 
condensation nuclei in the atmosphere. To 
decide whether the cloud and fog nuclei 
captured by us were sea-salt particles, small 
sea-salt particles with a diameter of about 
0.5 4 were produced by spouting a water 
solution of sea-salt. Following the same ex- 
perimental procedure as described before it 
was found that these salt particles had images 
of comparatively well defined crystalline 
form. Furthermore it was found that they 
were not only soluble in water but also did 
not lose their hygroscopic nature by electron- 
microphotographing. Thus the sea-salt parti- 
cles have different form and nature from the 
greater part of cloud and fog nuclei. 

As already described, however, a few nuclei 
of cloud and fog did not lose their hygroscopic 
nature by electron bombardment. But hygro- 
scopic substances contained in these nuclei 
were composed of those with thin filmy 
appearance, which were far from salt crystals. 

Recently Kurorwa and TADANO (Kurorwa 
1951, KUROIWA and TADANO 1951) observed the 
nuclei of clouds on Mt. Niseko and those of 
sea fog on the coast of Hokkaido by using 
an electron microscope. They found no sea- 
Salt nuclei on Mt. Niseko, and even in the 
sea fog only a small fraction of the nuclei 
were sea-salt particles. These results lead one 
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to the conclusion that sea-spray is not an 
important source of condensation nuclei of 
cloud and fog, at least not in these cases. 

It is well known, however, that a consider- 
able amount of sodium chloride is contained 
in cloud and fog water (KOHLER 1922, HoucH- 
TON and RADFORD 1938, MIYAKE 1948), and 
the sea-salt nucleus theory has been based on 
this fact. Therefore it is necessary to explain 
the origin of sodium chloride in cloud and 
fog water. From some preliminary experi- 
ments sea-salt particles were found to have 
equal hygroscopic activity to combustion 
particles of the same size. Observations of 
sea-salt particles on the coast of Honshu (Ocı- 
WARA and Oxia 1951) showed that the num- 
ber of particles larger than cloud and fog 
nuclei was about $ per c.c. and was far less 
than that of cloud and fog droplets in the 
atmosphere. Therefore it seems reasonable to 
conclude that the greater part of cloud and 
fog droplets are formed on combustion nuclei, 
and that small sea-salt particles cannot act as 
condensation nuclei in the atmosphere. When 
cloud or fog water is collected by wire screen 
or other kinds of collectors, small sea-salt 
particles floating in the atmosphere should be 
captured together with the cloud or fog 
droplets. Large particles, which would be- 
come condensation nuclei, but are not detected 
in the cloud and fog droplets owing to the 
smallness of their number, would be added to 
the sodium chloride content in the cloud 
and fog water. Quantitative explanation for 
the sodium chloride content requires further 
investigation. 


2. Comparison with combustion products 


As it was found that cloud and fog nuclei 
are not sea-salt particles but are similar in 
form and nature to hygroscopic haze par- 
ticles, the most probable source of the nuclei 
is considered to be combustion. Smoke par- 
ticles produced by combustion of many kinds 
of fuel, such as coal, lignite, leaf of Japanese 
cypress, pine tree and grass were observed 
following the method described in chapter II. 
Some of these smoke particles are shown in 
photos (21) to (24). It was found that they 
too are a mixture of hygroscopic and non- 
hygroscopic substances and that they may 
lose their hygroscopic nature by electron- 
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microphotographing. Cloud and fog nuclei 


are therefore similar in form and nature to 
combustion products. From the results describ- 
ed it may be concluded that most nuclei of 
cloud and fog originate in combustion. 

It remains a question, whether the non- 
hygroscopic particles contained in cloud or 
og nuclei are those captured by the cloud or 
fog droplets which have formed on the 
hygroscopic nuclei. As cloud and fog nuclei 
have similar form, nature and size to combus- 
tion products, the greater part of these non- 
hygroscopic particles are considered to have 
been adhering to the hygroscopic substances 
from the beginning. 


3. Chemical composition of combustion nuclei 


Generally combustion products are com- 
posed of many kinds of substances, therefore 
it is not easy to study their chemical compo- 
sition in detail. 

It has been considered that sulfur dioxide 
produced by combustion of coal etc. becomes 
sulfuric acid droplets by oxidation, and that 
these hygroscopic droplets can act as effective 
condensation nuclei of cloud and fog. To 
decide whether the hygroscopic substance in 
the nuclei is sulfuric acid the following experi- 
ments were made. 

When a droplet of sulfuric acid with a dia- 
meter of about 0.5 4 was put into the chamber 
of constant humidity it increased its diameter 
by 2 to 3 times. Thus the increment of the 
droplet was of quite the same order as that 
of cloud and fog nuclei. Therefore the quan- 
tity of hygroscopic substances in cloud and 
fog nuclei cannot be considered to be far 
smaller than that of sulfutic acid. The droplet 
was then exposed to the electron beam of 
the microscope. Even after the electron bom- 
bardment, it was found that the droplet 
did not change its size nor its hygroscopic 
nature. From the above results it is obvious 
that the greater part of cloud and fog nuclei 
do not contain sulfuric acid. Chemical engi- 
neers have found that many kinds of sulfuric 
compound are contained in the smoke of 
coal and that some of them are hygroscopic. 
Therefore it may be considered that the sul- 
furic ions dissolved in cloud and fog water 
would originate in these sulfuric compounds 
but not in pure sulfuric acid droplets. 
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Next, to study the chemical composition of 
the nonhygroscopic solid particles in the 
nuclei, they were exposed to the vapour of 
hydrogen fluoride. Some particles were attack- 
ed by the vapour as shown in photo (5), while 
others were unchanged. By the same tech- 
nique it was found that the solid particles in 
the smoke of light oil, leaf of Japanese cypress 
and grass were proof against hydrogen fluo- 
ride, while some of those from smoke of 
coal and lignite were attacked by hydrogen 
fluoride. It is presumed that the former is a 
pure carbon particle and the latter is a silicate 
or carbonate, and that the solid particles in 
the nuclei are composed of these substances. 


V. Conclusion 


From the investigations set forth above it is 
concluded that the greater part of cloud and 
fog nuclei are composed of hygroscopic and 
non-hygroscopic substances which originate 
in combustion and that sea-spray is not a 
significant source of condensation nuclei. It 
was also shown that the greater part of the 
hygroscopic substances contained in the nuclei 
are not pure sulfuric acid. It is important to 
bear in mind that, contrary to the opinion 
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which has been vigorously advocated so far, 
the greater part of cloud and fog nuclei are 
composed of neither sea-salt nor sulfuric acid. 

A question will arise: If the combustion 
products are the most important nuclei, what 
was the most effective nucleus before man used 
fire? Even at that time combustion nuclei 
would be supplied by the eruption of volcanos 
or by forest fires. Or it may be more probable 
to consider that cloud and fog droplets would 
be formed on other kinds of less effective 
nuclei such as small sea-salt nuclei at a slightly 
higher atmospheric humidity than at present. 
As it has been assumed that by the appearance 
of life the chemical composition of the air 
has changed its proportions, it is perhaps no 
too bold to assume that by the appearance o 
man the most effective nucleus had changed 
from sea-salt to combustion products. 
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Sunspot Areas, Flares and Filaments observed at the Stockholm 


Observatory in the Years 1950 and 1951 


By Y. OHMAN and L. O. LODEN 
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(Manuscript received 15 June 1952) 


Abstract 


Tables are presented giving sunspot areas, flares and filaments according to the solar 
routine observations carried out at the Stockholm Observatory in 1950 and 1951. The 
material has been distributed before as current reports to a number of Swedish institutions. 
A brief discussion is given and a description is made of the methods used for the ann- 


ouncements of the general solar activity. 


A. Sunspots and faculae. In an earlier publi- 
cation! a description has been given of the 
routine solar work started at the Stockholm 
Observatory in 1947 in order to furnish 
Swedish geophysical institutions and the 
Royal Swedish Board of Telegraph with solar 
data in addition to the information given in 
the Ursigrams. The photographic observations 
of sunspots have been made by us in the same 
way as before, that is by using a Zeiss objective 
with an aperture of only 30 mm and a focal 
length of 800 mm and by enlarging the solar 
image to a diameter of 56 mm. A passband of 
200 A is selected in the H and K region of 
the spectrum by means of a Schott and Gen. 
UG 2 filter with a thickness of 1 mm. The 
photographs are made on Ilford Lantern 
contrasting plates. 

The solar photographs obtained in this way 
have a remarkably uniform quality as the small 
aperture reduces the influence of variations in 
seeing-conditions. Moreover the ultraviolet 


ı Tellus, 2, 255, 1950. (Stockholms Observatorium 
Meddelande No. 65). 


7—203045 


filter brings out a high contrast. The spots are 
very dark and outstanding objects and the 
bright faculae can be traced not only close to 
the limb but all over the surface of the sun. In 
fact, the accumulated line absorption in the 
violet seems to produce the characteristics of 
the low chromosphere though with less con- 
trast than on real spectroheliograms obtained 
in the metallic lines and in the outer wings of 
the H and K lines. 

In table I and II a presentation is given of 
sunspot areas for 1950 and 1951 according to 
our measurements. These tables are arranged 
in the same way as those for 1948 and 1949 
given in the previous publication‘. The 
sunspot areas are thus measured in 1075 of the 
half sphere of the sun as unit and the figures 
refer to the first day-sector after the sun’s 
central meridian observed at about 8h 30™ 
U. T. The upper number in the tables 
refers to the north hemisphere and the lower 
number to the south hemisphere. When a 
number is given within parenthesis no ob- 


Wore, Coty 


~ 


Y. ÖHMAN AND L. O. LODEN 


242 


-| —o me 

—| —o [Do 

— 126) = 0 

5 1): (0) = 0 

_ — 0 — 0 

= — 0 — 0 

- | @ | —o 
=), (0) — © 

(+) (or) = = 0 
(+) (zt) u — I 
— 0 — — 0 
— 0 — — 6I 
—s — [D 0 
—0 - [D o 
— 0 — — 0 
— 0 == — 0 
ee) ©) = 6 
IE) a) (0) —= Gi 
— 0 — 0 — OI 
— 0 = 6 — 0 
— 0 — 0 (—) 0 
pase ss a) 
(=) © = == 
(—) (Sz) = = © 
=) S — — © 
(=)=0 = = € 
(—) (oz) = 1) 


(—) (0) ale) 


--— sdnos$ saAn9e uou JO + S 


— 0 — 0 Pe ie [Cease ro 
DT A0 JET (—) © (—) 0 (+) (82) 
(-)9 | —o — 0 — 0 —o [Do 
(ee) 0 — 0 — 0 + Sp | (+) o£ 
(+) (11) — 9 (+) (£1) — 0 —o — 0 
(—) (0) + 1% (—) (0) — 0 — 0 — 0 
70 ao (=) (6) 27) a (—) (o) 
"#0 9 (0) EIN © + 0 (—) (0) 
— 0 — 0 — oO (—) [0] Mi 25 ST 
en) 10 — 0 (—) 0 — 0 — 0 
(—) (0) Fa =» Ki 0 m (—) () 
(—) (0) ie ae = u (—) (2) 
—o —o — 0 — 0 + 7 — 0 
is = © — 0 — 0 — 0 — 0 
= © + 07 — oO (—) 0 — 0 — 0 
— 0 — 0 — 0 SE az Ar M 
(—) 0 — 0 — 0 — 0 — 0 (—) 0 
(—) 0 — oo + 1+ — 0 — 0 (—) 8 
at | NE EN —o |) 
— 0 — 0 — 0 (—) 0 — L (—) (0) 
— 0 + OI — 0 (—) o (—) 0 + TE 
— 0 — gi — oO ©) © (+) oz + LL 
—o | —o = & Hits) | Oops 
+ TI — 0 — 0 (—) 0 + OT (—) 0 
— 0 — 0 — 0 (—) 0 — 7 (+) (az) 
— 0 — © — © (—) of — .0I (—) (0) 
(—)1 — 0 — 0 (—) 0 + tz (—) 0 
ER To =e SR: ns —) 4 


-sısayyuased uryyIm u9AIS 91e sonjeA parerodioquy 


-0S61I ‘I °JquL 


ysıeW 


(+) (rE) 
(—) (01) 
(=) (0) 
(-) (©) 
(—) (0) 
(—) (0) 
(—) 4) 
(—) (©) 
(—) 0 
(—) 0 
(—) (0) 
(—) (©) 
— 0 
— 0 
(—) (©) 
(—) ©) 
(—) (©) 
(—) ©) 
(—) (0) 
(-) @ 
(—) (0) 
(—) ©) 
(-) © 
(-) © 


(—) (2) 
2.0) 


MO) 
(+) (©) 


‘1491 


(—) (0) 
(=3) 45) 
6) 
— 0 


(=) 
(=) Xo) 


(=) 40) 
(—) (0) 
(97020 
(+) zz 
(—) (0) 
Far) 
(—) 6) 
(—) (0) 


(=) Ie) 
(—) (0) 


VI 


(ep 


II 


OI 


dnos3 oanoe Jo suSis yyIm 1919807 uelpriolu jenyuas sung 24} 19178 10}935-Aep oy} ur seaie yodsung 


SUNSPOT AREAS, FLARES AND FILAMENTS IN 1950 AND 1951 


—o —o (—) o 
a (rm) 1, C)70) ZUNG 
Zul I.) Re) (18 
— 0 — 0 UT (=) © 
= © ==) © == 8 GeO 
— 9 No (—) (0) re 
== 6 a 0 (=) 0) RO 
(—) (0) (+) (+) + OI — I 
(—) (0) © + & —26e 
— 7 +7 == © =) © 
+ 0 == © — IT (+) 65 
GERT) + | (=) © ie 
(—) (0) = Ë (2) Gt == 
— 6 (—) 0 (=) © ae BE 
el me ey = 
Ze (—) (©) T0 Ce) 
Re (—) (0) 8 (+) (of) 
—. 6) — 0 (—) 0 — £ 
+ #E fife (Ge) 8 = 0 
+ or — 0 — 0 (+) (9) 
— 0 — SI — 0 (=) (0) 
— 0 — 0 — 0 — 0 
(0) + 08 — 7 — 0 
Jo) | =o | en 
ee) =o (20 =a) 
—o — 0 de — @ 
— 9 — 0 — 0 — 0 
— 0 + OI — 0 — 0 
a) 0% ne] EN O! BO) 
— 6 + 7 —o —o 
+ ze + £ı — 0 — 0 


| 


ENS 


| + 


— Ne mn 
ES TT En 


| | + 


na PNCRS TES 


| 


| 


+ 


| 


TO Vw 


+ 


| 
| 


re en 
NN TS ee ee 


a 


| 
| 


ne 
a par 


nn 


Y, OHMAN AND L. O. LODEN 


244 


=) ON) One sy =)" © + 9 COTE eg: re oe (>=) OD 1G) (0) en 
NONE) IE) ne = sg ee (+) TE + VI 0 EN + vi 
(—) (0) ay ee Oe Hoe AC nn me IHN Er ke) Ae EN = 
(—) (0) es en) SG a. Be Fr 6 TO HR PE eI 
(—) (0) Fe IO re =) © a ag ao at ay hog = = 
(=) @) — |(—) (0) — I (—) 0 — 6 + $8 — 0 — 0 (—) 11 — 0 = ZI 
=: — | —o | —o [Ho | —-o [oo | —o (+) en] —o — |) ©) 
— I — 0 — 0 (+) (67 — 0 (+) (Sr) — 0 (—) (0) — 0 — (—) (0) II 
) 0 — — 0 + op | (+) oz — 0 — 0 — 0 (—) (0) — 0 — — 0 
=) © — 20 =. 6 (+) zx a ai 0 0 (—) (0) 20 = = 6 OI 
(—) (o) — — 0 (—) 0 — 0 (—) (0) —o — 0 — 0 — 0 — — 0 
A) o (—) 0 — 0 (+) (S) — 0 — 0 — € — 0 -- — 0 6 
(—) (o) = at 49 RQ (—) (o) + 97 Fi (ne ME (JE € = tro 
10) = ag =o Ac (—) (o) ee peat Gd Ma (—) (0) =e) 8 
== 10 — | (—) (Lx) — 0 — 0 — 0 + 87 — 0 — 0 — 0 _ — 
— 0 — |(—) (0) — 0 +L — 0 — 0 — 0 + 8 — 0 — — L 
= ze 0 + À 0 =o CO Ga Bier ID) = = 
ve a = an FUIT ar = (=) ET GE) EI) — = 9 
(—) (0) = OF 1.0) We) 2 (0) sage Re re =. = aan 4 
(—) (2) Oe =) AO) ee) =. ae ae Er ee = = $ 
SAUER 0 pelo 0 Sane aan) AD rg EEE = = 
ON) =. no ar mee) CES + 0 ROUE =e = L 
— 0 — 0 — 0 — 0 — 0 — 0 — 0 — 0 (—} (0) EZ (—) (0) — 
—= @ — oO — 0 — 0 — 0 — 0 = © = à (—) (0) —= 6 (—) (0) — £ 
= a Fr ER = rd Ge WA) ea |B Ce ene Con Byte) 
= + 8 = 0 — EN mu) 0 GEHN, 7) US CON ENT z 
(=) (oy Bh) (Ol). 8 zu So (—) (0) = 70 (=) NIS EINE) = =? 
(+) @ |.) Cr atts =e (—) (0) ED DIOR E00) He: = | I 
22 | AON | ‘PO | das | any | Amf | sunf | AE | dy | ysıeW | "1q2J | uel | 


sısoyyussed uryyım UdAIS 91e sonjeA parefodısyuy] 


-— sdnos3 sAnse uou 10 + sdnoiS sAr9e Jo suSis yITM 1943030} uerprsaur Jesyuss S,uNg ou} 499372 103995-4ep ayy ur seaIe Jodsung 


“IS61 ‘II PGR L 


245 


[o) 


++ ++ 


a 
a 


RTE Se ER 
er 
= 

no on 


SUNSPOT AREAS, FLARES AND FILAMENTS IN 1950 AND 1951 


mn 
SR ST 
Se Ras 


° 


+ 
io} 


(=) 6) 
(=) 30) 


(+) (11) 
(—) (0) 


(—) (@) 
RO) 


==. “0) 


+ VI 


en 70 
(Ge) ee 


+ of 


=) AG) 


N) 
=) 0) 


O 
O 


O 


(—) (0) 
ES) (0) 
(—) (o) 
(—) (0) 


(=) 
(—) (2) 


Oo 


of 


67 


87 


Le 


97 


ST 


124 


£z 


CE 


he 


OT 


61 


QI 


LIT: 


91 


St 


246 Y. OHMAN 
G 
240 

| 
220; 
200) 
180 
160 
140 
120 
100 

80 
60 

<o i 

Se 
Pr 

20 L 

E 

00 726 20 60 80 100 720 740 760 wo 200 5 
Fig. 1. Comparison between Greenwich Observatory 


sunspot areas G and Stockholm Observatory sunspot 
areas S. 


servation was possible that day but the value 
has been obtained by interpolation. If this is 
the case for several days in succession the 
interpolation has been extended over several 
days. Zero means that no spot is situated in 
the day-sector after the central meridian 
and — that a long period of cloudy weather 
has made the determination impossible or very 
uncertain. 

When the numbers in table I and table II 
are followed by a + sign, this means that the 
sunspot group is surrounded by bright Hx 
faculae according to spectroheliograms ob- 
tained the same day with our patrol spectro- 
heliograph. A + sign within parenthesis means 
that no spectroheliogram has been obtained 
the same day but that the activity of the group 
is judged from spectroheliograms obtained be- 
fore or after that date. A — sign after the 
numbers means that no noticeable bright Hx 
faculae are seen on our spectroheliograms and 
a — sign within parenthesis that the low 
activity of the sunspot group is judged from 
earlier or later observations. If neither a + sign 
nor a — sign is accompanying the number 
a long period of cloudy weather has made 
even an interpolation of the activity character 
impossible. 
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The spectroheliograph used for the patrol 
work is the same as that described in the 
previous work. It has been equipped with two 
rotating Anderson prisms so that the instrument 
can be used not only as a spectroheliograph 
but also as a spectrohelioscope. 

As the sunspot areas determined by us are 
based on photographs obtained with a small 
instrument and in the. violet light it has been 
interesting to compare our determinations 
with those made with a considerably larger 
instrument at the Greenwich Observatory. By 
the great courtesy of the Astronomer Royal 
we are able to present in fig. 1 a comparison 
between our sunspot areas (S) and those 
determined at the Greenwich Observatory (G). 
For the comparison a number of well defined 
groups appearing near the sun’s central 
meridian in April and May 1948 have been 
selected. 

According to the comparison presented in 
fig. I our sunspot areas are systematically 
somewhat smaller than the Greenwich values. 
The reason for this my be the smaller scale of 
our photographs and the different wavelength 
regions used. The scattering of the individual 
points is moderate, however, which indicates 
that our determinations may be used as 
representative values of sunspot areas. In fact 
some of the greatest deviations seem to be 
due to rapid changes in the size of the spots. 

B. Solar flares. During 1950 and 1951 a 
number of flares of intensity 3 have been 
observed at the Stockholm Observatory with 
the patrol spectroheliograph. All these flares 
have been recorded after immediate announ- 
cements by the Royal Swedish Board of 
Telegraph of fade outs in the short wave 
communications. These observations are listed 
in table III together with some smaller flares 
observed by us accidentally. 

C. Filaments. A special attention has been 
given to the appearence of large filaments as 
such objects sometimes produce ionospheric 
disturbances three or four days after the transit 
over the sun’s central meridian. The large 
filaments are generally very outstanding 
objects on our He spectroheliograms. Those 
recorded by us in or near the first day-sector 
after the central meridian are listed in table IV 
and table V. The tables should be compared 
with tables I and II in order to judge which 
periods may be complete and which may not 
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, Stockholms Ob forium 
be complete. For such days when no + or Desum Nov.2 1951 Eu Pal i 
no — sign is given in table I and table II no 
information about the presence of filaments 


has been available. 


D. Polar maps for announcing the solar activity. 
ior the distribution twice weekly of reports on 
the general solar activity it has been found 
very appropriate to use polar maps of the sun 
of the type presented in fig. 2. The sun’s north 


a“ 
Oversikt av soloktiviteten 
Solen med dogsektorer sedd fran norra polen 


Table III. Flares observed at the Stockholm 
Observatory 
First Blast: 
Day | che Max he Imp | Lat |Long 


1950 hm|hm | hm 


May 3 9) SQ) (0756 | 70:28) 3 I4N IE 
SUIMTONOMIREONTO N TOs 4A 2 I2N | 29W 
6G) | 23738] 73 391 13-45 | 2 12N | 47W 
10 | 10 50] II I4| II 22 | 2—3 6N | 26E A. 
Dees 26 | To 30] To 37T | 1038| x 12N | 18E Mindre aktivitetscentrum e 
P Storre oktivitetscentrum ,,* 
1951 Teckenforklaring Flare (red datum) red. © 
4 5 
May 18 | 10 45 | 10 47 | II 19 | 3 ısN | 30W Storre filament i. 2 AU 
June 8 | 15 46 1.2 | I4N (STE Stork gron koronalinje gron a 
EURE a3 à 5 Fig. 2. Stockholm Observatory solar map for Nov. 2, 
s 27 Ss ; N 20E 1951 giving the distribution of activity centres (spots) 
SPEARS O8ilpr srs Io 


and filaments. The arrow is pointing towards the earth. 


Table IV (1950) 


Jan. | Febr. March | April | May | June | July | Aug. | Sept. | Oct. | Nov. Dee: 
5 (16) I 4 5 6 I 4 2 2 (2) (1) 
29 25 10 7 8 8 2 9 5 (4) (9) 9 

28 IS (11) Io (11) 5 II 7 5 12 (15) 
16 20 (15) IS of 12 8 (14) 14 (29) 
23 21 18 16 8 14 10 18 (17) 

30 22 25 17, 13 20 (16) 23 
(24) 18 15 24 27 
30 23 18 29 

19 

(23) 

25 

26 

29 

31 

Table V (1951) 

Jan. | Febr. | March | April | May | June July | Aug. | Sept. | Oct. | Nov. | Dee 
9 19 3 17 5 12 6 6 2 2 (2) 6 

(28) (4) (30) 14 (14) (9) 9 4 (13) 15 (14) 
(8 28 12 (12) 5 17 26 22 
(10) 18 19 (9) 22 29 (26) 

21 (26) 15 27 
30 28 29 
30 


ee eee ee eee 
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pole is used here as origin and 27 equal sectors 
are given, each one corresponding to the 
synodic rotation of the sun in 24 hours. The 
first zone is used for characterizing the solar 
phenomena north of the equator and next zone 
is used for giving the solar phenomena south 
of the equator. Different phenomena such as 
spots, flares, filaments and a strong intensity 
of the green corona line are indicated by means 
of different symbols. In this way a very clear 
picture is obtained of the whole sun and it 
is easy to see when old phenomena may be 
expected to appear again at the east limb. 

During such periods when the Stockholm 
Observatory coronagraph has been in regular 
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use the regions of high corona intensity have 
been obtained from our own recordings. For 
the information of different geophysical insti- 
tutions in Sweden we have also indicated high 
corona intensity values obtained from the 
Ursigrams. The postcard contains predictions 
for ionospheric disturbances as well. 

For our own work a big map of the same 
kind is used. This map can be turned around 
its axis so that the proper setting is obtained 
for each day. The device is in a way more 
simple than the globe used in some astronomi- 
cal institutions for the same purpose. 

Stockholm Observatory, Saltsjöbaden June, 


1952. 
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Abstract 


The paper describes the 10 m high power recorder with a rotating antenna that is used 
since May 1951 for the localization of aurorae at the Radio Wave Propagation Laboratory 
of the Kiruna Geophysical Observatory (67.8° N, 20.5° E). Continuous observations during 
the time May 1951—March 1952 have disclosed periods of auroral activity. The 
preliminary results from these observations indicate that there is a good correlation be- 
tween the auroral activity, the magnetic activity, and the appearance of the Nr-layer,.a 
special type of sporadic E ionization often appearing in connection with magnetic bays and 
supposed to be caused by the same ionizing agent as the aurora. The distribution in range 
and bearing of the recorded aurorae agrees with the simple theory that most of the radio 
wave scattering comes from those points where the radar beam is perpendicular to the 
surface of the auroral discharges. The calculated height distributions of the reflection centres 


have maxima around 120 km. 


1. Introduction 


In.order to extend the study of the general 
features of the polar ionosphere and other 
associated phenomena it was early planned 
by the Chalmers Research Laboratory of 
Electronics to study long distance scattering of 
radio waves from the Aurora Borealis (Ryp- 
BECK, 1949). Since May 1951 the Radio Wave 
Propagation Laboratory of the Kiruna Gco- 
physical Observatory (see fig. 1) has therefore 
been equipped with an aurora recorder of 
radar type operating on 10 m wavelength 
(frequency 30.3 Mc/s) using a rotating an- 
tenna. A great number of very interesting 
echoes have been received. Their distribution 
in range and bearing and the correlation with 
visually observed aurorae strongly support 


the theory that echoes are caused not by an 
ordinary type of ionosphere (ordinary or 
sporadic E) but by reflections from auroral 
discharges. . 

Scattering of radio waves from the Aurora 
Borealis has earlier been reported by HARANG 
(1940) on a frequency of 40 Mc/s and by 
Lover, CLEGG, and ELLYETT (1947) on 46 
Mc/s. During 1949 ASPINALL and HAWXKINS 
(1950) identified a number of aurora forma- 
tions in the range 400—900 km by the charac- 
teristic echoes which were produced on 72 Mc/s. 

This paper describes our 10 m aurora re- 
corder and presents the observations and the 
preliminary results obtained during the period 
May 1951—March 1952. 
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conventional self-pulsed tuned-grid-tuned-plate 
oscillator with too kW peak pulse power, 
40 ps pulse width, and a pulse repetition rate 
of so c/s. The rotating antenna, a three- 
element Yagi, is mounted on the roof of one 


IN 
NARVIK << 
ee 


\ of the observatory buildings (see fig. 2). It 


RUA, is elevated 30° above the horizontal and has 


) Sur ae = = 2 
LME LA an estimated gain of 9 db compared with that 


AER m 


of an isotropic radiator. The receiving unit of 
the recorder (fig. 3) consists of a modified 
Hallicrafter receiver, type R44/ARR-5, and a 
10 cm radar SO-13 indicator. The smallest 
signal detectable on the cathode ray tube is 
equivalent to 3 x 10-14 W at the receiver input. 
The antenna and the synchronized beam on 
the scope make two rotations per minute. The 
camera is open for exposure for 30 seconds, 


T1 SrocKHOLMNS i.e. the time of one complete revolution of 

; Foo - the antenna. It can be run automatically either 

pag continuously, taking one picture of each rota- 

Cae tion, or intermittently, taking only 12 pictures 
every full hour. 


Fig. 1. Location of Kiruna Geophysical Observatory. 


2. Recording equipment 


The equipment is principally the same as 
that used at our Rad Observatory in the 
vicinity of Gothenburg for the recording of 
meteor trails. Power is transmitted and re- 
ceived by the same antenna. The received 
echoes are shown on an intensity modulated 
cathode ray tube display of PPI-type, which 
is photographed. The transmitter is a quite 


Fig. 2. Rotating antenna of 10 m aurora recorder. Fig. 3. Receiving unit of 10 m aurora recorder. 
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Fig. 4. Theoretical radio echo from an auroral curtain. 


3. Types of echoes 

Roughly speaking, the great variety of 
auroral echoes can be classified into one or the 
other of two different types, diffuse or discrete. 
The diffuse echo must be due to reflections 
from an extended irregular surface or a large 
heterogeneous region with a random distribu- 
tion of electron density. The discrete echo, on 
the other hand, must originate from a regular 
surface containing a relatively homogeneous 
electron distribution. Therefore, as already 
done by AspINALL and Hawkins (1950) it is 
reasonable to associate a diffuse echo with an 
auroral glow or a curtain, and a discrete one 
with a streamer. In the following we will use 
the terminology, curtain, for a diffuse echo, 


and streamer, for a discrete one. Most of the 
reflected energy must come from those points 
of the curtains and streamers where the radar 
beam is perpendicular to the corresponding 
surfaces. For a certain locality the most 
favourable direction of reflection is the direc- 
tion of the magnetic north, if one assumes the 
auroral discharges to take place on surfaces 
generated by the magnetic field lines which 
pass through curves on the earth’s surface of 
constant horizontal force. 

Applying the theory of radio wave scattering 
from a turbulent medium (BOOKER, GORDON, 
1950) ‘to the case of reflection from an auroral 
curtain we get a diffuse echo on the indicator 
screen like that illustrated on fig. 4, when the 
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Ok: 16 


Fig. s. Typical radio echo from an auroral curtain. Fig. 6. Typical radio echo from an auroral streamer. 
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Fig. 7. Hourly variation of the auroral and the magnetic activities, the intensity of the Nr-reflections, 
and the appearance of polar blackouts during October—December 1951. 
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antenna rotates. We have here assumed the 
curtain to be an infinitely extended plane sheet 
with a thickness of 0.25 Ry, where R, is the 
minimum range to the sheet, and the antenna 
to have constant gain throughout an angle of 
90°. Further, we have assumed that the scale 
of turbulence is large compared with the 
wavelength, and that the mean square devia- 
tion of the dielectric constant of the layer is 
small. If, instead, the reflections had come from 
an auroral streamer, supposed to be a circular 
cylinder (with a diameter of the order of 
2 km) containing a uniform electron density 
the echo picture would have been a uni- 
formly bright, very narrow 90° arc of a 
circle. Figs. 5 and 6 show the two different 
echo types, as obtained by our Kiruna recorder 
on September 22nd, 2012 15°EMT and on 
October 16th, 2001 15°EMT respectively. The 
arrow on the top indicates the north, east is 
to the right and west to the left. The difference 
in range between two concentric circles is 


25 


Kiruna, 2202 15° EMT, Febr. 13, 1951. 


Virtual height 
km 
500 
400 
300 
200 


100 


Frequency Mc/s 


Fig. 8. Record showing simultaneous reflections from 
sporadic E-layer (100 km height) and Nı-layer 
(115 km height). 


so km and the maximum range 1,000 km 
(this could be increased if required up to 
2,000 km but is ordinarily quite sufficient). 


4. Preliminary recorded results 


During the time May 1951—March 1952 the 
aurora recorder ran day and night with the 
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Fig. 9. Day to day variation of the auroral and the magnetic activities, the intensity of the N1-reflections, 
and the appearance of polar blackouts during May, June, and July 1951. 
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Fig. 10. Day to day variation of the auroral and the 


magnetic activities, the intensity of the Nr-reflections, 


and the appearance of polar blackouts during October, November, and December 1951. 


exception of a period of relatively small auroral 
activity in August and September 1951 when 
part of the installation had to be made more 
permanent. Normally, in order to save film, 
only 12 consecutive pictures were taken during 
every hour, but in the case of special auroral 
activity, and at least once a week, each sweep 
was photographed for a period of 24 hours, 
when, also, the ordinary panoramic ionospheric 
recorder was in continuous operation. The 
hourly recordings are quite sufficient to permit 
statistical investigations of the aurora, but the 
continuous films, of course, give a better and 
a dynamic picture of the whole phenomenon, 
and are of special value for a detailed com- 
parison with the ionospheric soundings and the 
magnetic records. 

To interpret the hourly records one charac- 
teristic echo picture out of the 12 obtained has 
been analyzed with regard to the appearance 
of diffuse echoes (curtains) and discrete echoes 


(streamers). According to its strength each 
echo has been given a weighted number in 
the following way: 


= aurora of small intensity 
aurora of moderate intensity 


aurora of strong intensity 


D — 


wae, CD 


The maximum auroral activity, or frequency, 
of either kind is consequently during 24 hours 
24X3 — 72). Of course, many objections 
could be raised to such a procedure; the per- 
sonal errors are important in the reading of 
the records, no consideration has been given 
to the influence of distance on the echo inten- 
sity and so on. The method is justified only by 
the fact that it forms a basis for a preliminary 
statistical investigation of the observations. 
The variation throughout the day of the 
auroral activity, as shown by the hourly 
records, indicates a maximum activity around 
midnight and a minimum around noon. Fig. 7 
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Fig. I. 
and the appearance of polar blackouts during January, February, and March 1952. 


Aurorol frequency 


shows the curves obtained for October and 


31 


Day to day variation of the auroral and the magnetic activities, the intensity of the Nu-reflections, 


November 1951, the months of highest auroral 


activity, and for December 1951, a month 


with relatively low activity. The curves are 


typical. On the same graphs the hourly varia- 


tion of the magnetic activity as well as the 


intensity of the so called Nr-reflections and 


the number of polar blackouts has been plotted, 
as given by our Kiruna magnetograms and 


ionospheric soundings. yzmh is here defined 


i eee 


NO OBSERVATIONS 


as the maximum departure in y of the mag- 


netic z-component from undisturbed value 


during a certain hour, and y;mm is the 


maximum value of y;mh for the whole 


month. foNimm is the highest frequency for 
a certain hour during the whole month on 


1951 Moy June July Aug Sept Oct Nov Dec 1952 Jon. 


———— CURTAIN 


which echoes from the Nı-layer have been 
observed, and n is the number of times that 
Nı-echoes have been obtained during the 
month (n < 30). ZB is the total number of 


——'— STREAMER 
TOTAL 


Fig. 12. The monthly variation of the auroral activities 
during the time May 1951—March 1952. 
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Fig. 13. The reflection centres of auroral curtains during October 1951. 


blackouts observed during a certain hour for 
the whole month (2 B < 30). 

The Nr-layer is a special type of sporadic 
E-ionization often appearing in polar regions 
and always in connection with magnetic bays 
and supposed to be caused by the same 
ionizing agent as the aurora. The hf-trace of 
Ni looks quite similar to that of the ordinary 
E-layer except for the fact that the Nı-layer 
seems to be thicker and often shows compli- 
cated splitting near the penetration frequency 
(see fig. 8). For more detailed information 
about the Nr-reflections the reader may be 
referred to LINDQUIST (1950). 

The polar blackout is a certain type of “no 
echo” condition, which occurs in polar regions 
and shows no apparent connection with solar 
flares. The disappearance of reflections from 


all ionospheric layers during a blackout is 
believed to be due to an abnormally high 
absorption in a region below the 100 km 
level. This increased absorption may be caused 
by the impact of some ionizing agent in the 
ionosphere. The Nr-reflection appears before 
and after almost every blackout (Lrypqutst, 
1951). 

As can be seen from the graphs there exists 
a good correlation between the auroral fre- 
quency, the magnetic activity, and the in- 
tensity of the Nr-reflections. It is very inter- 
esting to note that the curve of the blackouts 
has the same form as the curves of the other 
phenomena but is delayed about 6 hours. 

The day to day variation of the recorded 
auroral activity is given on figs 9, 10, and ıı 
for the whole period of observation, May 
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Fig. 14. The reflection centres of auroral curtains during November 1951. 


1951— March 1952. On the same graphs we 
find the corresponding curves for the magnetic 
activity, the intensity of Nr-reflections, and 
the appearance of polar blackouts. yzmd is 
here the maximum departure from undis- 
turbed value of the magnetic z-component 
during the day, and foNimd is the highest 
frequency during the day on which echoes 
from the Ni-layer are observed. n is the 
number of times during the day that foN1 
has been obtained (n < 24). ZB is the total 
number of hours during the day on which 
polar blackouts occurred (2B < 24). There is 
a good correlation between all the curves and 
especially between those of the auroral activity 
and the appearance of polar blackouts. No 
significant difference exists between the two 
kinds of radar echoes, curtains and streamers. 


8— 203045 


The monthly distribution of the recorded 
auroral frequencies has a very marked maxi- 
mum in October and November 1951 and a 
second but smaller maximum in Februari 1952 
(see fig. 12). This is in good agreement with 
what is known about aurorae from visual 
observations during many years. 

During the time of observation the highest 
auroral frequency always appeared in the 
bearings northeast to northwest and very 
seldom in the south. This well agrees with 
the simple theory that the most probable 
direction of reflection is the direction of the 
magnetic north. To show typical cases we 
have plotted in polar coordinates on figs 13, 
14, 15, and 16 all the centres of reflections 
obtained during October and November 1951 
for auroral curtains and streamers respectively. 
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Distribution In direction 
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Fig. 15. The reflection centres of auroral streamers during October 1951. 


It is interesting to note that the optimal 
direction of reflection was constantly north to 
northeast during all the months except 
November when it changed to north-north- 
west. As for the distribution in range, the 
reflection centres of curtains have a maximum 
in the interval s00—600 km, and those of 
streamers, in the interval 550—650 km, as can 
also be seen on the same figures. 

From the direction and range measurements 
it is possible to calculate the height distribution 
of the reflection centres if one assumes the 
maximum reflected energy to come from 
those points where the radar beam is per- 
pendicular to the reflecting surface. These 
surfaces are supposed to be circular cylindrical 
surfaces parallel with the lines of magnetic 
force (streamers), or surfaces generated by 
the lines of force passing through the curves 
on the earth’s surface of constant magnetic 
horizontal intensity (curtains). Such calcula- 


tions have been made for all the reflection 
centres obtained within an angle of + 45° 
from the direction of the magnetic north (at 
Kiruna: magnetic declination 2°E, inclination 
77°). The results for October and November 
1951 are given in fig. 17. The height distri- 
butions have maxima of about 120 km and 
thus they fall in the region of maximum 
luminosity of visible aurorae. (CHAPMAN, 
BARTELS, 1940). 

The continuous radar recordings of the 
aurora give an indication of the comparatively 
rapid appearance and disappearance of the 
discharges, as well as the swift movements of 


the aurora, or rather the swift displacement of : 


the ionizing agent causing the aurora. This 
rapid variation is illustrated by the collection 
of pictures reproduced on fig. 18. The pictures 
show an aurora of a complex type with both 
curtains and streamers, lasting from 2228 
to 2245 I5°EMT, October 26th 1951. A 
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Fig. 16. The reflection centres of auroral streamers during November 1951. 


study of the continuous recordings gives times 
of duration of about 240 and 100 seconds for 
diffuse and discrete echoes respectively. The 
longest time observed during a representative 
period in October 1951 was 420 seconds for 
curtains and 240 seconds for streamers. Echo 
times shorter than 30 seconds occurred very 
seldom. 

From the amount of the reflected energy 
ASPINALL and HAWXKINS (1950) have calculated 
the typical electron density of a streamer to be 
approximately 6X 105 electrons per cc, sup- 
posing the streamer to be a circular cylinder 
(radius 1,000 m) enclosing a region of uniform 
electron density. An estimate on the same 
basis that we have made in some cases where 
the amount of the reflected energy could be 
determined with enough accuracy, gives den- 
sities of the order of 2x 10% electrons per cc 
as typical for streamers. This density value is 
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Fig. 17. Height distribution of reflection centres within 
an angle of + 45° from the direction of the magnetic 
north of Kiruna during October and November 1951. 
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Fig. 18. A series of records from an aurora appearing between 2228 and 2245 15° EMT, October 26, 1951. 


sufficient to explain the rapid fluctuation of where 


the echo pattern mentioned above. Solving Nısselestrons per cc 
the inhomogeneous Riccati equation, de- CR production of ions per cc and second 
scribing the simplest type of production and & = effective recombination coefficient 
recombination of electrons t = time, 
for the case of a sudden disappearance of the 
dN MAN ionization agent (q(t) = t < to; q (t) = 
dt 1 Fs | 


= 0, t > fg) gives as a solution 
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N | 
._1+aN(t—t) 


Taking N,, the initial value of the density, 
equal to 2 X 108 and choosing the value of the 
recombination coefficient found by RypBEck 
(1946) for the ionosphere at the height of 
100 km, = 1.1x10 we get an idea of 
the relaxation time t by putting «Nyt = 1. 
This gives us t & 45 seconds, which is in 
good agreement with the observations. 


PUY 


5. Concluding remarks 


Concurrently with the radar observations of 
the aurora visual observations were made. 
Because of unsuitable weather conditions the 
material is incomplete and a comparison 
between the results of the two observations is 
therefore very difficult. However, it turns out 
that during the months of highest auroral 
activity radar echoes were practically always 
seen on the scope when visual observations 
of the aurora were made, but often a simul- 
taneous radar and visual aurora did not agree 
completely with respect to direction and 
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character. There is little doubt, however, that 
the observed radio scatterings are due to 
reflections from auroral discharges. The distri- 
bution in direction and range of the echoes is 
the strongest evidence for the fact that the 
echoes could not come from the ionosphere 
or be any form of backscatter from ground 
via the ordinary or sporadic E-layer, the 
Nı-layer, or the F,-layer. 
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Abstract 


The 40-1, years of mean montly sea level pressure maps based on the Historical Map Series 
have been typed according to numerical indices that indicate to what extent each particular 
monthly map resembles a typical or normal winter or summer map on one hand, or a typical 
spring or fall map on the other hand. These indices are easily determined from the average 
pressure of certain key areas and are quite effective in separating the winter from the summer 
‚months and the spring months from the fall months. The persistence in these indices and their 
relation to the sunspot cycle has been investigated. 


I. Introduction 


In 1947 the writer (BRIER 1947) in reporting 
upon an investigation of the pressure data 
from the HistoricaL Dairy Map SERIES (1943 
—44) found indications that there had been 
a considerable decrease in the total mass of 
air over the Norther Hemisphere between 
1899 and 1939. The trend toward lower sur- 
face pressure was most pronounced in the 
higher latitudes and a further investigation 
(Brier 1948) disclosed that the trend existed 
in all months of the year. AHLMANN (1948), 
KINCER (1946) and others have reported that 
during this same period there seemed to be a 
general warming of the Northern Hemisphere, 
particularly in the higher latitudes and in the 
winter months. Because of general interest 
in the question of recent climatic fluctuations 
and the suggestion of Witetr (1949) that 
irregular solar activity is the most probable 
cause of climatic fluctuations, it was decided 
to examine the variations in the pressure 
data in relation to known solar effects, such as 


the seasonal change in the position of the sun, 
and to other possible effects that might be 
associated with an index of solar activity 


such as sunspots. 


2. The Seasonal Changes 


Regarding the seasonal changes in the distri- 
bution of pressure, there are two well-known 
and established facts. First, in the Northern 
Hemisphere winter there is an excess of air 
relative to the summer. This is shown by the 
annual cycle of monthly values of the average 
sea level pressure north of Latitude 20° N 
when plotted for each month from January, 
1899 to June, 1939, in Figure 1. These data 
are based on the series of MEAN MONTHLY 
Pressure Mars (1946 b) for the period. Since 
they are based on reduced sea-level pressures 
rather than on actual station pressures they do 
not represent the best possible estimate of the 
total mass of air. However, some computa- 
tions made by the author indicate that the 


FORTY-YEAR SEA LEVEL PRESSURE AND SUNSPOTS 


29 


203 


20 

” 

18 

17 

16 

15 

14 

3 4 

12 _02 où 04 8 06 L OT na 


2 os 19 | u 12 | 2 | 14 | is m iz 18 
= OBSERVED AVERAGE NH 
PRESSURE NORTH OF 20°N 


[ince at ER ty 


Fig. 1. Monthly values of the average sea-level pressure (in mb minus 1,000) north of latitude 20° N for the period 
1899 to 1939. Numbers along the abscissa are years (without the first two numbers) plotted in the middle of the year. 


effect of the artificial reduction to sea level 
is not a serious one for the purpose of this 
study and the qualitative picture given by 
Figure I is essentially correct. 


The second well-known fact is that the 
relative distribution of air over the land and 
water surfaces of the hemisphere is, in general, 
opposite in winter and summer. This diffe- 
rence, a seasonal effect, is illustrated in Figure 
2, which shows the normal difference in 
pressure between January and July. This map 
was obtained from the Normal Weather 
Maps (1946 a) which also were based on the 
40 ¥% years of sea-level pressure data. An addi- 
tional chart, Figure 3, shows the difference 
in pressure between October and April and 
indicates that for most of the hemisphere the 
seasonal march of pressure is not symmetric 


to January and July. 


3. Computation of Indices 


Since Figure 2 indicates that the average 
difference in pressure between January and 
July is of greater magnitude, and of opposite 
sign, in some regions than in others, it would 
seem reasonable to suppose that a few critical 
or “key” areas could be chosen from which 
numerical pressure indices could be computed 
that would show a good discrimination 
between individual months of January and 
July. The author has shown (BRIER, 1948) how 
monthly mean pressure maps might be typed 
or classified according to a few simple pressure 
indices chosen on the basis of areas showing 
maximum monthly change in surface pres- 
sure. An index, say W, might thus be com- 
puted for each of 486 months in the series, 
indicating whether a particular map was more 
like a winter map or a summer map. A 
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Fig. 2. Isopleths of difference in normal sea level pressure 
in mbs. between January and July for the Northern 
Hemisphere. 


statistical analysis of such an index might 
disclose interesting facts. Likewise, since Fi- 
gure 3 indicates differences between October 
and April on the average, it would appear 
possible to determine another index V which 
would indicate the extent to which a parti- 
cular monthly map resembled a spring map 
compared with a fall map. 

On the basis of Figure 2, four regions were 
selected where the average pressure was com- 
puted for each month of the series. The points 
selected are shown as crosses on the figure 
and the areas are labelled X,, X,, X, and X,. 

When the average pressure indices X, and X, 
were plotted against each other the points 
for January fell in the upper left side of the 
graph and the points for July fell in the lower 
right hand section. The points for the other 
months fell in between but with considerable 
overlap between adjacent months. The means 
for each month are shown in Figure 4, but 
the points for the individual months are not 
reproduced here. Using a graphical regression 
technique described elsewhere (BRIER, 1946) a 
series of arbitrary lines Y, indicated as curves in 
Figure 4 were drawn so that all the points had 
values of Y, ranging between o and 100. The 
spacing and curvature of these lines was deter- 
mined by the data for the individual months 
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Figure 3. The difference in normal sea level pressure 
between October and April for the Northern Hemi- 
sphere. 


which are not plotted here. The values of Y, 
were then read off for each individual month 
of the 40 year series and plotted against the 
index X, on a chart similar to that shown in 
Figure 5. Lines Z, ranging between o and 100 
were similarly constructed on this figure as 
shown. For each month, the value of Z, was 
read off the chart and a final index W was 
computed by the formula 
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Fig. 4. Scatter diagram showing the separation of the 
twelve monthly normal maps with respect to variables 
X, and X,. Months are numbered in order from 1 for 
January to 12 for December. The curves are indices of 
resemblance to the July normal chart. 
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Fig. 5. Scatter diagram showing the separation of the 
twelve monthly normal maps with respect to variables 
X, and Y,. The slanting lines are indices of resemblance 
to the January normal chart. Months are numbered in 
order from 1 for January to 12 for December. 
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Thus, the index W has been scaled to have a 
mean of zero and standard deviation of unity. 
The average values of W for each of the twelve 
months are shown in Table 2. An attempt was 
made to use variable X,, but it was found to 
contribute no additional discriminating infor- 
mation. 

For the spring-fall index, V, the regions 
chosen are shown in Figure 3. After a little 
experimentation it was found that the best 
discrimination was obtained when Y, was 
plotted against Y, where 


Ya = do 
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Fig. 6. Scatter diagram showing the separation of the 
twelve monthly normal maps with respect to variables 
Y; and Y,. The slanting lines are indices of resemblance 
to the April normal chart. Months are numbered in 
order from 1 for Januarv to 12 for December. 


After Y, was plotted against Y,, the lines 
Z, were drawn as shown in Figure 6. When 
the individual monthly values of Z, were 
determined from their position on Figure 6 
and tabulated they were found to have a 
slight correlation with index W. In order to 
have an index independent of W the linear 
regression of Z, on W was determined and a 
final index V was computed by the formula 


where 


== 486 
2, = 2 Z,,/486 


1 


and b is the regression of Z, on W and o, is 
the standard deviation of the residuals from 
the regression line. Thus, index V has zero 
mean and unit standard deviation and is 


Table 1. Average values of index W for each month. 


Month | Jan. | Feb. | Mar. | Apr. | May | June 


| 
July | Aug. | Sept. | Oct. | Nor 
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Fig. 7. Scatter diagram showing how the W and V in- 
dices separate the January maps from the July maps and 
the April maps from the October maps. Months are 
numbered in order from 1 for January to 12 for December. 


uncorrelated with W since it is a residual 
from the regression line of Z, on W. 

The effectiveness of these two indices in 
separating the winter from the summer maps 
and the spring maps from the fall maps is 
indicated by Figure 7. On this figure the mean 
values of W and V are shown for each of the 
twelve months of the year but the individual 
monthly values are plotted only for January, 
April, July and October. If the remaining 
individual months were shown they would 
indicate considerable scatter and overlap bet- 
ween adjacent calendar months. 


4. Statistical Analysis of the Indices 


A number of correlations were computed 
between W and V and other factors. The 


Table 2. Linear Trends in W and V. 


Correlation Regression 
Month Coefficient Coefficient 
w V RS nz 
January .... 0.45 0.35 0.0068 | — 0.029 
February ... 0.59 0.44 0.0090 | — 0.032 
Nfarchie ne 0.60 0.23 0.0058 | — 0.016 
END Talley) Are 0.64 0.07} — 0.0079 | — 0.004 
Mayr: 0.55 0.21 0.0071 | — 0.015 
RER RE — 0.47 | — 0.32] — 0.0072 | — 0.019 
TOP eee 0.49 0.18 0.0119 | — 0.009 
ENTS TS Beever 0.2 0.30 0.0046 | — 0.017 
September .. 0.63 0.05| — 0.0065 | — 0.002 
OCtODer iia 0.44 0.08 0.0051 | — 0.005 
November .. 0.56 0.39 0.0072 | — 0.033 
December .. 0.45 0.33 0.0084 | — 0.023 
Annual: 0.75 0.36 0.0732 | — 0.070 
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Fig. 8. Auto-correlations of index W for one and two 
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Fig. 9. Auto-correlations of index V for one and two 
month lags. 


persistence in these indices was investigated by 
computing the one and two month lag auto- 
correlation for each month of the year. The 
results of these computations are shown 
graphically in Figures 8 and 9. The one month 
lag correlations were highest for W with 
values greater than 0.50 for April with May 
and for May with June. The average of the 
one month lag correlations was 0.35 for W 
and 0.27 for V. Ordinarily, these would not 
be considered as very useful for forecasting 
purposes. 

The trends in the indices during the period 
were investigated by computing the regression 
on time for each index and for each month. 
The results, shown in Table 2, indicate that 
during the 40 years there were pronounced 
trends in both of these indices and that the 
sign of the trend was the same in all seasons. 
Since positive values of W means more 
winter-like indices, the negative correlations 
mean that on the average the monthly mean 
patterns were more like normal winter pat- 
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terns in the earlier part of the period and 
more like the normal summer pressure pat- 
terns toward the end of the period. The most 
pronounced trend was during the month of 
April and Figure ro shows a plot of the indi- 
vidual values of W for the 41 months of April 
and the regression line on time. This tendency 
for the maps to show more summer-like 
patterns during the 40 year period appears 
to be consistent with the findings of AHLMANN 
(1948), KINCER (1946) and other investigators 
reporting a general warming up of the hemis- 
phere during this period. The negative corre- 
lations for index V are not so easy to interpret 
but indicate a tendency for the monthly mean 
pressure patterns to become more like normal 
fall patterns during the 40 years. 

Index W was determined on the basis of 
areas showing the greatest change in pressure 
from January to July. Since the basic cause 
of this seasonal change is undoubtedly radia- 
tional in character, due to the changing posi- 
tion of the sun relative to the earth, it would 
seem reasonable to hypothesize that real 
changes in solar activity might affect the same 
regions. Therefore, the index W was corre- 
lated with sunspots (Munro 1948) the easiest 
available measure of solar activity; the re- 
sults for each of the twelve months are shown 
by the solid line in Figure 11. Correlations 
between V and sunspots were also computed 
although there was no a priori reason to ex- 
pect any relationship. The solid line in Figure 
12 shows a graph of these correlations. 

For index W, all twelve of the correlations 
turned out negative, ranging in value from 
—0.07 to — 0.44. Using the standard test of 
significance, two of these would be considered 
significant at the 1 % level and another one 
at the 5 % level. The algebraic average of the 
twelve monthly correlations is — 0.26 and on 
the basis of the ¢ test this is a significant de- 
parture from zero. For index V the situation 
is quite different. None of the twelve monthly 
correlations is significant at the 5% level, 
and the average of the group does not depart 
significantly from zero. 

Since the analysis indicated a considerable 
trend in W and V during the 40-% year 

eriod, the question might be raised as to 
whether the observed negative correlations 
between W and sunspot numbers might only 
be the result of a trend in the sunspot numbers 
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Fig. 10. Time series plot of index W for the month of 
April. 


opposite the trend in W. For example, the 
trend in W during the period might have no 
physical reality, resulting from some bias in the 
data introduced by the analysis of the original 
synoptic maps which were based on less 
complete data in the carlier years. Since sun- 
spots showed a slight upward trend during 
the period, a “spurious” correlation between 
them and W would be found even though no 
physical relationship existed between the two. 
For this reason the correlations between suns- 
pots and the indices were recomputed after 
removing the linear trends in W, V and suns- 
pots. The results are shown graphically by 
the dashed lines in Figures 11 and 12. 

The removal of the trends in the series re- 
duced the correlations between W and sunspots 
only one of them (October) remaining signi- 
ficant at the 5 % level. However, 10 of the 
12 correlations are of the same sign as before 
and the average for the group of — 0.10 de- 
parts significantly from zero. The fact that 
the correlations of W and sunspots are nega- 
tive after the trends are removed and that the 
trend in sunspots is opposite to the trend 
in W strongly suggest that the correlations 
found between sunspots and W presented 
by the solid line in Figure 11 are not spurious. 
Furthermore, it is the author’s opinion that 
the trends in index W during the period were 
real and cannot be explained by asserting 
that insufficient data in the carlier years 
biased the results in a considerable way. The 
correlations of sunspots with V show essen- 
tially the same thing as before the removal 
of trends, the group average of 0.05 not being 
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Fig. 11. Correlations of index HW” with sunspots, before 
and after removal of the trend during the period. 


significant although the January correlation 
of 0.34 is significant at the 5 % level. 

The validity of some of the results reported 
here might be open to some question due to 
the fact that some of the points used in 
determining the indices W and V are located 
in regions of less reliable data and weather 
analysis. For example, variable X, includes a 
point located in the mountains of Mexico and 
another point located in the Pacific off Lower 
California. This question was investigated by 
substituting for variable X, the mean monthly 
pressure data for Corpus Christi, Phoenix, 
and San Diego. The correlations of these 
stations with variable X, averaged 0.92. 
Likewise, the correlations between variable 
X, and a new variable X,’ which eliminated 
questionable points located in China and Tibet 
ranged between 0.87 and 0.94. These changes 
were too small to have an appreciable effect 
on the W and V indices and their relationships 
to the sunspot series, and indicate that the 
significant correlations reported here cannot 
be simply explained by the assertion they 
were produced by faulty data or map anal- 
ysis. 


5. Conclusions 


A method of map classification has been 
developed which can be used to type monthly 
mean sea-level pressure maps according to 
the degree to which they resemble a normal 
winter or summer circulation pattern on one 
hand and a normal spring or fall map on the 
other. Index W based on the average pressure 
of points chosen from three regions measures 
the resemblance to the January normal map 
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Fig. 12. Correlations of index F with sunspots, before 
and after removal of the trend during the period. 


and index W is a measure of the resemblance 
to the normal spring map. 

Indices W and V show considerable per- 
sistence from one to two months with auto- 
correlations ranging from — 0.03 to 0.58 for 
one month lags and from — 0.09 to 0.51 for 
two month lags. Some of these correlations 
appear to be large and consistent enough to 
suggest a possible usefulness in long range 
forecasting. 

Both indices show rather consistent negative 
trends during the period, indicating that the 
monthly mean patterns were more like sum- 
mer and fall normal pressure patterns toward 
the end of the aed and more like winter 
and spring during the earlier part of the period. 
These results seem to be consistent with the 
findings of others investigating the climatic 
flucatuations of the past half-century. 

For each month the correlation of index W 
with sunspots was negative, indicating that 
periods of high sunspots tended to be associated 
with summer type patterns, at least so far as 
these are determined by the distribution of 
pressure in regions chosen for this study. 
Since most of the correlations remained nega- 
tive when the trend was removed from both 
the sunspots and W series, there appears to be 
the definite suggestion that part of the varia- 
tion in the monthly mean sea level pressure 
patterns during the period 1899 to 1939 was 
Tiers by variable solar activity. The corre- 
ations are not large enough to indicate that 
a very substantial portion of variation was 
solar in origin and they are not high 
to be of much forecasting value, even if they 
are real. However, the results appear to give 
some slight supporting evidence to the con- 
clusion oS | y WILLETT (1949) that irre- 
gular solar activity is sible for some of 
the changes in the circulation of the atmos- 
phere. 
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